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I* INTRODUCTION 



Although many thin-wire computer programs have been 
developed for the purpose of analyzing antennas and 
scatterers, few of these programs have been directed toward 
the student of electro-magnetic theory. The majority of the 
programs are directed to the engineer or advanced student 
for the purpose of analyzing designed structures or 
verifying experimental data. 

The purpose of the study is to develop a computer 
program by modifying an existing computer code which can be 
utilized as an educational method to develop insight into 
radiating structures by the beginning student of 
electro-magnetic theory. 

The modified Ohio State University Antennas-Scatterers 
Analysis Program (OSUMOD or ASAP) is directed toward the 
beginning student who does not yet have the expertise 
necessary to manipulate the input data for proper execution 
of the larger more comprehensive analysis program. Even 
though ASAP is small in core requirements and is fast in run 
time, it is capable of analyzing structures to assist the 
engineer with design problems. 

Since the resulting program, ASAP, is primarily directed 
toward students, the program has been limited to structures 
which contain less than 50 monopoles (segments) , no longer 
than one-fourth of a wavelength, and which have less than 51 
nodes (intersections and endpoints) . If a ground plane, 
either perfect or finite is present; the stated limits above 
are halved due to the generation of an image structure. 



7 



II* ORIGINAL PROGRAM 



A. THEORY 

Reference 1 presents the electro-magnetic theory for the 
analysis of antennas and scatterers in an isotropic, linear, 
and homogeneous ambient medium. The analysis is performed in 
the freguency domain with an excitation caused by either a 
generator or an incident wave. 

In the analysis, a piecewise-sinusoidal expansion is 
used for the current distribution. The matrix equation 
Z I = V is generated by enforcing reaction tests with a set 
of sinusoidal dipoles located in the interior region of the 
wire. Since the current distribution has the same form as 
the expansion mode, this formulation is known as the 
"sinusoidal reaction technique". 

B. COMPUTER PROGRAM 

Reference 2 presents the computer program corresponding 
to the theory presented in Ref. 1. 

1. Input Format 

In the program, the input data must specify the 
frequency, wire radius, wire conductivity, the parameters of 
the exterior medium, coordinates of the points to describe 
the shape and size of the wire configuration, a list of the 
wire segments, and the indicators for the various outputs. 
Table 1 is the input data necessary to analyze a half-wave 
dipole. 

2. Output Format 

In the original form, the only outputs which could 
be requested by the input data stream are the following: 
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a. 


Antenna Problems 




(1) 


Current Distribution on the Structure 




(2) 


Input Impedance. 




(3) 


Radiation Efficiency. 




(4) 


Near-Zone Field. 




(5) 


Far-Zone Field. 


b. 


Backscatter ing Problems 




(1) 


Absorption Cross Section. 




(2) 


Scattering Cross Section. 

• 




(3) 


Extinction Cross Section. 




(4) 


Complex Elements of the Polarization 


Scattering 


Matrix 


c. 


Bistatic Scattering Problems 




Echo 


Area. 



Table 2 is an example of the output data available for 
data of table 1. 



3. LIMITATION 

Although the program can analyze a structure with up to 
50 segments, 55 points and 60 dipoles modes; it can not 
analyze a structure in the presence of a finite ground 
plane. 
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AN EXAMPLE OF THE OUTPUT DATA FOR THE ORIGINAL PROGRAM 



Ill* MODIFIED COMPUTER PROGRAM 



A. Input Format 

As illustrated in table 1 the format for the input data 
cards is not self explanatory. This format can be determined 
by referring to the FORMAT statements of the program of Ref. 
2. Since the modified program is directed toward the 
student, the input data format was changed to allow free 
format. Reference 2 was written in a form which permitted 
modifications to allow flexibility in specifying input data 
for the analysis program. Appendix B, titled "User's 
Manual", discusses the input data cards necessary for proper 
execution of an analysis problem. Appendix B is self- 
contained and may be used independently of the remainder of 
this document. 

B. Output Format 

In the original computer program, the absence of labels 
encumbered the output data and lessened the usefulness of 
the program. To improve the usefulness of the modified 
version, detailed labels were added to the output data. As 
with the input data. Ref. 2 was written in a form which 
enabled modification to allow more specific output data for 
the analyzed problem. With the addition of the polar 
plotting package, the far-zone electric field intensity 
polar radiation and reradiation patterns can be plotted. A 
sample problem can be found on page 120 in Appendix B, 

User's Manual. 

C. Finite Ground 

To enable the student or the engineer to have an 
improved analysis program, the finite ground effects were 
added to ASAP. The theory corresponding to the ground 
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effects, which utilize Fresnel reflection coefficients, is 
discussed in Appendix A, titled "System Manual". Also 
discussed in Appendix A is the modified computer program and 
the corresponding theory. The electro- magnetic theory was 
developed in Refs. 1, 2, and 3; and it is restated with its 
corresponding computer code to assist in the understanding 
of the methods applied. Appendix A is self-contained and 
may be used independently of the remainder of this document. 



IV* CONCLUSION 

The addition of ground effect technigues to the original 
program did not alter the accuracy or the computational 
capabilities of the program. The ground effect technigues 
utilized the results of the original program and modified 
these results to account for the effects of the presence of 
the finite ground. 

To verify the numerical results of ASAP, the input 
impedances of both a horizontal and a vertical dipole were 
compared to the solutions of the exact form of the 
Sommerfield 1 s eguation. As can be seen in table 3 the 
finite ground treatment of ASAP agrees favorably with 
Sommerf ield 1 s solutions. The ASAP finite ground results are 
also in excellent agreement with the previous computer 
solutions of Refs. 4 and 5. 



V* RECOMMENDATIONS 



Although the program is a general analysis tool for 
students, several future modifications will enhance the 
program as a design tool for engineers. These items include: 
varying the wire radius on the structure; incorporation of 
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non-radiat ing elements such as transmission lines; varying 
the wire insulation radius, conductivity, and dielectric 
constant; and a geometry generation package such as dipole 
array or helix. One manor change that would both improve the 
speed and reduce the core requirement is that of symmetry. 

No attempt was made to utilize the symmetry in the 
admittance matrix when the ground plane is present. If 
symmetry were applied, the structure size limit with the 
ground plane present would be approximately that of the 
structure without the ground plane. 



VERTICAL 0 1 POL E 
FREQUENCY 3MHZ 
LENGTH .5 WAVELENGTH 

RA01US .005 METERS 

* DIELECTRIC CONSTANT (RELATIVE) 
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CONDUCT LV1TY HE 1 GHT/ WAVEL E NGTH 



.1 



.00 1 



.00001 



.25 
.30 
.35 
.45 
.25 
.30 
.35 
.45 
. 25 
.30 
.35 
.45 



108 

92 

82 

74 
98 
86 
80 

75 
97 
86 
80 
75 



ASAP 

. 43* J 
► 43* J 
. 5 9* J 
.0 1* J 
.74*J 
. 15*J 
. 1 5* J 
.42*J 
.94* J 
.1 2*J 
. 82* J 
.93* J 



63.22 
37.73 
34.76 
40.11 
63.52 
39.02 
38.08 
41.79 
64 • 86 
40.85 
38.89 
41.55 



EXACT* 



113, 


.7 *J 


56. 


22 


93, 


,97* J 


32. 


62 


83, 


, 56* J 


29. 


74 


74, 


,87* J 


35. 


34 


99, 


,86* J 


46. 


43 


87, 


.40* J 


34. 


,04 


81 , 


. 1 1 * J 


33. 


,16 


76, 


, 32* J 


37. 


,03 


96, 


.68* J 


49. 


, 1 1 


8 7 , 


. 37 * J 


35, 


,94 


81, 


,79* J 


33. 


,99 


76.83* J 


36, 


.78 



HORIZONTAL DIPOLE 
FREQUENCY 3 MHZ 
LENGTH .5 WAVELENGTH 

RADIUS .001 METERS 

OIELECTRIC CONSTANT (RELATIVE) 
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CONDUCTIVITY 


HEIGHT/WAVELENGTH 


ASAP 




EXACT* 


. 1 


.5 


7 2 . 9 5* J 


28.53 


78 . 25 * J 


23.91 




. 3 


103. 70*J 


58.27 


117.5 *J 


52.05 




. l 


30. 16*J 


69.54 


33.20* J 


73.65 


.001 


.5 


74 • 2 3* J 


35.84 


80.66* J 


31.84 




. 3 


94. 16*J 


51 . 32 


105.3 *J 


46.01 




. l 


54.52* J 


62.89 


61.31* J 


64.33 


• 00001 


. 5 


76. 13*J 


35.99 


82 . 79* J 


31.73 




. 3 


9 1.8 5* J 


53.21 


103.0 *J 


48. 56 




. 1 


54. 3 0* J 


56.81 


61 .03* J 


57.27 



* COURTESY OF LAWRENCE LIVERMORE LABORATORY 



TABLE 3 



APPENDIX A 
SYSTEM MANUAL 
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SYSTEM MANUAL 



INTBODUCTION : The Antennas-scatt erers Analysis Program 

(ASAP) for thin wire structures in a homogenous conducting 
medium performs a frequency domain analysis of antennas and 
scatterers. The program is applicable in the presence of a 
ground either perfect or finite. This appendix will describe 
the computer program which accomplishes this. Although the 
program was written for the IBM 360 computer system it can 
be executed on another system with minor modifications. 

A piecewise-sinusoidal expansion is used for the current 
distribution. The matrix equation ZI = V is generated by 
enforcing reaction tests with a set of sinusoidal dipoles 
located in the interior region of the wire. Since the test 
dipoles have the same current distribution as the expansion 
modes, this may be regarded as an application of Galerkin's 
method. Rumsey's reaction concept was most helpful in this 
development, and therefore the formualtion is known as the 
"sinusoidal reaction technique". 

The main routine and each subroutine is discussed 
separately in this appendix. The writeups for the 
subroutines are arranged alphabetically by subroutine name 
after the main program. Each of the discussions includes the 
purpose of the subroutine, brief description, and a listing. 
After the subroutine writeups is a table of the more common 
symbols used in this program. 

The input data and program limits are discussed in 
detail in the next appendix titled "USERS MANUAL". 
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GROUND EFFECTS: In the modified antenna analysis 

computer program finite and infinite ground effects were 
added by using the reflection coefficient technique. The 
method in which tnis technique was used required the 
generation of an image structure. In this section the 
reflection technique will be discussed in detail. 

In order to apply ground effects to the electric field, 
the field for the image structure was first calculated as if 
a ground were not present. Then, the field was decomposed 
into parallel and perpendicular components. (A parallel 
component is the component which is parallel to the plane of 
incidence. A perpendicular component is one which is 
perpendicular to this plane. The plane of incidence is the 
plane containing the normal to the reflecting surface and 
the incident ray.) 



Z 1 
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Consider an image monopole with the electric field in 
the 1 direction. The ray, S, is a vector which is 
perpendicular to 1 and passes thru the point of interest. To 
apply reflection technique, the plane of incident must be 
found. It is advantageous to define a new coordinate system 
(a»/3rZ) where and are parallel to the xy plane with a in 
the plane of incident and /3 perpendicular. 



If the direction cosines (cosx, cosy, and cosy) are 
known, it can be shown that the components of the field in 
the a ) 3 (xy) plane have the following relationship: 



II 

E 












— 




COS 


4> 


sin <t> 




E 

X 


E 1 ‘ 




sin 


4> 


-cos $ 




E 

_ y _ 



where <) = arctan (cosy/cosx) . 



Now the reflection coefficients for the interface can be 
applied as: 



II (R) n J 

E = R E 

I I 



-L ( R) L 

E = R E 



where R and R will be defined later in this section. 

I I 1 

Applying the matrix equation above yeilds: 
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(R) 








1 1 (R) 


E 




cos <|> sin <|> 




E 


X 










(R) 








_L(R) 


E 




sin $ -cos <(> 




E 


y 












(the square matrix is unique, in that, the inverse is equal 
to the original matrix) . Since the image direction is 
opposite to the original monopole, that is. 



(1 x z) . . = - (1 x z) . , 

original image 

the z component of the field, which is in the plane of 
incident, is given by: 



E 



(H) 



z 




E . 

z 



From electro-magnetic theory the reflection coefficients 
for the fields in medium (1) at the interface with another 
medium (2) are defined as: 

for parallel 

cos 0 - \/ S' - sin 2 0 

R = ; 

H cos 0 + Y&* - sin 2 0 

and for perpendicualr 

- &' cos 0 + \/ &' - sin 2 0 

R — - / ■ ~~ 

V &• cos 0 + v S' - sin 2 0 

where 0 is the angle of incident as measured from the normal 
to the interface and 



&• = (S + o/joj)/(& + a /jo>) 

2 2 1 1 

where the subscripts correspond to the mediums above. 
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To determine the relationship between R , R and R , R 

I I J_ V H 

a perfect ground (£ =0, o = oo ) was investigated. 

r 

limit R = - 1 
H 

limit R = +1 
V 

But, for a perfect ground the contributions to the field 
from the image monopole would be equal to the field of the 
original monopole but opposite in sign due to the chosen 
reference direction, 

I I (R) II II 

E = R E = - E 

I I 



±(H) 



= R 



J_ -L 

E = - E 



therefore 



R = R 

I I H 

R = - R 



In summary, the contribution to the electric field of a 
monopole over a ground plane at a given point is given by: 

(R) (R) <R) (R) 

E = E cosx + E cosy + E cosz 
x y z 

where 



(R) 



E = R E cosx + (R - R ) E cosx cos* $ 

* J_ II _L 



+ (R - R ) E cosy sin $ cos $ 

I I _L 
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E = R E cosy - (R - R ) E cosy cos 2 $ 

y ii 1 1 J_ 



+ (R - R ) E cosx sin <t* cos $ 

I I J_ 



E = - R E cosz 

z I I 



where E is the field without the ground plane present and 



- &' cos 9 + v/lT* - sin 2 9 

R = - - — — — 

|| &' cos 9 + v - sin 2 9 

cos 9 - v/&' - sin 2 9 

R = ~ > t 

JL cos 9 + y &' - sin 2 9 



&• = & - j (a/& oj) 

r o 
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MAIN 

PURPOSE: to control the input, output, and the flow of 

calculations. 

METHOD: The main program controls the flow of the 

required calculations by calling only a few subroutines. 
These subroutines in turn call other subroutines which 
actually do the required calculations. The order of the 
calling sequence is diagramed after the listing for the main 
program. 

The DIMENSION statements at the beginning of the main 
routine provides the required storage for a wire structure 
with up to 50 segments, 60 nodes and 60 dipoles without the 
presence of a ground plane. If a ground plane is present 
one-half of the reserved storage is required for the image, 
therefore a wire structure with up to 25 segments and 30 
nodes can be analyzed. 

NM denotes the actual number of monopoles (segments) , 

INM is the corresponding dimension, and the dimension for 
CG,. VG, and ZLD is twice INM. The second subscript for MD 
always has a dimension of 4 to correspond to the number of 
segments meeting at a given node. 

N denotes the number of simultaneous linear equations 
and ICJ is the corresponding dimension. The dimension for C 
is (ICJ * ICJ + ICJ)/2. 

In the statements above statement 4, the initial 
conditions and defaults are established. After calling 
subroutine READ to determine the input parameters, the IF 
statements output the parameters to be used for the 
calculations. In the DO LOOP ending at statement 7, the 
the input data of the structure geometry is stored in order 



23 



to recall if the structure is to be moved for ground plane 
calculations. 

After the image structure is generated and structure 
location is moved, subroutine SORT is called to determine 
the dipole modes. Prior to calling SG ANT , the load and 
generator information is established. 

Subroutine SGANT is then called to calculate the 
elements of the impedance matrix. If FEEDS or GENERATORS 
are specified by the input data stream, subroutine GANT1 is 
called to solve for the current distribution due to these 
forcing functions. 

In the DO LOOP ending with statement 29, subroutine 
GNFLD is called to calculate the near-zone field for the 
current distribution of the subroutine GANT1. 

The subroutine GFFLD is called for the far-zone field of 
the current distribution of the subroutine GANT1 in the DO 
LOOP ending at statement 35. The subroutine GFFLD is called 
again in DO LOOPs ending at statements 42 and 51, if 
bistatic and backscattering calculations are requested by 
the input data stream. 

CALLS TO: GANT1 

GFFLD 
GNFLD 
POLPRT 
READ 
SGANT 
SORT 
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D I ME NS I UN X(60). Y(60), Z<60>, XG(60), YG(60), ZGI60) 

D I Mb NS 1 UN (11601, 12(60), 13(60) f J A i 63 ) t J0(6O), KF LAG I 30 ) 

0 I M f N S I UN C PH( ( 500 » . C I HE T ( 500 ) 

UIMINSIUN i) A T Y I ( 36 0 1 » DATy 2<366), 0ATY3I360). DATY*(360) 
niMI NSl (iN 0( SO) * ( A ( 5 0 ) f ( B ( 50 ) t ISC(50), M0(50,*|, No! 50 ) 

DIMlNSIUN L ZD < 60 ) » KGtN ( 60 ) 

OlHtNSION XNP(5J)» YNP ( 50 ) * ZNP150) 

COMPLEX c< I B 30 ) 

COMPLEX COAT I (500 ), COAT 2 ( 500), COAT 31 503) .COAT* (500) 

COMPLEX CJ<60)»EP(60) ,EPP(60),ET( 60 ) « E T T (60 ) 

COMPLEX CGO( 50) * SGD( 60) .CGI 100) ,VG( i 0 0 ) , ZLO(lOO) 

COMPLEX VOLT ( 60) »ZLL0(60) 

COMPLEX EPPS. EPT$.ETPS,ETTS ,EX,EY.EZ 
COMPLEX EP2,EP3,Er*,ERR,ETA,GAMyYlltZll , ZS 
DATA PI , TP/3. 1*159.6. 20318/ 

DATA EO,UO/8.05*E-12,1 .25666-6/ 



'1 



-1 

1 



NGEN 
1GRD 
LOAD 
8 M * ■ 
ICARD 
AM = -I 
IFLAG * 0 
VOLT(l) - 
HOT • 0. 
NM * 0 
NP * 0 
MSG - 0 
SIG2 - -1 
F02 - -1. 
S1G3 * -I 
ER3 



( 1 .,0. I 



T 03 * 
CMM * 
ER2 - 
FMC * 
1NH = 
(CJ - 
WRITE 



0. 

50. 

) . 

300. 

50 

60 

(6,7*) 



00 2 1 7 1 , 30 
2 KFLAGI 1 ) 



-1 



00 3 J«l , 1NM 

iscuT - o 

VG(J) * t.0.,.0) 
ZLOM) - 1. 6, .3) 



2 

3 

* 

5 

6 
7 
0 
9 

10 

11 

12 

13 

I* 

15 

16 
17 
10 

19 

20 
21 
22 
23 
2 * 

25 

26 
27 
20 

29 

30 

31 

32 

33 
3* 

35 

36 

37 
30 
39 
*0 
A 1 
*2 
*3 
** 
*5 
*6 
*7 
*0 



J J 



VG ( J J ) 
ZLD( J J ) 



J* 1NM 



( . 0, .0 » 

* l.0,.0) 



* NFFP ■ 0 
N 8 I P * 0 
NPAP a 0 
AFFP - 1000. 

AFFT - 1000. 

ARIP - 1000. 

ABIT c 1000. 

ABAP * 1000. 

A0AT - 1000. 

STFP » l . 

KNM = 0 

CALL RE AO ( 1A, 10, (0 ISC .ICARD, IGA IN. 1GRD, (NEAR, l NT , (SCAT. IWR, IFLAG. 
|KFLAG,KG6 N,l HAD, L Z 0 , MSG, NB A P, NB I P.NFFP. NGE N ,NM. N P , A0 AP. A0AT.AFFP.A 
2FFT,AB(P.A>3I T 1 AM,0M,CMM.fcR2,eR3,ER*.FMC I HGT ,PHAF»PHAI.kH1F,PH( 1,*H 
3SF.PHSI ,)HAF,THA(,TH1F,TH11 . T HS F , T HS 1 , S I G2 , 5 1 G3 , S 1 G* , T 0 2 , T03 , VOLT , 
*X, XNP,Y, YNP. Z,ZILD,ZNP, STEP) 

WRITE (6,56) 

IF (MSG.LT-1 I GO TO 5 

IF (MSG. tO. I) WRITE (6.70) KFLAGI30) 

IF (I FL Ao.fcO.* ) GO TO l 

5 IF (IFLAG.FQ.5) STOP 

IF ( AM. L T .0 » WRITE (6,127) 

I F (AM. L T .0 ) GO TO 6 

IF ( (NM.GT .0) .AND. (NP.GT.O) ) GO TO 7 
WRITE (6,116) 

6 IF ( lFLAG.EO. 1 ) GO TD 1 
MSG - 2 

GO TO * 

7 WRITE (6,11*) 

WRI TE (6,113) 

WR 1 T F (6.112) 

(F ( KF L AG ( I ) • EO. 1 ) WRITE 16,03) F MC 
IF (KFL AG( 2) . EQ. I I WRITE (6,0*) AM 
(KFLAGO).EO.L) WRITE (6,05) CMM 



I F 
(F 
IF 
1 F 
IF 
IF 
IF 
IF 
IF 
l F 
IF 



b,6DI 

(KFLAGI 20I.NE . ( I WRITE (6,07) 
( KF L AG i *1 .EO. 1 J W® 1 TE (6 ,06) 



(KFL AG *) . EO. 1 I 
(KFL AC( 5) .t 0. I ] 
(KFL AG( 6) .EO. 
(KFLAt,( M .( 0. 
(KFL AG( O) .Nt . 
(Kf LAC.(9) .E 0. 
(KFLAG( 10). tO. 



WRITE 
WRI TF 
WR ( Tt 
WRITE 
WRITE 
WRITE 
> WR ( Tfc 



( 6 , 00 ) 
(6,09) 
(6,93) 
<6,91 ) 
(6 ,92 > 
(6,93) 
6,9*1 



0 M 

S IG2 
E R 2 
TD2 

S ( G 3 
ER3 



(KFL AG( II ) .EO. 1 ) WRITE (6,95) T03 



*9 

50 

51 

1 3 
5* 

55 

56 

57 
50 

59 

60 
61 
62 
63 
6 * 

65 

66 
67 
60 

69 

70 

71 

72 

73 
7* 

75 

76 

77 
70 
79 
00 
01 
02 
0^ 
0* 

05 

06 
07 
0B 
P9 
9 0 

91 

92 

93 
9* 

95 

96 



25 



i I 2$ ! f E<L 1) ) WRITE 



) 6, 1 2 3 I 

16.125) 

( 6 , 12*) 

16.126) 



ER9, SIG9 
HGT 



IF (KFL Abl 26 ) .ME . 1 ) WRITE , 

IF I I 1GP0.GT . ) ) . AND. IKFLAGf _ 

IF IUGPD.FQ. ) ) . ANO. (KFL AGI 25) .EQ.IM WRITE 
IF I I It PU .GT . I ) . AND. I KFL AG ( 25 ) . E U. I ) | WRITE 
IF lUGRD.GT .0) . AND. I KFL AGI 25) .EO. I » 1 WRITE 
IF (KFLAGI2I I.EQ.I ) WRITE 16,121) INT 
WRITE 16*1111 

IF I KFL ACI 12 ) .E3 . 1 ) WRITE (6,96) I I , I A I I I , X I I A( I)) , Y I 1 A I 1)1 , Z 1 I A (I 
))), 1811), X|I8(I)),Y(|BII)), Zimi)), 1*1, NM) 

WRITE 16,11)) 

IF IKFL A(,) 29 ) .GT .0 ) WRITE 16,1)9) 1 1 ZO I I I , Z ll 0 1 I ) « 1 ■ 1 , 10 AO ) 

IF IKFI AGI 19 | .GT .0) WRITE (6,118) ILZOI I ) • ZLLUt 1 ) , I* ltLOAOl 

WRITE (6.11)) 

) F JKFLAtl 23) .GT .0 ) WRITE 16,120) I KGENl I I , VOL Til) » 1 -1 » NGE N) 

IF (KFL AGI 13) .GT .0) WRITE (6,97) I KGENI 1 ) , V 01 T 1 1 ) , I • 1 1 NGE N ) 

WRITE I 6, I I ) ) 

WRITE 16,11**) 

WRITE 16,93) 

WHITT (6,112) 

IF (KFLAGI22 .NE.l) WP1TE 
IKFL Ag( I 5 ) .fc 0. I I WRITE 
lKFLAOIiCl.EQ.il WRITE 
IkFLAGI 1 n .EU.I ) WRI IE 
(KFLAGI 10). E3. I ) WRITE 
IKFL AGI 19 ) .EQ. I ) WRITE 



(6,110) 

( 6.99 ) 

(6.100) PH A | t PHAF , THAI ,THAF,$TEP 

16.101) PHI I ,PHIF ,THI 1, THIFtSTEP 

16.102) PHSi .PHSF tTHSI ,THSF,STEP 

16.103) (XNPi I I.YNPI 1 | , ZNPll I , l-I i I NEAR) 
AFF P 
AFF T 
A 0 AP 
ABAT 
AB IP 



I Oi iu 

I AFF P.c T .500. ) WRITE I (5 , I 6 5 ) 

IF (AFF T. IT. 500. ) WRITE (6, 109) 

(F ( ASAP. IT. 500. ) WRITE (6,109) 

IF I A8AT.I T.500. ) WRITE (6,108) 

IF ( A8 I P. L T. 500. ) WRITE (6,107) 

IF ( A8) T.L T .500. ) WRITE (6,106) A8IT 
IF ( ) 18 ISC.GT.O) .ANO. I ISCAT.LT.O) ) WRITE (6,73) 
F HZ • FMCM.E6 
OMEGA * T P ♦ F HZ 

IF IS1G2.ET.0.) EP2-E«2*E0*CMPLX( I .,“T02 ) 

IF (T02.LT.0.) E P 2 ■ C MPl x I ER 2 • E 0 . - S 1 G2 / OME GA ) 
IF (SIG3.LT. 0.) E P3»ER3*E0*CMPLX ( I..-T03) 

IF (T03.LT. 0.) E P 3 ■ CMPLX( ER3*EO,-SIG3/OMEGA) 
IF (1GR0.GT.M EP9 ■ C MPL X I ER** E 0 ,-S I G9 t OME GA ) 

I F ( IGPO.GT. 1 | ERR - f P9/EP3 

IF IKFL AG( 21 ) .GT .0 ) WRITE (6,121) INT 

ETA - CSQR MU0/EP3) 

GAM - OMEGA*rSORTl-UO*EP3> 

If IKF L AGI 12 I.NE .1) GO TO 9 
NPG - NP 
NMG ■ NM 

00 8 I » I » NPG 



9 7 

98 

99 
1 00 
101 
102 
103 
109 
105 
I 06 
107 
1 08 

109 

110 
I I I 
112 
113 
119 

112 

I 1 7 

lie 

119 

120 
121 
122 
1 23 
I 29 
125 
I 26 

127 

128 
1 29 
1 30 

m 

) 33 
I 39 
135 
1 36 
I 37 

138 

139 

190 

191 

192 

193 
199 



C 

C 



c 



c 



c 



c 



c 



c 

c 



c 





X G I I ) ' 


XII) 




Y G 1 1) = 


Y I I ) 


0 


Z G I I ) ■ 


Zl I ) 


9 


00 )0 1 

X I ) ) ■ 


r 1 t NPG 
X G 1 I ) 




Y 1 1 ) - 


YGI 1 ) 


10 


Z 1 1 1 - 


ZGI I ) 



NP « NPG 
NM * NMG 

If ( IGPO.LE.O) GO TO 15 

SFT UP (MAGE FOR GRUUNO PLANE 

Z M ) N ■ Z( 1 ) 

IWL ' 0 
K e 0 



IF (Z)il!lKZM)N) zmin-ziii 

I F 1 ! Zl ) ) 1 Git V ?E-60) GO TO 11 
IWL ' I WL ♦ I 

11 continue 

IF (ZMI N.LE.HGI ) GO TO 12 
WRllE (6,117) ^ , 

IF ( IFL Ag.EQ.1 ) GO TO I 
IF ( IFL AG. EQ. 2) STOP 
MSG « 2 
GO TO 9 



12 00 13 J=1,NM 
K « J ♦ NM 

) F * * | A I J ) t GT 1 IWL ) 1 Al K ) ■ I A I J ) ♦NP- I WL 

13 I 8 1 K ) ' I 8 ( J > ♦NP- X WL 

1WLP - 1 WL ♦ 1 



00 19 1 -1WLP.NP 
J - | ♦ NP- | WL 

xui ; *}}| 



YU. 

19 Z I J ) 



-Zl 1) 



KNH - NM ♦ I 
NM - 2 • NM 



195 

196 

197 

198 

199 

150 

151 
) 5 2 
153 
159 

155 

156 

157 

158 

159 

160 
161 
) 6 2 
163 
169 

165 

166 

167 

168 
169 
I 70 
171 
1 72 
173 
179 

175 

176 
I 77 
I 70 
I 79 
1 80 
181 
182 
183 
1 09 
185 
I 86 
107 
I 88 
109 

190 

191 

192 



26 



NP = 2 * NP - | W L 

15 CALL SHUT ( IA, Id, l 1 , 12, 13. JA* J0 , MO ,N0, NM , NP , N,NAX , MIN, 1CJ«1NH1 
IF lHAX.lt. 41 30 TO 16 

WRl Tt I6.MI 
IF HFLAo.LO.il GO TO 1 
IF (lFLAG.EU.zl STOP 
HSG * 2 
GO TO <i 

16 IF 1H1N.GE.11 GO TO 17 
WRITE 16,72 1 

IF 1 l FL AG .FO. 1 1 GO TO l 
IF HFLAG.Eu.2l STOP 
MSG * 2 
GO TO 4 

17 WRITE 1 6,561 

IF (MAX . GT .4.0R.MIN.LT. 1 . OR.N.GT • 1CJ 1 GO TO 54 
112*= 1 

IF 1L0A0.GT.01 30 TO 19 

00 18 1=1. NM 
10 2 L D 1 l I = l0.,0. ) 

19 IF INGEN. GT. 01 30 TO 21 

00 20 1=1 . NM 

20 VGUI ■ 10..0.) 

21 KN * NH 

IF 11GR0.GT.0) KN • NH/2 
J = 1 

ANTENNA CALCULATIONS 

[F UOAO.LE.OI GO TO 24 

IF (KFLAGI24l.GT.OI GO TO 22 

00 23 J = 1 , KN 



22 



00 23 1 *1 • LO A 0 
K ■ L 2 0 4 II 

IF l( IAI Jl.EQ.K) .ANO. IKFLAGI 14I.GT.0) I l LO I J I ■ l LL 0 ( 1) 
IF 1 KFL AG ( 241 . GT . 0 ) 2 L 0 ( K I * l L L 0 I I I 

IF 1 1KF LAGI 141.GT .01 .AND.I IGRO.GT.O)) 2 L 0( J *KN 1 - 2 LO ( J ) 
( (KFLAGI24l.GT.0l . ANO.I 1GR0.GT .oi 1 2 LO 1 K *KN ) - 2L0( K ) 
I JGRO.LE.Ol 30 TO 23 



IF 

IF 

IF 

IF 



l(K. LE. 1WL ) . ANO. (KFLAGI 14) .GT .01 ) 2 LO ( J ♦ KN I - ( 0 . , 0. I 
( 1AIJ1.LE.IWLI .AN0.IKFLAGI24 I.GT.O) ) 2L 01 K*KN ) * I 0. , 0 • ) 



23 CONTINUE 

24 IF (NGEN.LT.OI GO TO 27 



193 

194 
l 95 

1 96 

197 
100 
I 99 
200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 
217 
210 

219 

220 
221 
222 

223 

224 

225 

226 
22 7 
220 

229 

230 

231 

232 

233 

234 

235 

236 

237 
230 

239 

240 



KN = NH 

I F I IGRO.GT.O) KN - NM /2 
J ' 1 

IF IKFL AGI 23 ) . 3 T .01 GO TO 25 

00 26 J *1 , KN 

25 00 26 I = I .NGEN 
K - K Gt N 1 I ) 

IT 1 ( IA I J I .FU.Kl .AND. (KFl AGI 131 .GT.Oll VG ( J ) - VOL T ( I I 
IF IKFLAC.I 23 I. 3 T .01 VG ( K l • VOL T ( l ) 

IF ( (KF LAGI l 31 .GT.OI . ANO. I l GRO.GT .01 I V G I J ♦ KN I - V G I J I 

1 F I IGWO.Lt-.O j 30 TO 26 

• IF II K. I E . iWLi . ANO. IKFLAGI 1 31 .GT ,0 I I VGI J *KN I - 1 0 . , 0 | 

IF << lAUl.t E.IHL) .ANO. (KFLAGI 231 , GT.OI I VG IK*KN ) • ( 0 . « 0 . ) 

26 CONTINUE 

27 CALL SC* AN T ( I A t | p , i NH , 1 NT . 1 SC , I I , 1 2 , 1 3 , J A . Jfl , HO , N , NO *NM , NP , AM , 6 M , C 
1 .CGD.CHm, O.EP 2 .EP 3 .ETA.FH 2 , GAM,SGO,X,Y, l ,JlO, 2 S,ERR, 1 GR 0 ) 

IF (N.GI.O) GO TO 28 
IF (IFLAG.EU. 2 I STOP 
HSG = 2 

IF ( IFl AG.EO . 1 I GO TO 1 

GO TG 4 

20 IF (NGFN.LE.OI GO TO 36 
WRITE ( 6 , 75 ) 

WRITE ( 6,761 
WRITE ( 6,771 
WRITE 16.021 

CALL GA Nil ( I A , 18 . INH, IWR, I 1 , 12 . 13 , 112 , JA, J 0 . MO , N , NO ,NM , AM , C , C J ,CG 

I » CHM, 0 , tF F ,GAH. GG,CGO. SGO , VG* V 1 i « 2 i 1 « 2L 0 * 2S * I GRO ) 

WRITE ( 6,571 EFF.GG. 2 ll 
NEAR F I FI 0 

IF (INFAR.LE. 0 I GO TO 30 
WRITE ( 6 * 75 ) 

WRITE ( 6 , 70 ) 

WRITE ( 6 , 77 ) 

00 29 1 * 1 , INEAR 
XP - XNP( | ) 

YP • YNPI I | 

IP • 2 N PI I I 

CA( L ON FI 0 ( IA , I 0 , INH. I 1 . I 2 , I 3 . HD.N.NO, NM, AM.CGO, SGO,ETA,GAM, CJ, 0 , 
1 X.Y .2 ,XP,YP. ZP.EX.EY.eZ.lGRfi.EflRT 
WHITE ( 6 , 50 ) XR,YP, 2 p 
WRITE ( 6 , 59 ) EX,EY,EZ 
CONT INUE 



29 



241 

242 

243 

244 

245 

246 

247 
240 

249 

250 

251 

252 

253 

254 

255 

256 

257 
250 
2 59 
260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 
2 74 

275 

276 

277 
270 
279 
200 
281 
282 

203 

204 

205 

206 
2d 7 
208 



27 



30 



FAR F 1 F L0 

IF I IGA IN.LE .01 



00 31 |*|, J60 



GO TO 36 



31 



OAIYIII J 
0ATY7I I 
DA TY3I | | 
OAT Y* ( 1 ) 



WRITE 
WRITE 
WRITE 
W«| IE 
INC «= 
NPL 
I F 
IF 



16,75) 
<6, 79) 
(6,7 7) 

( 6,62) 

* 0 
= - 1 
(KFL AG( 16 I , 
(NFFP.EO.I ) 



.EQ.l) WRITE (6,691 
. : GO TO 32 

N PH A * ( PhAF-PHA 1 I /STEP»l 
NTHA • ( THAF-THAI J/STEPM 

GO TO ?* 

32 IF ( AFF I .GT.500. ) GO TO 33 
NPL • I 
NPHA * 360 
NTHA * | 

PHAI « 0. 

THAI * AFFT 
STEP ■ I. 

GO TO 



33 NPL 
NPHA 
NTHA 
PHA | 

THAI 
STEP 
3<i PH 

DO 35 r, = |.r,rn 
PH ■ PH ♦ S I E P 
TH » THAI-STE P 
00 35 1 -I , NTHA 
PHSPH * 0. 
PHSTH - 0. 



2 

■ 1 

360 
AFFP 
0. 

I . 

PHAI -STEP 
K = I . NPHA 



TH ■ TH*STEP 

IF me.RO.GT.O) .ANO.( 1 TH.GT .90) .ANO.( TH.LT .270) ) I GO TO 35 

CalT - ‘ ‘ * * 1 ‘ ~ - - - 



- - . - - - - - . . _ _ __ 

L GFFIU I I A, | 0. INC, INN. I WR. M, 12.1 3, | 12 t H0,N,N0.NH, AH.ArSP.ACST 
I ,C,CG0. CG.CJ.CHH.dfECSP.EdST ,tP ,E^P, E TT.EPfrb, E Pt S .£ T P S , E T T S.f.G. 

2GPPiGTT.PH,SG0,SLSr,SCST,SPPM,SPTM,STPH,STTM,TH,X#y,Z,ZL0,ZS,ETAfb 
3 AM , ERR , iGRD) 

ETMAG • CA8S1ETTS) 



269 

290 

291 

292 

293 
29* 

295 

296 

297 
290 

299 

300 

301 

30 2 
303 
30* 

305 

306 

307 
306 

309 

310 

311 

312 

313 

31 A 

315 

316 

317 
316 
319 
3 20 

321 

322 

323 
32* 

325 

326 

327 
326 

329 

330 

n\ 

333 

33* 

335 

336 



E PM AG = CABS! EPPS ) 

IF (ETHAG.GT . 1 . E - 3 2 ) PHSTH = 5 7.29 5 7 79*ATAN2U 1MAGCE TTS).REAL(ETTSI) 
I F ( E PHA G. G T • I . E - 32 ) PHS PH = 5 7. 29 57 79*AT A N2( A IHAG1 E P PS I , RE AL1 E PP S ) ) 
IF (NPL. LW. I) DATY1 (K)-EPMAG 
IF (NPL. tO. I) DATY2IK j*ETHAG 
IF (NP1 .tO. 2) OA T Y l I | ) =EPHAG 
IF (NPL.EQ.2) 0ATY2I 1 )*ETHAG 
1 f (KFl AG ( 16) .ME • I ) GO TO 35 

WRITE (6,60) TH,PH, GTT, GPP, ETTS,EIHAG, PHSTH, EPPS, EPHAG, PHSPH 

35 CONT 1NUE 

WRITE ( 6,561 

IF (NPL .1 E .0 ) GO TO 36 

CALL POLPRT (|,UATYI) 

CALL POLPRT (2.0AIY2) 

0 A( K SCATTERING 

36 IF (ISCAT.LE.O) GO TO 5* 

WHITE 16,75) 

WP HE (6,00) 

WPITE (6,77) 

WRITE (6,02) 

L * 0 
NPL = - I 
INC * I 

IF (NBAP.EO. I ) GO TO 37 
NPH I * (PHIF-PHI | I/STEP’ I 
NIHI * ( THIF-TH1 1 l/STEPM 

IF r I WR.LE.O) WRITE (6,62) 

GO TO 39 

37 IF ( A8A T. GT.500. I GO TO 38 



N °L 
N PH I 
NTHI 
PHI I 
THI 1 
STEP 



= 1 

■ 360 

* I 

■ 0 . 
A8AT 
I . 



GO TO 39 



30 NPL 
N PH | 
N TMl 
PHI 1 
TH| I 
STEP 
39 PH ■ 



■ I 

- 300 

- ABAP 

■ 0 . 

■ I . 

PHI l-STFP 



00 *2 K - 1 « NP M 1 
PH - PHfSIEP 
TH - THI l-STEP 



337 
330 
3 39 
3*0 
3*1 
3*2 
3*3 
3 ** 
3*5 
3*6 
3*7 
3*8 
3*9 

350 

351 

352 

353 
35* 

355 

356 
? 5 7 
350 

359 

360 

361 

362 

363 
36* 

365 

366 

367 

368 

369 
3 70 

3 71 

372 

373 
37* 

375 

376 

377 

370 
379 

300 

301 

302 

303 
30* 



23 



on 



DC 42 I • I .NTHI 
TH - TH*SltP 

If (I 1CRD.GT .0) .ANO. I I TH.GT.90) .ANO.( TH.LT.270) ) I GO TO hi 
i • L ♦ I 

CALL G F FLO ( U, 18, INC. INM. 1WR.I I .12,1 3, 1 1 2 » MO , N , N 0 , N H, AH , AC SP , AC S T 
I .C.CGO.CG.C J.CHH.O.ECSP.ECSTteP.f T,tPP, ETT.EPPS.EPtS.FTPS.ETTS.GG. 
?&PPtG r I .PHAiD, StS^tSCSTtSPPM, S^TH t STPH, STTM,TH,XtYti,2L0f 2StETA t 6 
3AM,tHH, IGRD) 

IF I 1 wR . 0 T . 0 ) GO TO 40 

IF (NPL.LI.O) WRITE 16*631 PH ,TH , SPPM, SP TH, STPM, STTM, AC SP , ACST,ECS 
I P»fCST . SC SP » SC S T 
40 C FHI I L ) ■ PH 

CTHUIU • TH 
COATI (L) • EPPS 

COAWU I • E P I S 
C D A f 3 | l I ■ E TPS 
C U A I 4 ( Li • fMS 
IF ( NPL .NE . I I SO TO 4 l 
0 A r v I IK I - CAB SI EPPS) 



C ABSl £ PIS) 
CAbSI E TPS ) 
CAbSI E TTS) 

CABS! EPPS) 
CA0SI EPTS ) 
CABSI E TPS) 
CAbSlE TTS) 



DA nr? I * ) 

0ATY3IK ) 

D A T Y<* I K ) 

GO TO hi 

41 0AIY1I I ) 

0AIY2II) 

0ATY3I I ) 

0 A I Y 4 I I ) 

42 CONTINUE 

WRITE 16,02) 

IF I N PL • L E • 0 ) GO TO 43 
CAtL PCLPRT I7.DATYI) 

CALL POL PR f ( 8 , OA T Y 2 ) 

CALL PCLPRT 19.0ATY3) 

CALL POL P R T I1U.0ATY4) 

1 F 1KFLAG1 I 7) .N§ . 1 ) GO TO 45 

43 WRITE I 6,64) 

00 44 I M , L 

44 WRITE 16,65) C PH I II ) , C THE T ( I ) , C 0 A T 1 I I ) , C OAT 21 I ) , C 0AT3 II) , COA T4 1 1 ) 

BISTATIC SCATTERING 

45 IF (I8ISC.LE.0) GO TO 54 
WRITE 16,75) 

WRITE (6,81) 

WRITE (6,77) 

WRITE (6,82) 



305 

306 

307 
3 00 
309 

390 

391 

392 

393 

394 

395 

396 

397 
390 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 
410 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 



ssifi 

L 1 0 
INC * 2 
NPL • -I 

IF (N01 P.EO. I ) GO TO 46 
NPHS ' ( PH SF -PH S 1 I /STEP* I 
N IMS - ( THSF-THSI J/STEPM 
IF ( IWK.LE.O) WRITE (6,67) 
GO TO 40 

46 If (A0I T.GT.500. ) GO TO 47 
NPL = 1 
NPHS • 360 



NTHS 

PHSI 

THSI 

STEP 



GO TO sb 



1 

0. 

A0|T 

1. 



47 N PI 
NPHS 
N THS 
PHSI 
THSI 
STEP 
40 PH - 



■ 1 

■ 360 

• A 0 | P 

■ 0. 

* 1 . 

PHSI-STFP 



on 51 KM .NPHS 
PH ^ PHfSlEP 
TH c TH SI-SUP 

IF ( ( ILRO.f.T .0) .ANO.l 1 TH.GI .90) .ANO.I TH.LT.270) ) ) GO TO 51 

00 51 I =1 ,NTHS 
TH = THKSTEP 
L * LU 

i C ? L rrn F fr°r 1 r 1 » i N ?1 J * *1 I 1 1 I * » 1 3 » 1 1 2 , MO , N , NO^ NH, AM . ACS P . AC S T 

UnA C / G 5 , l C Si C ^;£ H c 7 9 A ECSP l E f ST fcpp 7 ET T .EPPS, EPTS.E TPS, ETTi.GG. 

3 AM?ERh! IgROI ^ 0 ,SCS*,SCST,SPPM,SPTH,STPH, $TtM,TH,X,Y,i,ZLO, 2 S,EfA ,4 

IF ' ( 1WM .GT .0) GO TO 49 

aq rLwfV LT ^2 ) -,|TE ) 6,63) PH,TH,SPPH,SPTH, STPH,STTH 
49 C PHI ( L ) ■ PH 

CTHEML ) - TH 
* COATI (L I - EPPS 
C OA T 2 ( l ) - EPTS 
C0AT3I L ) • E TPS 
C U A T 4 ( t ) - £ TTS 
IF (NPL . NE . 1 ) GO TO 50 
OATYI 1 K ) - CA0S1EPPS) 



433 

434 

435 

436 

437 
430 

439 

440 

441 

442 

443 

444 

445 

446 

447 
440 

449 

450 

451 

452 

453 

454 

455 

456 

457 
450 

459 

460 

461 

462 

463 

464 

465 

466 

467 
460 

469 

470 

471 

472 

473 
4 74 

475 

476 

477 

478 
4 79 
480 



29 



OAT Y2(K ) « CABS! F P F S * 

0 A I Y 3 ( K ) * CABS! E TPS) 
OA1YMK) - CABS1ET IS) 

50 IF ( NPl .NE .2 I GO TO 51 
OAIYI (l ) * CABS! EPPS) 

UA I y <M I ) * CASSt F PIS) 

DA I Y3tl) - CA0SI E Pf S I 
0ATY4I ) ) • CA0SCETTS) 

51 CONTINUE 

WR I TP (6,02 1 

IE (NPl. IE. 0) GO TO 52 

CALL PCLPRT (3.0ATYI) 

C A i l POLPRI ( 4 , D A T Y 2 ) 

CAH PJLPRT (5.0ATY3) 

CALL PfJLPRT I 6 » 0 A T Y4 1 
IF (KFLAG( I8I.NE.11 GO TO. 54 

52 WRITE (6,66) 



00 53 IM ,L 

53 WRITE (6,65) C PH l ( 1 ) , C T HE T ( 1 1 , C 0 A T 1 ( 1 ) , C 0 A T 2 I l ) , C OA T 3 I l I , COA T 4 I 1 | 



54 IF (lFLAG.FQ.il 
IF IIFLAG.EQ.2) 

00 55 J =1 , 30 

55 K FI AG ( I 1 » - I 



I 



GO TO 1 
STOP 



K F L AG ( 0 ) * . 

KFIAG(20) * I 
KFLAG(?6» = 1 
IF UFLAG.EQ.3I WRITE (6,60) 

IF ( IFl AG.l : Q.6I WRITJE 16,1151 
GO TO 4 

56 FORMAT ( 1 HO 1 

57 TORMAT I I OX , * THE R AO 1 A T I ON EFFICIENCY IS • , FI 5. 7//I0X. • THE TIME-AV 
IFMCEjfpwER INPUI IS ••FI5.7//10X.-THE 4NIENN* IHPE04NCE IS 

5 ^ra^;ET 5 :;^■:5^^ 5 IlG/;, J ;r ,5 • 7/^ox •' Ey •'• F,5 * 7 •' * j ’ * f > 5 - 7 '*° 

60 1 FORMAT 1 3 X , F 5 . I , 2 X , F 5 . I ,3X«E10.4,2X,E 10. 4,2 ( 3X , 3 l E10.4,2XI,F6.1,IX 

WII*HSs“///i ,NG ™ e inc10ent, ' t 'H- , ' , '-*ne wave 



<.81 
<.82 
4 R 3 
4 8 4 
<.85 
<.86 
<.87 
<.88 
<.89 
<.90 
<.91 
<.92 
<.9 3 
<.94 
<.95 
<.96 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 
52? 

523 

524 

525 

526 

527 

528 



62 FORMAT (• INC IOEN T • , T 2 7 , * ECHO AREA S I GM A ' , I 66 , ' A 8 SPH P T I ON • . T 9 0 . • E X 
1 T INCT 1 ON' , I I 14, * SCAf T ERlNG' /• PLANE* ,T26, • ( INCIOENT-SCATTEREOl*, l 
24X,3( 5* ,* CROSS SECT ION' ,6X» /* WAVE * , 5 2 X , 3 ( l OX , * F nR • , l I X ) / * PH| 

3 T HI I A ' . 3X . • PHI -PHI • , 3X , * PH | -THE TA» ,4X, * THE TA-PH1 * ,2X, • THETA- THET A 
4' , 3 ( 6 X , * PH i * • 7 X , • THFTa* , 4 X ) ) 

63 FORMAT ( I X ,2 i F5. I . I XI , 10( FI 0.4, 2X ) ) 



64 FORMAT IT5<..'0ACK$LAITFR[NG*/* INCIDENT* ,T37, ’ELECTRIC FIELO POLAR 
IWATIUN SCAtTFUNG MATRIX*/’ P l A NE ' , T 4 9 , * ( I NC I 0 E N T - SC A T T E R t 0 ) * / 3 X 



, T 4 9 , * PHI ■ 
, 3 X , 4 ( 3 X . 

65 FURMAI ( I X ,2 ( F 5 . I , I Xj . 2X , 41 F I I 



< . ' WAVE* , T23. ' PHI-PHI • 
3 I Hf T A ' / * PHI I 1 1 F I A 



» I *» V I l _ 

THETA* ,T 74, ' THFIA-PHI *,Tl02,* THETA- 
RFAL • .RX , * IMAS* ,0X ) ) 

5 , 2 x * t I I . 5,3X1) 



j rur i \ I a i c if ~ * I i I a | | r A | i i ri l • D r c At t 1 1 • ;i JAI I 

66 FORMAT ( T64.*BISTAT IC* /T3/, 'ELECTRIC FIELD POIARUATIUN SCATTERING 
| MATRIX*/' OBSERVATION* ,150,* ( INCTDENI-SCATTEREO)*/' POINT', I4X. 
2 • PHI-PH I ' , T 49, 'PHI -IMF T A' , f 76, ' THE! A-PHl * , T 101 , •THETA-THE TA ' /• P 
3H| THETA* ,4X,4( 3 X , 'REAL* • 8 X , ' 1 M A G * , 8 X ) 1 

67 FURMAT (» Or,, RSVAT ION* . T2K TCHO ARLA S 1 GHA ' / • PO I NT • . I 2 5 , M INC l 

IDF NT- SC AT TIRED)*/' PH I THETA* ,T14, 'PHl-PHl* ,T24,*PH|-THETA»,T37. 
2 * THE T A-PHI • , T4 0. » I HE T A- THE T A’ ) 

68 FORMAT I IHI y 5X, •CONTINUE EXECUTION WITH THE FOLLOWING AODITIONS AN 
IO/IJR CHANGES'//) 

69 FORMAT (54X, • El FCTR 1C FIFLO l N T E N S I I Y ' / 5 X , ' OE GR t E S • , I l X , ' POWE R G A I 
IN' ,?8X, ' T H ( I A ' , 4? X i ' P H l ' / 3 X * ' THET A • , 3 X v 'PHI * . 7X,» THETA', 8X,'PH|»,1 
2 X , 2 ( 8 X , ' R t A l '»0X, ' |naG',8X, 'HAGN'.SX, 'PHASE')) 

70 FORMAT ( IOX. ' HROR IN OA I A CARO NJMBFR '.12,* EXECUTION STOP 

I PL 0* ••• ** * ' ) 

A WIRE SEGMENT MAYNOT 8F SHAREO HY MORE THAN FO 

DIPOLE MOOES CHECK DESCRIPTION OATA CA 

E XECUT ION STOPPEO 



71 FORMAT ( 4 0 x 

I UR • ' /40X , • • 

2RD * ' / 40X , * * 

3 «' ) 

72 FORMAT ( 4 0 X , * • 

IS * ' / 4 ox , • « 

2R0 • ' / 40 X , • • 

3 • ' ) 

73 FORMAT ( 3 0 X , * A BACKSCAT TF RI NG CALL MUSI 
I CAl L *//50X, 'RE QUEST I GNOR L 0 ' // / / / I 

74 FORMAT I ' | ' ,T50 . 3 7( * •' ) /T50 , **’ , T86, •♦•/ 

1 T 50 , ' * OHIO STATE UNIVERSITY 

2 T50,'* ANTENNA ANALYSIS PROGRAM 

3 T 50 , ' « MOO l F ICO FOR USE AT 

4 T 50 , ' * NAVAL POSTGR AOUAT E SCHOUL 

5 T 50, ' ♦ 17 JULY 1974 

6 T50, '• ’ , T86.»«* /T50, 3 7 ( • I ) 

75 FORMAT ( M ' , t50, 2 9 ( • ♦ ’ 1 /T 50 . • ♦' , T 78 ,' •• 

J 6 FURMAT ( 150, •♦* , l IX. * ANTENNA* .W0, ••• I 
7 7 FORMAT U50,*^' , 9X, * C A L C UL A T 1 6 n S • , I 78 , ' ♦ * / T 50 , * ♦ * ,T78, •♦'/T 50,291 

l ♦' I) 

78 FORMAT I T 50 , ' ♦ ' * 9X , ' NE AR F l EL Cl • , T 7B, • • • | 

79 FORMAT i T 50 , * ♦ * , 9 X , ' F A R FIELO' , T 7 0 , • • • ) 



AN ISOLATED WIRE MUST HAVE AT LEAST TWO SEGMENT 

ANO THREE POINTS CHECK DESCRIPTION OATA CA 

EXECUTION STOPPEO 

BE INCLUOEO FOR A BISTATIC 



♦ •/ 
♦•/ 
*•/ 
♦*/ 
♦ '/ 



1 I 



529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 
56| 

562 

563 

564 

565 

566 

567 

568 

569 

570 

W 

573 

574 

575 

576 
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00 r or hat ! T 50, ' • ' , 7X, • 0ACKSCA t t er ING' , T JB . ' *• ) 

61 FURMM I T50, ' *• , 4X , ' 01 STAT ( C SCATTERING* , T7 fl . ' • ' ) 

02 f UR maT 1////1 

03 F HAT ( T30. ' FREQUENCY 1 MHZ ) • , T 8 1 . E l l . 5 ) 

04 F CK HAT I T 30 1 'MIRE RADIUS 1 M E T E R S ) 1 , T fl l , t 1 I . 5 > 

05 FURMAT (130, 'WIRE CONDUCTIVITY ( HECAHH J S /HE TE R I ' . T fl 1 , E 1 1 • 5 ) 

06 EUrMAT ITJO,»W|RE INSULATED ( N U / Y E S ) • , I 0 5 , ' Y E S * > 

0 / fopmat it3o,»whe insulated inu/yesj • ,t as. ' nu' ) 

00 FORMAT t r 30 . • insulation radius 1 METERS ) ' ,T 0 I , E 1 1 . 5 ) 

09 rOR«AT 1 T30, » INSULATION CONDUCTIVITY MHO S/ HE T E R » • , T 0 1 , E 1 1 . 5 I 

90 FORMAT I T30 t 1 INSULATION DIELECTRIC CONSTANT I REl A T I VE ) ■ , T 8 1 . E 1 1 . 5 I 

91 FORMAT ( T30, ' INSULA I ( ON LOSS TANGENT' ,T0 1,E II .51 

92 FdhMAT IT 30. 'EXTER liH MED I UM* , T fl l , ' FR F E S^ACE') 

93 FORMAT 1 T 30, ' EXTER lUH MfcOlUM CUNOUCTIVITY 1 MHOS/ ME T E R) * . T 6 1 . E 1 1 . 5 ) 

94 FORMAT (T30, 'EXTERIOR MEOIUM DIELECTPIC CONSTANT I RE L AT I VE ) * , T B 1 . 

1 tl I .51 

95 FORMA! (T30, *f XTERIOR MEDIUM LOSS T A NGC N T ' , T 0 1 , E 1 1 . 5 > 

96 FORMAT I T 50 » 1 W I RE S I RUC T UR E ' / / T 2 0 . ' SE G * , 4 X , 2 1 ' NOO E * , 1 9 X , * L OC A T 1 ON * 

I i I 8 X ) / T 2 1 * ' NU • ' , 3 X » 2 1 * NU.',9X,'X l ? l3X, l Y*,13X f 'Z' t 7X)/(T21,12,5X, 
221 i2,5X|F11.5,4X,E II. 5, 4X, Ell. 5, 1 X i ) ) 

97 format (t so. 'Antenna feeds' /t<»o, 'node' , i bx, *vol ts' / t s l . *no.' . I2x, 

I * Rf A L ' , TX,' IMAGINARY'/! T41 ,I2,6X,2(4X,E11.5))) 

98 FORMAT i T 50 , * • ' , 63, 'OUTPUT RE QUE S T E 0 ' , T 7 0 , ' * ' > 

99 FURMAT I T30, 'ST3UC TURE CURRENTS') 

100 FORMAT 1 T 30 * ' FAR F1ELUS FJR PHI VARYING FR 0 M ' , 1 X , F 5 . 1 , ■ TO *^5.1, 
1 'AND THETA VARYING FRUM '.FS.I.' TO ',F5.1/ 

2150,'IN STfRS OF '.F5.1,' OEGRE E S . * ) 

101 FORMAT 1 T3J, '0ACKSCAT TEftlNG FOR PHI VARYING FROM , .F5.1,* TO *,F5. 
11,' AND THETA VARYING FROM ',F5.),' TO «,F5.1/ 

2150,'IN STEPS OF • , F 5 - I • • OLGREES.') 

102 FORMAT 1 T30, ' 01 STA T 1C SCATTERING FOR PH I VARYING FRUM •.FS.l.* TO 
l'.FS.l,' AND THETA VARYING FROM ',F5.1,* TO ',F5.1/ 

?T50.'1n STEPS OF ',F5.I,' OEGREES.') 

103 FURMAT 1T30,'NEAR FlELOS FOR FOLLOWING POINTS I X , Y , l ) ' / 5 01 T40 , 3 1 £ 1 
11.5,5X1 )) 

104 FURMAT IT30,'P10T FUR FAR F I ELO T HE T A = • , F 5 . 1 ) 

105 FORMAT I T 30 , ' PL 0 T FOR FAR F1EL0 PHl='.F5.1) 

IOo FURMAT I T 30 , ' PL UT FOR 01STAT1C SC A T T E ft I NG- F OR THE T A- • , F 5 . 1 1 

107 FORMAT 1 T 30 , 1 PL 0 T FOR 61STAT1C SCATTERING FOR PH 1 - • , F5 . 1 ) 

100 FUR MAT I T 30 , * PL 0 T FU9 0 AC K SCA T T E R ING THE TA« ■ t F5 . 1 ) 

109 FORMAT 1 T 3 0, 'PLOT FOR 8ACKSCAT T ER I NG PHl-'.F5.1) 

110 FORMAT I T30» ' ND OUTPUT OR PLOTS REQUESTED') 

111 FORMAT (//I 

112 fURMAT 1 T50,'*' ,T70, ••' /T50 .29I ' ) ) 

113 FORMAT 1150,'**, 0X, 'INPUT DATA ',T78,**'> 

114 FORMAT IT50.29I * ♦' >/T50, »•' , T70.'*' > 

115 FORMAT (10X, 'SINCE THIS OATA 8 LOCK. OOES NOT HAVE A TERMINATION CAR 
10 A CHANGE CA*0 IS ASSUMEO') 



577 

5 70 
579 

500 

501 

502 

503 

504 

505 

506 

507 
500 
509 

590 

591 

592 

593 

594 

595 

596 

597 
590 

599 

600 
601 
602 

603 

604 

605 

606 

607 

608 

609 

610 
61 1 
612 

613 

614 

615 

616 

617 

618 

619 

620 
621 
622 
623 

6 24 



116 FOPMAT I // 10X .401 ' *' ) /I OX . • THE OESCRIPTION AND THE GEOMETRY OF THE 
1 STRUCTURE '/ lOx, ' MUST 0E STATrO IN THE FIRST OATA 8L OCK . ' / I OX , ' •• ♦ 
2* FXtCUT ION STDRREO •**» ) 

117 FORMAT 1//IJX,»N0 PART OF THE WIRE ST R JC JURE CAN LIE 8EL0W THE GRO 
l UNO PL ANE . ' / IOX , ' ♦♦♦♦EXECUT I JN STOPPED****') 

118 FORMAT I T 50 , * ST RUC T'jKf LOADS' / T 40 ,« NODc ',) 6X. ' OHMS ' /T4 1 , 'NO. 12X 
1 i * AL ' i 7X, * IMAGINARY'/! T4l, 12 »6X,2(4X,cl 1 .5) ) ) 

119 FORMAT IT 50,' SIR UC TORE LOADS’ /T 39,' SEGMENT * ,14X, '0HMS'/T41f , N0' ,12 
lXj'PtAL'.TX,' IMAGINARY* /l T41,I2,6X,2(4X,E11 .5))) 

120 FORMAT (750, 'ANTENNA F FED S' /T39, 'SEGMENT' ,14X,'V0LTS'/T41. •NO*' .12 
IX, 'REAL ', 7X, ' IMAGINARY' /(T4| , 12 ,6X,?(4Xj El l .5)) ) 

121 FORMAT I//T30,' THE NJM0ER OF INTERVALS FOR CALCULATING THE ELEMENT 
lS'/TJO, ' IN T HE IMPEOANCE MATRIX WITH SI MP SONS-RULE INTEGRATION IS* 
2,/T30»l3»*. IF CLOSEO FORM INTEGRATION IS REQUIREO SET INT-0'///| 

122 FORMAt I I 30, ' GRUUNO PLANE I NU/ Y E S ) • , T 0 5 , ' NU ' ) 

123 FORMAT (T30,'GR0UN0 PLANE I NO/Y E S ) ' , T05 ■ * YE S* ) 

124 FOPMAT (T30, 'GROUNO OIELECTPIC CONSTANT I RE L A T l VE 1 • , T8 1 , E 1 1 . 5/ 

1 T30, 'GROUND CONDUCTIVITY I MHOS /ME TE R ) ' , T 8 l , E 1 1 . 5 > 

125 FORMAT IT30, 'GROUNO P L ANE ' , T 83 . ' PERFEC T 1 ) 

126 FORMAT I T30, 'ANTENNA HEIGHT I ME T E RS » ■ , T 8 1 . E 1 1 . 5 ) 

127 FURMAT 1//10X.601 '*')/! OX,' THE WIRE RAOlU* MUST 8E S TA TEO • / 1 OX ,4 0 I 
l'*')) 

ENO 



625 

626 
bZJ 
620 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 
6 39 

640 

641 

642 

643 

644 

645 

646 

647 
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DLNK 



PURPOSE: to compress data to the left by removal of 

the blank spaces on the input data cards. 

METHOD: A (I) character is compared to the blank; and 

if it is true, the A (1 + 1) character is shifted to the A (I) 
position. 

CALLED BY: READ 



CALLS TO: NONE 



c 



c 



SUHROUT INE OLNK 1 A) 
D l HE NS l ON A ( 80 ) 

DATA BLANK/' '/ 

K - 0 



00 l 1-1,80 
J « I ~ K 
A ( J J ■ AIM 

IE (A< I I . EQ. BLANK) K-KO 



IE (K.EO.O) RETURN 
A181-K) - BLANK 

return 

END 



1 

2 
3 
A 

5 

6 
/ 



9 

10 

11 

12 

13 

\^ 
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CBES 



PURPOSE: to calculate the quantity E01 where 

B0 1 = J (z) / J (z) . 

o 1 

METHOD: If the absolute value of the argument for the 

Bessel function is less than 12, B01 is calculated via the 
power series expansion for the Bessel function in the DO 
LOOP ending at statement 3. If greater than 12, the 
asymptotic expression is utilized at statement 4. If the 
magnitude of the complex part of the argument for the Bessel 
function is greater than 20, B01 is set to (0.,-1) . If the 
complex part of the argument is negative, the sign of E01 is 
changed prior to returning to the calling program. 

CALLED BY: SGANT 



CALLS TO: 



NONE 



suhrooi iNe tees li.eoi i 

CTHPLex ARG.CC.CS.fX 

CGMPUX 601 ,i,URHJ,ieRHN t Hi2A,JN(2) 
0A1A Rt n. IM59/ 

IF ICABSI2) .GE. 12.01 GO TO 4 
FACTOR * 0.0 
TEKMN * I 0. ,0. I 
M224 * -0.25*2*2 



TEPMJ * 1 l .0,0.01 


9 

10 


00 3 N P = l , 2 


11 


N = NP-I 


12 


jn(np) « re«Mj 


13 


M « 0 


1* 


* = MO 


15 


TEPMJ *= TERMJ*M22WFLOAT(M*IN*MI 1 


16 


JNINPI * JN l N P I OERMJ 


17 


IF INP.NE.I 1 GO TO 2 


IB 


FACTOR = FACTOR* ) .0/F10ATIMI 


19 


TOWN = TEPMN*TE RMJ*FACT0R 


20 


ERROR = CABS(TFRMJ) 


2) 


If (ERROR. or . 1.0E- 101 GO TO l 


22 


TEPMJ = 0.5*2 


23 

zs 


B01 « JN ( ) 1 / JN 1 2 1 


25 


RE TURN 


26 


Y « A l M AG ( 21 


27 


IF IABSIYI.GT.20.) GO TO 5 


2 B 


APG ■ ( .0, 1 . 1*2 

ex • cexp(aagi 


29 


30 


CC *= EX *| ./EX 


to 


CS - 1 . 0,-1 . J * J EX-1 ./EX) 
boi - ici*cc>/ics-cci 


33 


RETURN 


3* 


BUI - I.O.-I.I 
l F 1 Y .1 T • 6 • 1 BOI - (.0,1.1 


35 


36 


RE TURN 


37 


ENO 


3B 
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DSHELL 



PURPOSE: to calculate the mutual impedance term 

contributed by the dielectric insulation on the surface of a 
thin wire. 



METHOD: The contribution to the impedance matrix is 

calculated utilizing the equation below 



z 

ran 



(6 - &) In (b/a) 

2n jw&& 

2 



/ 



m, n 



?• (1) F' (1) dl , 
ra n 



where z is defined in subroutine SGANT, & is the 
mn 2 

dielectric constant of the insulation, b is the outer radius 

of the insulation, a is the inner radius, & is dielectric 

constant of the external medium, and F is the sinusoidal 

expansion function. 



CALLED BY: SGANT 



CALLS TO: NONE 



SUBROUTINE OSHELL UM,8M,OK,CGOS,SGOS.EP2,EP 1 ETA,CAM,PU.PI2) 
COMPLEX C&OS, SGDS.EP2.EP, ETA, GAM, PU.fr l 2 .Go.tST 
OAI A PI / 3 . 1 <» l 59/ 

GO - GAM^OK 

C ST - ( €P2-tP)*eTA*ALOG<BH/AM)/U.*Pl*EP2*$GDSA$GDS) 

Pll - -CSr*IGO*$GOS*CGO$) 

Pie * C ST * 1 G0*CG0 $♦ $G0 S ) 

PET URN 
ENO 
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EQUAL 



PURPOSE: to determine position (location) of the equal 

symbol on input data card. 

METHOD: The character search begins in the column 

passed to the subroutine. On returning to the calling 
program, the argument passed is the column following the 
equal symbol. 

CALLED BY: READ 

i 

CALLS TO: NONE 



SUBROUTINE EQUAL IN) 

COHMON /A/ A(0O) 

DATA EQUAL/*- 1 / 

K * N 

00 l l *K » 80 
N * lfl 

IF (All ). EQ. EQUAL) GO TO 2 

1 CONTINUE 

N • l 

2 RETURN 
ENO 



1 

2 

3 

4 

5 

6 
7 
0 
9 

10 

IS 
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EXP J 



PURPOSE: to calculate the exponential integral with 

complex limits. 



METHOD: 



The exponential integral is defined as: 




V*1 



dv = E ^ (VI) - E ^ (V2) + j2nn , 



where the integration path is the straight line from V 1 to 
V2 on the complex v plane and 



The integration path is a horizontal line in the w plane 
or an inclined straight line from VI to V2 the v plane. The 
integer n is zero unless this path intersects the negative 
real v axis at a point between VI and V2. When there is 
such an intersection. 




,00 



E t (Z) 



Z' 



a) n = 1 if Im (VI) > Im (V2) 



b) n = - 1 if Im (VI) < Im(V2). 



The term j2nn is caluclated below statement 12 



CALLED BY: 



GGMM 



CALLS TO: 



NONE 
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SUBROUTINE FXPJ ( V l , V 2 , W 1 2 ) 

COMPLEX f C » t 1 5 , S » T . UC * VC , V 1 , V 2 , W 1 2,7 
0 IMF NS ION V(2(), W ( 2 1 ) , L) ( l 6 ) , E(lo) 

OAT A V/ 0.2? 20465 7FOO.O. 1 188932 l FOl , 0 . 299 2 73 63E 0 l ,0.5 7751436E01 ,0. 9 
l8374674t01 l 3.15982d/4t02 I 0.933J78l2E-01 ,0.492691 74 EOO, 0 . 1 2 ( 5 595 4E 0 
2 l .0.2 26 994 95E01 , 0. 366 76 22 7E01 .0.54 25336 6E01 .0.75659162E01 ,0. 1 0( 202 
3 2 8E 02,0. 1 31 3 0232 £02,0. 1 66 54 40 8 E 0 2 , 0 . 20 7 J 6 4 7 9£ 0 2 , 0 . 2 562 3 094EO 2 , 0 . 3 I 
4407519E02,0. 38530683E02 .0.48026086E02/ 

OAT A W/ 0.4 5 8964 60 E 00.0. 4 1 70 00 8 3 E 00 , 0 . 1 1 3 3 7 3 38EOO . 0 . 1 0399 1 97E - 01 ,0. 
126101 72 OE -03, O.89854(91E-O6»O.21823487E00, 0.3422101 7E00 , 0 . 26 302 7 58 
2E00.0. 12642582E00.0.40206865E-01 , 0 . 8 56 3 8 7 7 8 E - 02 ,0. 1 2 124 36 1 E -0 2 . 0 . I 
3 1 (6 7440 E- 03, 0. 64 59926 7E -0 5, 0.2 2 2631 69E- 0 6,0.4 22 74 304E-08 ,0.39 21 89 7 
43E-10,0.14565152E-(2,0.(4830270E-15,0.16005949E-19/ 

DATA 0/0. 224 9584 2E 02, 0. 7441 1568 E 02 ,-0.4 ( 43 1 5 76E 03 , -0 . 7 8 7 543 39 E02 ,0 
1. 11254 / 44 E 02,0. (6021 76 IE 03, -0.2 386 2l95E 03,-0. 5009 460 7E03, -0.6840 70 
25 4E 02,0. 1225477 8E 02, -0.10161 97 6E02. -0.47219 59 1E01 , 0 . 79 7 2 960 1 E 01 ,-D 
3.2(O69574EO2.O.22O4649OEO(jO.09720244EOl/ 

OAT A E/0. 211 53(0 7E 02, -0.3 79 59 78 7 E 03, - 0 . 9 74 092 20E 0 2 , 0 . 1 29006 7 2 E03 , 0 
1 . 1 7949 2 26E 02,-0. 1 2 9 1 093 1 E 03 0 . 55 705 574E 03 . 0 . 1 3 52 4 & 01E 02 ,0 . 1 46967 2 

21 E0 3,0. 17949 528E 02 . - 0 . 3 29 8 1 6 1 4E 00 ,0 . 3 1 0 2 88 3 6E 02 , 0 . 8 1 65 76$ 7E0 1 ,0 . 2 2 
323696 l £02,0.391 2 4 89 2E 02 , 0 . 0 1 63 6 7 99E0 1 / 

2 = VI 

00 (2 J 1MM , 2 
X * REAL! 1) 

Y = A 1 M AG ( l ) 

E 15 - ( .0..0) 

A 8 - CAOSm 

IE (A8.EQ.0.I GO TO 11 

IE (X .GE .0. .AN0.A8.GT.10. I GO TO 10 

Y A - A8 S ( Y ) 

IE (X.IE.O.. AND .YA.GT .10. I GO TO 10 

IE (YA-X.GE.l 7.5.0R.YA.GE .6.5.0R.X*YA.GE .5, 5* OA. X . GE . 3 • I GO TO 2 

IE (X.LE.-9.1 GO TO 6 

(F (VA-X.GE.2.5) GO TO 7 

IF (X»YA.GE.(.5I GO TO 3 

N = 6 • ♦ 3 • * A8 

E 15 « 1 ./( N-l. )-Z/N**2 

1 N - N-l 

E 1 5 * 1 ./(N-l • )-2»E15/N 
IF (N.GE.3) GO TO 1 

£15 * Z*E15-CHPLX( ,577216*AC0G(A8 ),ATAN2(Y,X11 
GO TO 11 

2 J1 - 1 
J2 - 6 
GO TO 4 

3 J1 * 7 
J 2 - 21 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 
( 3 
(4 
(5 
16 

I! 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 
3 Z 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



S * 
YS 



( . 0 , . 0 ) 

• Y* Y 



00 5 ( = J 1 , J2 

xi = vi t ) *x 

CF * W( I )/( X l*X( ♦YS) 

5 S = S*.CHPLX(X1*CF ,-YA*CE) 

GO TO 9 

6 T 3 - X • X- Y • Y 
T 4 - 2 • ♦ X • YA 

T 5 = X*T3-YA*T4 
T 6 = X*T4*YA*T3 

UC - CNPLXOm ) »0( 12) 4XA0I 13)*T3»T5-Ell2l*YA-E(131*T4,E(in*E(12l 
1 • Xt£ ( ( 3 I • T 3 ♦ T 6 *0 ( l 2 ) * Y A *0 ( ( 3 > * T 4 1 
Vt = CNPLXI D( 14) »0( 15)*X*0( (6)*T3*T5-EI 151*YA-E( l6l*T4 f E(14UEI151 
(*XfE((6)*T3*T6*0(15)*YA»D(l6)*T4) 

GO TO 8 

7 T 3 = X • X- Y • Y 
14 - 2 . * X* YA 

T 5 * X* T3-YA* T 4 
T 6 * X*T<,*YA*T3 
T J r X* T5-YA«T6 
T 0 = X* T6fYA*T5 
T9 * X* T 7-YA* T 8 
T 10 = X • T 8 f Y A * T / 

UC = CMPL X ( 0 ( 1) ♦0(2)*X*0( 3)*T3*D(4)AT5*0(51*T7*T9-(E(2)*YA»E(3)*T4 
1*E(4)*T6*-E( 5 > * T 8 I ,E ( I ) ♦£( 2)*X*E(3)*T3AE(41*T5*E( 5 ) * T 7a T (0MD(21*YA 
2*0(3)*T4fr0(4)*T6A0(5)*T8M 

VC = CMPIXO(6)»0(7)*X»0(8)*T3»0I9I*T5»0(10)*T7aT9-(E<7I*YA*E(8I*T 
(4fn9)*T6fE(l0)*T8),E(6)^E(7)*X#E(8)*T3»E(9)*T5»E(i0)*T7AT10f(0(7) 
2 • YA * D( 8 ) * T4 *0( 9) * T6 aO( 1 0 ) * T 8 ) ) 

8 EC * uc/vc 

S = EC/CMPLX ( X , Y A ) 

9 EX * E X P( -X) 

T « EX^CMPLX (COS ( YA) ,-S(N( YA) ) 

E 1 5 * S * T 

(F (Y.IT.O.J E 1 5 » CONJGi E l 51 
GO TO 1 ( 

10 E 1 5 - . 409319/( 2 *. 1930441 a. 421031/(ZM.O2666)A. 147126/( 2*2. 56708)f 
( .206335E-l/( L *4. 90035 )*. 107401E- 2/(2 *0. 1 82 1 5 ) ♦ . 1 5 8654E - 4 / ( 2 * 1 2 . 73 4 
22) ♦. 317031E-7/( 2*19.3957) 

E 1 5 - E15*CEXP(-2) 

11 IF (JlH.EQ.17 Wl 2 • E l 5 
(2 2 * V2 



2 * 
TH 



V2/V1 

• ATAN2(A(MAG(2),REAU2) I -AT AN2( A( HAG( V2 ) , REAL ( V2 M * AT AN2 ( AlHAG 



49 

50 

51 

52 

53 

54 

55 

56 

57 
50 

59 

60 
6 ( 
62 

63 

64 

65 

66 
67 
60 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 
00 
8 ( 
82 

83 

84 

85 

86 
07 
88 

89 

90 

91 

92 

93 

94 

95 

96 



1 (VI ) , PE AL ( VI ) ) 

AB - A8 S ( TH) 

(F (A8.1T.1.1 TH 



.0 



p J Ju'f I ? l , ™ * 6 * 2831853 
IF (TH.CT.-1.) TH « -6.203183 
W12 - W (2-E 1 5*CHPLX( . 0, TH) 



RETURN 

ENO 



97 

90 

99 

100 

m 

103 

104 
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GANT1 



PURPOSE: to consider the wire structure as a 

transmitting antenna and calculate the input impedance and 
current distribution. 

METHOD: If a wire antenna is driven by a voltage 

generator v located at one of the current sampling points 
i 

1 and if displacement currents are neglected. Ampere's law 
i 

yields 



V = v F (1 ) 
m i m i 



where F is the sinusoidal expansion function. Thus, the 

excitation voltages V will vanish everywhere except where 

m 

v is not zero, 
i 



The DO LOOP ending with statement 50 uses the delta-gap 
model defined above to determine the excitation voltage 
CJ (I) for all the dipole modes. These are stored 
temporarily in CG (I) . Then subroutine SQROT is called to 
obtain a solution of the simultaneous linear equations. 

SQROT stores the solution (the loop currents) in CJ (I) . 

In the DO LOOP ending at statement 80, the complex power 
input and input impedance (s) are calculated. The 
time-average power input (PIN) is the real part of the 
complex power input. 

Subroutine RITE is called to make the transformation 
from the loop currents to the branch currents. If IWR is a 
positive integer, RITE will write out the list of branch 
currents. 
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Finally, GANT 1 calculates the radiation efficiency by 
calling subrouinte GDISS to obtain the time-average power 
dissipated in the lumped loads and the imperfectly 
conducting wire. 



CALLED BY: 



MAIN 



CALLS TO: 



GDISS 



RITE 



SQROT 



SUBROUTINE GANT l ( 1A , I B , I NH , I U R , I 1 , | 2 . 1 3 . I I 2 , JA . J8 . MO* N f NO * NM , AH , C 
I ,CJ»CG»CMM,D,EFF iGAHf GG , COD * S GO , VG , Y i t , l 1 1 , ZLD , l S , l GRD I 
COMPLEX YY 

COMPLEX Cl 1 1 ? C J ( l ) > CGD1 l ) , SGD( 1) . VG( 1 ) , Z ID ( I) , Y 1 1 , l l 1 , l S .GAM » CG( 1 ) 
DIMENSION 01 II. IA(|). 16(11. JAU), JB ( 1 I 

DIMENSION lili), 1211). I 3 ( 1 } * MO ( 1NM, 4 ) , ND( 1 ) 

. DO 3 1M ,N 

CJ(L) • ( .0, .0) 

K • JAt II 

DO 2 KK = I • 2 
KA • |A(K) 

K 8 ■ 18(K 1 

JJ ■ K 
FI • |. 

IF (K8. EG. 12(D) GO TO l 
IF 1K8. EO. 11(111 F 1 ■ - 1 . 

CJ( 1 ) • CJ( 1 1 fF I • VG( JJ ) 

GO TO 2 

1 IF (KA.FQ.13I D) F | ■ - i , 

JJ • K ♦ NM 

CJ( 1 ) - CJ( 1 ) ♦F I • VG( J J ) 

2 K - J61 I ) 

3 CONTINUE 



DO * I * I * N 
CG ( 1 ) • CJ1 1) 



CALL SQROT lC.CJf 0.112. N) 
1 12 - 2 
V 1 1 - I .0, .0) 

NNN - N 

1 F 1 I GRD.GT .0) NNN • N/2 



00 6 I = l , NNN 
NN - lAlJH(l)) 

YY - CJ ( I ) ♦CDNJG(CG1 1 I ) 

IF (CABS1YYI .Lr.i.E-20) GD TO 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

Si 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



2 11 - l . /YY 
WRITE 1 6, 8 ) NN, 2 I I 

5 YU = VlhYY 

6 CONTINUE 

IF MHK.GT.OI WRITE (6,7) 

CALL RITE ( 1 A. 18. INM, IWR, I 1.12. I 3 , MO, NO , NM , CJ »CG , 1GR0) 
GG - RE AL ( Y I 1 1 
2 II - l/Yll 
PIN • GG 

CALL GDISS ( AM, CG , CMM .O.DISS.GAH.NM.SGO , ZLD.ZS) 

PHAO • PIN-DISS 

EFF « 1 00 • *PRAD/ P I N 

RETURN 



7 FDPMAT ( 5DX, * ANTE NNA BRANCH CURRENTS*) 

8 FORMAT ( 1 OX, * THE INPUT IMPEDANCE AT NODE 
I F I 5 • 7/ / ) 

END 



* • 13. • I S* , FI 5 . 7, * ♦ J*. 
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GDI SS 



PURPOSE: to calculate the time-average power 

dissipated in the imperfectly conducting wire and in the 
lumped loads. 

METHOD: The time-average power dissipated by the wire 

is calculated in the DO LOOP ending at statement 1 utilizing 
the equation below: 




where R is the surface resistance of the wire and a is the 
s 

radius of the wire. 

The power dissipated by the lumped loads is calculated 
by the DO LOOP ending at statement 3. If the wire is 
perfectly conducting, CMM <0, the first calculation is 
by-passed. 

CALLED BY: GANT1 

CALLS TO: NONE 
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no 



SUBROUTINE GOISS ( AM , CG , C MM , 0 , D 1 S S . G AH , N H . S GO ♦ ZL D ♦ l S I 
COMPLEX CG( 1 ) , SGDCU fZLOU) ,CJA t Cj6,GAM, ZS 
0IMFN5 ION Dili 
DATA PI /3. 1 1 5*5/ 

D ( S 5 * .0 

tr icmm.lc.d.i go to 2 

ALPH = RtAKGAMI 
BETA - A ( M AG ( G AM ) 

RH « REAL! 2S)/< 4.*Pl*AH) 

DO 1 K«1.NM 
OK - Dl K I 

OEN - C AO S ( S GO < K ) )**2 
EAU - EXP( Al PH * UK ) 

C AO - C EADM./EADI/2* 

CBO * C OS ( BE T A*DK 1 
5AU - OK 

IF ( ALP H . NE • 0 • ) SAO«l EAD-1 ./EAD )/ 12.*AIPH) 

580 = OK 

IF (BETA.NE.O.) SBD=SINCBETA*DK»/BETA 
FA = RH*(SAO*CAD-SBD*CBD)/DEN 
FB * 2. •RH*(CAD*SBD“SAD*CBDI/DEN 
CJA « CG(K) 

L - K ♦ NM 
CJB = CG(L) 

I 0CS5 - DC SS+FA*lCABS( C J A 1 **2 *C AB S 1C JB I * ♦ 2 1 ♦ FB* 1 RE At 1 C J A )*REALICJB I 
1 ♦ A1HAG(CJA)*A IMAGI CJB) 1 



2 00 3 J>l ,NM 
K - J*NM 

3 DCSS - DISS^REALIZLOI J 1 1* ( C ABS I CGtJ 11**2 ) + RE AU ZLO 1 K 1 1 ♦ 1 CABS I CG 1K> 
1 ) ♦•2) 

RETURN 

END 



1 

2 

3 

A 

5 

6 

7 

B 

9 

10 

1 1 

12 

13 

14 

15 

{? 

IB 

19 
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GFF 



PURPOSE: to calculate the far-zone field of a 

sinusoidal electric monopole. 



METHOD: 



If an electric line source has length d and 



endpoints at (x ,y ,z ) and (x ,y ,z ), then the coordinates 

111 222 

of any point on the source are 

X = X +1 cosx 

1 

y = y +1 cosy 

■I 

z = y + 1 cosz 

1 



where cosx, cosy, cosz are the direction cosines of the 1 
axis, and 1 is the distance along the source measured from 

the endpoint (x ,y ,z ). Let the current distribution on 

111 

the monopole be 



I sinh y(d - 1) + I sinh /I 
1 ( 1 ) = _! 2 



sinh yd 



where I and I are the endpoint currents. 

1 2 

field of this source is 



E = (cosx cos<5 cost ~ cosy cosQ sint 

$ 



E = (-cosx sint + cosy cost) E 
D 1 



The far-zone 



t cosz sinD) 



E 

1 



where 



E = 

1 



-y r 
ge 



4nr (1-g 2 ) sinh yd 



+ 




ygd y f<i> 

[ (e - g sinh yd - cosh yd) I^e 

yf<2 ) 

+ g sinh yd - cosh yd) 1^ e ] 
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f<i) = x sin© cos$ + y sin© sin# + z cos© 
111 

f< 2 > = x sin© cos# + y sin© sin<|> + z cos © 

2 2 2 

g = cosx sin© cos# + cosy sin© sin<|> + cosz cos© 

and (r,©,#) are the spherical coordinates of the observation 
point. 

In this subroutine the range dependence has been 
suppressed. The far field vanishes in the endfire direction 
where GK = 0. If a ground plane is present (IGRD > 0) the E^ 

equation above is decomposed into the x, y, and z components 

and the reflection coefficients are applied before E and E 

© <j> 

field components are returned to the calling program. 

CALLED BY: GFFLD 

CALLS TO: NONE 
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SUBROUT INE GF F ( X A . r A , Z A . X B , Y B , Z B , 0 , CGO , SGO . CTH , STH.CPH , SPH ,GAM, E T 
1 A , E T 1 , f T2.EP1 , E P 2 , I GRO , E H R ) 

COMPLEX ERR.RV.RH.RR.EX.EY, EZ.EE 

COMPLEX ETIjEtJ.EP1,EP2»GAM,ETA 

COMPLEX GO .L GO t bGO, EGO 

COMPLEX EGFA.EgFB fEGGOfESA, ESB 

COMPLEX CST 

FP - 12.56637 

XAB = XB-XA 

YAB • YB-YA 

ZAB - ZB-ZA 

CA » XAB/O 

CB “ YAB/O 

CG “ ZAB/O 

G * ( CA*CPHfCB*S PH)*STH*CG*CTH 
GK - l.-G*G 
EM - ( .0,.0) 

ET2 - ( .Q,.0) 

E PI > ( .0, .01 





EP2 * ( .0, .0) 


20 




IF (GK. L T . .001 ) GO TO 3 


21 




F A s (X A*CPM*YA*SPH )*STH*ZA*CTH 


22 




FB » (XB*CPHfYB*SPH)*STHMB*CTH 


23 




EGFA • CEXP(GAM*FA) 


2A 




EGFB • CEXP(GAM*FB) 


25 




EGGU ' CEXPI GAM* G* 01 


26 




CSf = E T A/ ( GK • S GO * F P 1 


27 




ESA =' CST*FGF a*( EGGO-G*SGD-CGD) 


2B 




ESB - CST*EGFBMl ./EGGO*G*SGD-CGO) 


29 




(F (1GR0.LE.0) GO TO 2 


30 




RV • (- 1 . ,0) 


31 




RH = ( - 1 . .0) 


32 




IF UGRO.fO.l) GO TO 1 


33 




RR - CSCRT <ERR-STH*STH) 


3* 




RV - -( ERR*CTH-RR»/(ERR»CTHFRR) 


35 




RH ■ (C TH-RR)/(CTHFRR) 


36 


1 


EX - C A* E S A 


37 




EY ■ C B * E S A 


3B 




EZ ■ CG • E S A 
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EE - (tX*SPH-EY*CPH)* ( RH-R V ) 


AO 




EX • EX*RV*EE*SPH 


A i 




EY - EY*RV“ E E*CPH 


A2 




EZ ■ -EZ*RV 


A3 




ESA - EX*CA*EY*CB*EZ*CG 


AA 




EX - Ca*ESB 


A 5 




EY - CB *ESB 


A6 




EZ - CG • E SB 


A 1 




EE ■ f EX*SPH-EY*CPH)* (RH-RV) 


AB 




EX - E X • RV fEE • SPH 






EY - EY*RV-EE*CPH 


A9 




EZ - - E Z* R V 


50 




ESB - EX*CA*EY*CB*EZ*CG 


51 


z 


T - (CA*CPHfCB*SPH )*C TH-CG* STH 
P ■ -CA*SPHfCB*C PH 


52 

53 




EU • T • E S A 


5 A 




E T 2 - T *E SB 


55 




EP1 ■ P • E S A 


56 




EP2 - P*E SB 


57 


3 


CUNT 1NUE 


5B 




RETURN 


59 




ENO 


60 






61 
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GFFLD 



PURPOSE: to calculate the far-field for the thin wire 

structure. 



METHOD: The far-field for the structure is calculated 

from the loop currents. The loop currents are either the 
currents produced by the transmiiting antenna calculations 
of subroutine GANT1 or the currents produced by an incident 
plane wave. 



If the incident field is generated by a distance source 
with spherical coordinates (r^,0^,<|>^) • tlie exc i tat i°n 

voltages induced by a incident plane wave are 



where 



V 

m 



/ 



F E 
m i 



dl 



m 



E = E exp (y r r ) dl 
i 0 0 



where E is a vector constant, f is a vector from the 
0 0 

coordinate origin to the distance source, and r is the 
radial vector from the origin to the observation point. 



The field E is generated by test dipole m when 
M 

radiating in the homogeneous medium. Using the vector 

potential, the field at the distance point (r ,0 ,4> ) is 

0 0 0 



E 

m 



-y r 

jwu e 0 
4nr 

o 




exp(y r f ) 
o 



m 



dl 
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where the radial component is to be suppressed. From the 
above equations, 

4nr /r 

V = 2- e 0 e e 

m jwu 0 m 

If an antenna gain calculation is desired, INC is set to 
zero. PH and TH denote the spherical coordinate direction 
of the distance observation point. The phi-polarized (EPPS) 
and the theta-polarized (ETTS) components of the electric 
field intensity are returned to the calling program. 

If INC = 1, a backscattering calculation is desired. In 
this case PH and TH denotes the incident angles for the 
incident plane wave. These are also the spherical 
coordinates of the distance source. The outputs returned to 
the calling program include absorption, extinction, and 
scattering cross section for each polarization; scattered 
electric field; and echo areas. 

If INC = 2, a bistatic calculation is desired. In this 
case PH and TH denote the spherical coordinate of a distance 
observer. Since this calculation uses the induced loop 
currents (EP and ET) , a backscattering call must preceed 
this calculation. The outputs returned to the calling 
program consist of the scattered electric field components 
and echo areas. 

EPP(I) and ETT (I) denote the phi-polarized and 
theta-polarized far-zone fields of dipole mode I with unit 
terminal current. In a backscattering situation, the 
excitation voltages EP(I) and ET (I) are obtained by 
multiplying EPP and ETT by the constant CJI. Then calls are 
made to SQROT which stores the solution (the induced loop 
currents) in EP(I) and ET(I). RITE is called for the branch 
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currents CG (J) , and GDISS is called for the time-average 
power dissipated in the imperfectly conducting wire and the 
lumped loads. This power is denoted PDISS and TDISS for 
phi-polarized and theta-polarized incident waves, 
respectively. 

In scattering problems, the incident plane wave has unit 
electric field intensity at the origin. GGG denotes the 
time-average power density of the incident wave at the 
origin. ACSP and ACST denote the absorption cross sections 
for the phi and theta polarizations. 

PIN and TIN denote the time-'average power input to the 
wire structure, delivered by the equivalent voltage 
generators VP and VT at the terminals. PIN and TIN apply for 
the phi and theta polarizations, respectivity. The 
time-average power input is regarded as the sum of the 
timer-average power dissipated and the time-average power 
radiated or scattered by the wire. ECSP and ECST denote the 
extinction cross sections and SCSP and SCST denote the 
scattering cross sections. 

The distance field is calculated in the DO LOOP ending 
with statement 7 for scattering situations, and in the DO 
LOOP ending with statement 9 for the antenna situation. 

The radar cross sections (echo areas) SPPM, SPTM, STPM, 
and STTM, are defined as 



2 2ar 

o - limit 4nr e S / S 
r si 



where S and S denote the time-average power densities in 
s i 

the scattered and incident fields evaluated at the origin. 



For an antenna, the following definition is employed for 
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the power gains: 



2 2ar 

G = limit 4nr e S(r,0,<|>) / P 

p r i 

where P , GG, denote the time-average power input and 
i 

S(r,e,<|>) is the time-average power density in the radiated 
field. GPP and GTT denote the power gains associated with 
the phi-polarized and the theta-polarized components of the 
field, respectively. 

The use of the variables JFLAG and KFLAG are described 
in subroutine SGANT. 

CALLED BY: MAIN 

CALLS TO: GDI SS 

GFF 
RITE 
SQROT 
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on 



SUBROUT ( NE GF FLO I 1 A , 1 B , I NC . 1 NM, 1 WR. 1 1 . 1 2. 1 3. 1 1 2 , MD , N , NO ,NM. AH . ACS 
IP.ACST ,C ,CGD,CG,C J,CMM,6,EC$P,ECST,£P,£T , t PP » ET T .EPPS»EPts*ETPS,ET 
21 S.GG. GPP, Or I , PH, SGO.SCSP* SCSI, S PPM* SPTH»STPM,ST I H,TH, X ,Y, l ,ZLD, ZS 
3.F TA.GaM.EKR, IGRD I 
CUHPl.EX ERR 

COMPLEX CJI »ET1 . ET2 ,EP1 • EP2 ,EPPS,ETTS,EPTSfETPS»ZS*VP,VT 
COMPLEX CU I ,CJUI,fcPI l J ,ET1 1 1 .EPPll ) ,ETTU l,ZL01 M 



Cl I I ,CJi l I , EPI 1 J.ETl II .EPPll) , 
ETA.dAM.CGOl 1 l.sdoi l)fCGU) 

N (Mil, 10(1 I , lit 1 It 12(1) , 

VIII , zl 1 » 



1311), NO ( 1 ) » MO ( 1 NH ,4| 



COMPLEX _ 

DIMENSION 

DIME NS I ON Dill, XIII. 

DATA PI , TP/3.141 59,6.20310/ 

CJI * -4.*PI/IETA»GAMI 
GGG = REALI I ./ETA) 

THK * .01 74533MH 
CTH = CUSI THR| 

STH - SINI THR I 
PHR * .01 745 33*P H 
CPU - C OS I PH R } 

SPH * SINI PHR ) 

DO 1 1 = 1 ,N 
ETTI | | - ( .0 , .0 I 

l EPPI 1 I • 1 .0, .01 

DO 3 K= 1 t NM 

KA = I A I K I 

K0 = 1BIKI 

NGRD ' IGRD 

IF IK.LT.NM/21 1GR0--1 

CALL GF F I X I K A ) * Y I KA| * Z ( K A ) ,X(K0) t Y IK0I , Z1K0),D(K) , C GD 1 K ) » 5 GD I K ) , C 
1 TH, ST H, C PH, SPH, 6AM, E TA, ET ) ,£T 2, E Pi, EP2, I GRD , E RR) 

IGRD - NGRD 
NDK * NDIKI 

DO 3 1 1*1, NOK 
1 * MOI K, tl I 
FI = 1. 

IF IK0.EQ.121 II) GO TO 2 
IF IK0.EQ. 11(111 F 1 *- 1 • 

EPPI I I ■ EPPI 1) *F 1*EP1 
ETTI I I - ETTI 1 >»F1*ET1 
GO TO 3 

l IF IKA.EQ. 1311)1 F1--1. 

EPPI I ) * EPPI 1 1 ♦F 1 •£ P 2 
ETTI1I - ETTI 1)*F1*ET2 
) CONTINUE 



3 

4 

5 

6 
7 
0 
9 

10 

11 

12 

13 

14 

15 

16 
17 
10 

19 

20 
21 
22 

23 

24 

25 

26 
27 
20 

29 

30 

31 

32 

33 

34 

35 

36 

37 
30 

39 

40 

41 

42 

43 

44 

45 

46 

47 
40 



EPPS ■ I .0,.0I 
ETTS - I.Of .01 
IF I INC .EQ.OI GO TO 0 
IF ( INC .EQ.2I GO TO 6 

DO 4 I = 1 , N 

ETUI • ETT I I I*CJ1 

EPI1I - EPPI I I *C J1 

CALL SCROT IC.EP.O, 1 12, Nl 
112-2 

CALL SQROT I C , E T . 0 , 1 1 2 , N I 
IF IIWR.GT.OI WRITE 16,10) PH, TH 



IF (1WR.GT.0I WRITE (6,11) 

CALL RITE I l A, 1 0. 1NM, 1WR. I 1 . \£ . 1 J , HU.NU , NM. t K.Ll 
CALL GD1SS I AM,c6.CMM,D f PDl$,GAM,NM,$GDt ZL&.ZSl 
IF | [WR .GT .0 ) WR 1 T E 16,12) 

0. 1NM, IWR. I 1.12. 13, MO.NO, NM.ET.CG, IGRD) 

c&,cmm,d,tdis,gAm,nm,sgd,zld,zs) 



[12.13, MD.NDfNM.EP.CG, IGRD) 



CALL RITE (I A, 10 
CALL GDISS 1AM, 
AC SP ■ PD 1 S/ GGG 
AC S T - TD1S/GGG 
PIN - .0 
TIN - .0 



DO 5 1 " 1 » N 
VP - CJ I * E PP I 1 I 
VT - CJl*ETT(l) 

PIN ■ P I N ♦ RE AL IV P *CON JGIEPII))) 
TIN « T 1N»REALIVT *CON JGIE T I 1 I I I 

ECSP - P1N/GGG 
ECST - TIN/GGG 
SCSP - FCSP-ACSP 
SCSI - ECST-ACST 
EPTS - 1.0, .01 
E TPS - 1.0, .01 

DO 7 1-l.N 

EPPS * E P PS ♦ E P I II *EPPI 1 
EPTS - E P T S ♦ F P I 1MFTTJ1 
ETTS - E T T Sa E T ( 1 ♦ETTI 1 
E TP S - E T P S ♦ E T I I )*EPP( 1 




• J.*tp* $aoj 

" ?.*TP* CA0$ 

- 2.MP* CA0S 

- 2.*TP*(CA0S( 



EPPS|**2 

EPTS 

ETPS •♦2 
ETTS ♦♦Z 



0 DO 9 I- l ,N .... 

ETTS - ETTSK J( 1 I * E T T I 1 
9 E PPS - EPPS»CJI I I • fc PP III 

APP * CA0SIEPPSI 
ATT • CA0SIETTSI 
GPP • 4.*PI*APP*APP*GGG/GG 
GTT ■ 4.*P1*ATT»ATT*GGG/GG 
RETURN 



10 FORMAT 1 1 OX . * ©RANCH CURRENTS ASSOCIATED WITH PLANE-WAVE SCATTERING 
1 FOR THE 1 Nt 1 OE NT ANGLES, PH1* , ,F5.1, 1 AND THE T A ■ 1 , F 5. 1 / / ) 

11 FURMAT I44X, ^CURRENTS In6uCED BY THE PH 1 P ° LAR 1 ZE D WAV E» 1 

12 format 144 X, 'Currents induced oy the theta polarized wave ) 
end 
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GGMM 



PURPOSE: to calucate the mutual impedance between two 

filamentary monopoles with sinusoidal current distribution. 

METHOD: As stated in subroutine SGANT, the mutual 

impedance of coupled dipoles may be expressed as sum of four 
monopole-monopole impedances. This subroutine calculates the 
mutual impedance with closed-form expressions in terms of 
exponential integrals. 

For skew monopoles it can be shown that the 
monopole-monopole mutual impedance is given by: 



t -z z 

i+j n mm 

z = (-1) B[e (F -e G + e G) 
ij jl 12 22 



-t -z z 

n mm 

-e (F -e G +e G)1 
j2 11 21 



where m= 2/i, n = 2/j and 



16 n sinh d sinh d 
1 2 



The functions F are defined by: 
ik 



gz cos V 
i 

F =2 sinh d e E (R + gz cos ijf - gt) 

ik * i i 



k 

where g = (-1) , 
being considered. 



d and d are the lengths of the monopoles 

1 2 

The functions G are defined as follows: 
ik 
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G = E (R 

ik 2 


+ qz + q ' t 
2 


- jq”) + E (R + qz 
2 2 


+ q't 


+ jq”) 


- 


* qz + q’t 
1 


- jq") - E (R + qz 
1 1 


q q't 


+ jq”) 


i 

where q = (-1) , 


k 

q' = (-D , 


and q” = qb + q'c 


with 




b = c cos f and 


c = d/sin J. 


The angle p is the 


angle 


formed 



by the apparent intersection of the two monopoles. This 
will be discussed later in detail. 



ik 

In the above equation for G , t denotes the position of 

an observation point somewhere on monopole 2. R and R are 

12 

the distances from the endpoints of monopole 1 to this 
observation point. Finally, the E functions are defined as 
follows: 



jKq" f 
(a + jq ”) = e J 



a + jq” 

2 -y w 

£ dw 



a ♦ jq" 



where a and q” are real quantities with dimensions of 

length, a is a function of t, a = a (t ) , a = a (t ) and 

1 12 2 

y - jw \ZusT The integral above is evaluated by subroutine 
EXP J. 




PARALLEL SKEW. 



To explain the input data for GGMM , refer to the above 
figure. If the monopoles are parallel, then the new 
coordinate system is defined such that the new z axis is 
parallel to the monopoles. The coordinate origin may be 
selected arbitrarily. Si and S2 denote the z coordinates of 
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the endpoints of the test monopole, T1 and T2 are the 
coordinates of the endpoints of the expansion monopole, and 
D is the perpendicular distance (displacement) between the 
monopoles. The mutual impedance of parallel monopoles is 
calculates in the last part of GGMM below statement 5. 

For skew monopoles, let the test monopole s lie in the 
xy plane and the expansion monopole t in the plane z = D. (D 
is the perpendicular distance between the parallel planes.) 
If the monopoles are viewed along a line of sight parallel 
with the z axis, the extended axes of the two monopoles will 
appear to intersect at a point on the xy plane. Let s 
measure the distance along the axis of the test monopole 
with the origin at the apparent intersection. SI and S2 
denote the s coordinates of the endpoints of the test 
monopole. Similarily, let t measure the distance along the 
axis of the expansion monopole with the origin at the 
apparent intersection. T1 and T2 denote the t coordinates 

of the endpoints of the expansion monopole. Let s snd ^ be 

unit vectors parallel with the positive s and t axes, 

respectively. Then CPSI = s*^ = cos p. The monopole lengths 

are d and d . 
s t 

The output data from GGMM are the impedances P11, P12, 
P21, and P22. In defining these impedances, the reference 
direction is from Si to S2 for the current on monopole s, 
and from T1 to T2 for the current on monopole t. In the 
impedance P , the first subscript is 1 or 2 if the test 

ij 

dipole has terminals at SI or S2 on monopole s. The second 

subscript is 1 or 2 if the expansion dipole has terminals at 

T 1 or T2 on monopole t. The monopole lengths d and d are 

s t . 

assumed positive in defining the input data CGDS, SGD1 and 
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SGD2 



For parallel monopoles, CPSI = 1 or -1. Si, S2, T1, and 
T2 are cartesian coordinates for parallel monopoles and 
spherical coordinates for skew monopoles. For skew 
monopoles, the radial coordinates si, S2, Tl, and T2 tend to 
infinity as the angle tends to zero or n. Therefore, if 
the monopoles are within 4.5° of being parallel, they are 
approximated by parallel dipoles. 

CALLED BY: GGS 

SGANT 

CALLS TO: EXP J 



c 

c 

c 



c 



l *P?? G p| GGHH * Sl •S*>u .12, 0.CG0S.SG01 . SCO? ,CP S I , ETA , G AM , P 1 1 ,p 

l ^;;pl*Tp!^SO| l ?DPS| , :co;?O 0 ,S,S ’ TS, ' IS 2 ' STl,Sr 2 <CD ' BO >C* , SS,SK.TLl.l 

COMPLEX c SI^B, EC, ’ EK?EUE°L ’iczi'isf '«?•?}!• Pl2,P2l I P22 
COMPLEX E<?,?l,Fl2,?| tL ' e *l-,EGEI,ESI,ES2,ETl,ET2,EXpi,EXP 



KPB 



OSO * 6* D 
SGOS • SG01 
l s 2.LT.SI I 
SGOT • SGD2 
J I J J2 • l-T • T I I SGOT • 



S30S ■ -SG01 



• I s i * - S leSM cinf s f *®®* T0 5 

IH : EltflBSW 

So ? .. D CE xp ,G * HM2 ' 

OPSI • CPSI 
JOI - T| 

T 02 « T 2 

fPSS • DPSI*OPSI 

DO/OSQRTU. DD-CPSS 1 

BO - C 0*DP S I 
B - BO 

|B • CEXPIGAM*CMPLXI.O,BII 
EC • CEXPI GAM*CHPLX | .0 »CI I 

00 I K-1,2 

00 1 LM.2 
l EIK,LI - I .0, .01 

T s I ■ roi • roi 
TS? - TD2*TD2 
OPG • 00*00 
SI • SI 

00 4 1-1.2 

E I - I - ll**l 

SOI - SI 

JJS • S0l*S0l 

s T J • 2 . • SD | • TO 1 *0P SI 

ST2 • 2 . *S0I • T02*0PSI 

Rl • OSORTIOPO*SIS*TSI-ST1| 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

2 

3 

14 

5 

6 
17 



19 

20 
21 

22 

II 

IS 

27 

28 

29 

30 

11 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
40 
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R2 - DSQRT ( DPQ»$1$*TS2-ST2» 

EK » EP 

00 3 K * 1 , 2 

FK * l — 1 I * • K 

SK * F K • SO 1 
EL • EC 

00 
FL 

EKL = E K * E L 
XX = FK*PD*FL»CO 
Til «= F L • T 0 1 
TL2 = F L • T 02 
R R ( = R l»SK»TU 
RR2 = K2*SK*Tl2 

CALL EXPJ ( GAM^CM PL X( RRl , -XX ) , GAH*CMPLX ( RR2 , -XX) , EXPA) 
CALL EXPJ (GAM^CMPLXt PR 1 , XX) ,GAM*CMPLX( RR2 ,XX) , EXPB) 
E(K,L) = E ( K, L ) »F 1* (EXPA*EKL*EXPB/EKL ) 

EL » l./EC 



V-il^L 



EK 



1 . /EB 



ZO * SDl*DPSl 
ZC * ZO 

EGZ1 - CEXP1GAHUC) 

RR1 - R 1 ♦ Z 0- T 01 
RR2 > P2UD-T02 

CALL EXPJ ( GAM^RR l ,GAM»RR2 , EXPB » 

RRl • R l - Z 0 * T Pt 
RR2 «= R2-ZUMD2 

CALL EXPJ (GAM*RRi , GAM*RR2,EXPA ) 

F ( 1 » 1 I - 2.*SG0S*EXPA/EGZ 1 
F (1*2) * 2.*SG0S*EXPB*EGZl 
* SI - S2 

CSF - E TA/(16,*PI*SGDS*SGDT) 

p i 1 - CSf*((F(l I l)*E(2,2l»ES2-E<l«2)/ES2)»ET2»(-F(|.2>-E(2,l>»ES2» 
IEU.il/ES2l/EV2) 

P12 - CSl*U -F( | , 1 >-E (2.2) »ES2*E (1 ,2)/ES2> »ET1 HF (1,2 I *E( 2. II »ES2- 
IE ( 1. 1 I /ES2I/ETI I 

P21 - CST»( ( — F (2.1 )-E( 2,2)*ESl+E( l»2)/ESl)*ET2»(F ( 2 » 2 ) ♦ E ( 2.1 ) • E S I — 
1E(1.1)/ES1)/ET2) 

P2 2 « CST't (F<2. 1 )»E(2,2)»ES1-E( l,2)/ESl )*ET1H-F (2.2)-E(2.l J’ESlA 
lECi.il/ESil/EVil 
RETURN 

5 IF (CPS1.LT.0.) GO TO 6 
TA • T 1 
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50 
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53 
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57 

5 B 

59 

60 
61 
62 
63 
6A 

65 

66 
67 
6B 

69 

70 

71 

72 

73 
7 A 

75 

76 

77 
7B 
79 
BO 
B 1 
B2 
B3 
BA 
B5 

06 

B9 
90 
9 1 

92 

93 
9 A 

95 

96 



C 

c 



c 

c 



TB - T 2 
GO TO f 
TA * - T 1 
TB 

SCOT 



7 SI - SI 



T 2 

-SGDT 



00 9 

TJ • 



1 = 1,2 

TA 



00 B J- I * 2 

Z1J = TJ-SI 

R - SORT(OSO»Z1J*Z1J) 

W * R*Z I J 

IF (Z 1J.LT.D.) W • 0SQ/1P-ZIJ) 
V = R-Z I J 
IF (Z l J.GT .0. ) 

IF (J.fcQ. I ) VI 
IF ( J.EQ.l I W1 
E GZ ( 1 , J I • CE XP l GAH*Z l J I 
0 TJ - TB 

CALL EXPJ (GAN^Vl ,GAM*V,GP( 1 )) 
CALL EXPJ <GAM*Wl ,GAM*W,GM( 1 )) 
9 SI - S2 



V • 0SQ/(R»Z1J) 
V 
w 



CST 
Pt 1 



-ETA/I B.'PI'SGOS’SGOT } . 

CSTMGM(2)*EGZ t 2,2) ♦GP(2)/EGZ( 2,2 ) - CGOS*( GMT 1 1 • E G Z C 1,2) ♦GP(1 



l P12 & ^ * C i T i l-GHC 2 )*EGZ 12, 1 )-GP( 2 )/EGZ( 2, 1 UCGOSMGM| 1 ) *E GZ < 1 , 1 ) ♦GP< 

1 P21 E - Z csf MGM( l )*EGZC l,2)*GP< 1 )/EGZ< 1 * 2 ) “CGDSM GM ( 2 ) *EG Z (2,2 ) ♦GPI 2 

1 P22°- * I -GMI l )*EGZ (1 ,1 )-GP( 1 ) / EGZ ( 1 , 1 ) ♦CGOS*( GM(2)PEGZ(2.1 UGP< 
12)/ EG Z< 2,1)) ) 

RETURN 

ENO 
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101 
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1 0 A 
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107 
IDB 

109 

110 

Hi 

113 
1 1A 
115 
l 16 
117 
110 

iiS 
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12A 
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I 30 
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GGS 



PURPOSE: to calculate the mutual impedances between 

two filamentary monopoles with sinusoidal current 
distributions. 

METHOD: The monopole-monopole mutual impedance as 

defined by SGANT is calculated using the equations defined 
in subroutine GNF. The endpoints of the axial test monopole 
s are (XA,YA,ZA) and (XB,YB,ZB), and the endpoints of the 
expansion monopole t are (X1,Y1,Z1) and (X2,Y2,Z2) . DS and 
DT denote the lengths of monopoles s and t, respectively, 
CAS, CBS and CGS are the direction cosines of monopole s, 
and CA, CB and CG are the direction cosines of monopole t. 

The effects of ground for vertical co-linear monopoles 
are applied in a slightly different manner than mentioned 
previously. As with self impedance calculations, the test 
monopole and the expansion monopole are laterally displaced 
by the wire radius. - This lateral displacement is used to 
determine the angle of incident. This technique is applied 
at statement 8. 

If INT = 0, GGS calls GGMM for the closed form impedance 
calculations. Otherwise GGS calculates the mutual impedance 
via Simpson' s-rule integration with the following number of 
sample points: IP = INT + 1. If the monopoles are parallel 

with small displacement, GGS calls GGMM to avoid the 
difficulties of numerical integration. 

Since the point (X,Y,Z) of subroutine GNF lies on the 
expansion monopole t, T is the integration variable and is 
measured from (X1,Y1,Z1). Cl is the current at T for the 
mode with terminals at (X1,Y1,Z1), and C2 is the current at 
T for the mode with terminals at (X2,Y2,Z2) . C denotes the 
Simpson • s-rule weighting coefficient. 
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Below statement 7, GGS performs some analytic geometry 
in preparation for calling GGMM. The remainder of this 
section is concerned with this preparation. 



Let s denote a unit vector in the direction from 



(XA,YA,ZA) toward (XB,YB,ZB). Also let f denote a unit 

vector from (X1,Y1,Z1) toward (X2,Y2,Z2). Then s.f = 

cos 9 = CC where 9 is the angle formed by the axes of the 

two monopoles. Let monopole s lie in one plane P and 

s 

monopole t lie in another parallel plane P . CAD, CBD and 



CGD are the 
sin 9 which 
distance DK 

constructed 



direction cosines of the unit vector d = T x s / 

is perpendicular to both planes. To obtain the 

between the two planes, a vector T? is 

r n 

from ( XA , YA , ZA) to (X1,Y1,Z1) and take DK = 



k .a. 

1 1 



A line is constructed from (X1,Y1,Z1) to the test 
monopole, such that the line is perpendicular to the test 
monopole. SZ denotes the s coordinate of the intersection of 
this line with the test monopole, and the cartesian 
coordinates of this intersection are XZ, YZ, and ZZ. The 

direction cosines of s x d are CAP, CBP, and CGP. 

From the point (X1,Y1,Z1) in plane P^, a line is 
constructed perpendicular to the point (XP 1 , YP 1 , ZP 1 ) in the 
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plane P . This line is parallel with 3 and 
s 

Let "R represent a vector from (XZ,YZ,ZZ) to 

PI denotes U (s x 1) . Si and T1 are defined 
GGMM. 

CALLED BY: SGA NT 

CALLS TO: GGMM 



has length DK. 
(XP1 ,YP1,ZP1) . 
in subroutine 



c 

c 



SUBROUTINE CCS I XA.VA , 2 A , XB , VB , ZB , X l . V l , 1 l , X2 . Y2 . 22 , AMtOS ,CGOS .SGD 
IS.OT.SGOI . INI ,E T A, GAK.Pli , P12.P21 .P22.EKR. iGRD) 

COP PL E X EXl,EVl,EX2,EY2,E2l,Ei2 

COMPLEX PI 1 » PI 2 » P21tP22tEJA 1 EiBiEJl«EJ2*ETAfG4MfClfC2 v CST 

COMPLEX f GO,CGDS» SGDS.SGOT , ER1 »6R2,eTl,£T2 

COMPLEX ERR 

COMPLEX E E ? EXX » EY V 

COMPLEX PP.PX.PY.PZ 

COMPLEX RRl 1 RA2 1 RR3»RR4 f RHl ,RV1 , RH2 »RV2 ,RH3 ,RV3,RH4,RV4 
DATA FP/12. 56637/ 

CA - I X 2— X | I /0 T 
CB * ( Y2-YI I /DT 
CG ■ ( Z2-ZI I ✓ OT 
CAS * ( XB- X A ) / OS 

CBS ■ < YB-Y A ) / D S | ' 

COS - < ZB-ZAJ/OS 

CC - CA*CAS»CB*CBS»CG*CGS 

IF ( (CG.LE..003) .AND* (CGS.LE«.003).AND.( IGRD.GT.O)) GO TO 1 
IF (ABS(CC) .GT..997) GO TO 6 
l SZ * tXl-XAJ»CIS*<Yl-YA)*CBS*l21-ZA)*CGS 



IF ( INT.LE.O) GO TO 7 


21 


INS ■ 2* ( INF/2) 


22 


IF < INS-LT.2I Ins « 2 


23 


IP * l NS ♦ l 


24 


DELT * 0 T/ INS 


25 


T « .0 


26 


OSZ ■ C C* DEL T 


27 


Pll * < . 0 , . 0 ) 


2B 


P 12 - 1 .0,.D> 


29 


P 21 - 1 .0, .01 


30 


P 22 - ( .D.-DJ 
AMS ■ AM*AM 


31 

32 


SGN ■ -l. 


33 




34 




35 


00 5 1 N ■ l » IP 


36 


ZZ1 - sz 


37 


ZZ2 ■ SZ-OS 


3B 


XXZ • X 1 ♦ MCA-XA-SZ^CAS 


39 


YYZ ■ Y !♦ T»CB-YA-SZ*CBS 


40 


III ■ 1 1 ♦ T*CG-ZA-SZ*CGS 


41 


RS ■ XXZ**2»YYZ**2*ZZZ**2 


42 


R 1 ■ SORT ( RS *Z Z l • *2 ) 


43 


EJA ■ CEXP<-GAM*R1) 


44 


E J 1 ■ EJA/Rl 


45 


R 2 ■ SUR T IRS« , ZZ2**2I 


46 


EJB - C EX P < -G AH*R2 ) 


4 1 


EJ2 ■ EJB/R2 


4B 
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E R 1 - 


E J A • SOUS ♦ Z Z 1 •£ J 1 *CGOS - ZZ 2* E J 2 


49 




E R 2 ■ 


-fc J8*SG0S ♦ZZ2*E J2*CG0S-ZZ 1 ♦ E J 1 


50 




FAC • 


. 0 


51 




IF (RS.GT.AMS1 FAC - 1 C A* XX Z ♦ C 8 ♦ Y Y Z »CG* Z ZZ 1 / R S 


52 




E r i = 


CC*U J2-E Jl*CGOSUFAC*ERl 


53 




E T 2 = 


CC*I E J1-EJ2*CG0S l»FAC*ER2 


54 




IF 1 1GR0.L T.OI GO TO 4 


55 




R V 1 * 


( -l . .01 


56 




R Ml - 


1 -l . .01 


57 




RV2 * 


1 -1 . ,0) 


58 




RH2 * 


( - l . ,01 


59 




IF 1 1GRO.E0. 1 1 GO TO 2 


60 




XG1 ■= 


X 1* T • C A- X A 


61 




YG1 • 


1 1 ♦ T *C0- Y A 


62 




ZG1 • 


Zl»T*CG-ZA 


63 




XG2 * 


X 1 ♦ T • C A- X 9 


64 




YG2 « 


Y 1MH6-YB 


65 




ZG2 = 


ZIU*CG-Z8 


66 




RG1 * 


SOR T ( XG1*XG1»YG1*VGI 1 


67 




RG2 - 


SURT UG2*XG2*YG2*YG21 


68 




TT1 * 


A TAN1 RGl/ZGl 1 


69 




T T 2 - 


A \ ANIRQ2/ZG2I 


70 




C TH1 


* COSI 1 T 1 | 


71 




ssthi 


* S1NITT11*$IN(TT1I 


72 




CTH2 


- COSI T T 2 1 


73 




SSTH2 


« SINI T T2 1 * S I N ( T T 2 1 


74 




R R 1 - 


CSORTl ERR-SSTH1 1 


75 




RH1 ■ 


1 CTH1-RR1 1/ 1CTH tfRRl ) 


76 




R V 1 ■ 
R R 2 ■ 


- [F RR*l TH1-RR1 )/(ERR*CTMl4RRl) 
CSURT(ERR-SSTH21 


77 

78 




RH2 - 


(CTH2-RR21/(CTH2*RR2) 


79 


2 


RV2 ■ 
RG - 


-(ERR*CTm2-RR2)/< ERR*CTH2»RR2 1 
SQRTI < X8-XA1MX8-XA) ♦( yA-YAMyS-YA)) 


80 

81 




CPH ■ 


0 


82 




SPH ■ 


0 


83 




IF (RG.LT.l.E-321 GO TO 3 


84 




CPH ■ 


( X8-XA ) /RG 


85 




SPH ■ 


MT8-YAI/RG 


86 


3 


EXX ■ 


F T 1 • C A $ 


87 




EYY ■ 


E T 1 *C8 S 


88 




EE - 


(EXX*SPH-EYY*CPH|»(RH1-RV1) 


89 




EXl - 


E XX *R V 1 ♦ E E • $PH 


90 




EY1 ■ 


EYY*RVl-EE*CPH 


91 




EZ1 - 


-eti*rvi*cgs 


92 




ET1 - 


EX1*CAS»EYI *C8S»EZI»CG$ 


93 




EXX ■ 


E T 2*C AS 


94 




EYY - 


E T 2 • CO $ 


95 




EE • 


IEXX*$PM-EYY*CPM1MRM2-RV2I 


96 



E X 2 - EXX*PV2*EE* SPH 
t Y 2 ■ E VY • R V2-F E • CPH 
E Z2 * - FT2*CGS*RV2 
E T 2 - EX2*CAS*EY2*C8S*eZ2*CGS 

4 C - 3. ♦ SGN 

IF * IN.EU. I .OR. 1N.EQ. IP) OI. 
EGO = CEXP(Gam*10T-T| ) 

Cl = C* < EGO' l ./EGDl/2. 

EGO = ( ExPICAM* r » 

C 2 > CMEGD- l ./EGOI/2. 

PI l • PI UET l*CI 
P l 2 - P l 2 * F T l *C 2 
P 2 1 - P2WET2*Cl 
P22 - P22»ET2*C2 
t * T»OtU 
S Z - SZ^OSZ 

5 SON - - SGN 



CST 
Pll 
P 12 
P 2 l 
P22 

RETURN 

6 SZ1 - IX1-XA)*CASHY1-YA1*C8SHZI-ZA)»CGS 

J>QRM l Xl-X A-SZl »CAS1**2HY|-YA-$ZI *C8SI»*2H Z l-ZA-SZ 1 *C GS I » » 



-ErA»UELT/(3.*FP*SG0S*SG0T| 

C SI • PI l 

C SI *P12 

csr*p?i 

CW*P22 



ORl 
12 I 
SZ2 
OR 2 
121 
000 



SZl ♦0T*CC 

SQRTU X2-XA-SZ2*CAS1»»2H Y 2-Y A- SZ 2*C8 $ I ♦ • 2M Z2-Z A-SZ 2»CGS l • • 



< DRUDR2I/2. 

IF lOrjO.GF .20.*AH.AN0. INT.GT.OI GO TO l 
IF (000. IT. AMI 000 ■ AM 
CAU GGMM (.0,0S,SZI,SZ2»000,CGDS»SG0S»SG0T ,1 .,6TA ,GAH, P 11 , PIZ,P21 
1 » P22 1 

IF (lGRO.LE.O) RETURN 
IF ( IGPO.GT. I ) GO TO 8 
Ml - -Pll 
P l 2 ■ -P12 
P 2 1 ■ - P 2 l 
P 22 - -P22 
RETURN 

7 SS - SORT ( l . - CC • CC ) 

CAO - 1CGS*C8-C8S*CGI /SS 
C 80 - 1CAS*CG-CSS*CA)/SS 
CGO • ( C8S*CA-CAS*C9 l/SS 



97 

98 

99 
l 00 
101 
102 
10) 

104 

105 

106 
107 
l 08 
109 
1 10 
1 1 1 
112 
1 1 ) 
1 14 
1 15 
1 16 
117 
1 18 
1 19 
l 20 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 
1)1 
1 32 
133 
1 34 
135 
1 36 

137 

138 
1 )9 
1 40 

141 

142 

143 

144 
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DK * (Xl-XA)*CAO»m-YA)*C0DMZl-ZA»*C3O 145 

DK « ABS(UK) 146 

IF (UK.LT.AMJ DK ■ AM 147 

XZ - XA ♦ SZ *CA S Ue 

YZ ■ YA ♦ SZ*CRS 149 

ZZ - Z A ♦ SZ • CGS 150 

XPl * X1-0K*CA0 151 

YP1 * Y1-DK*C0O 152 

ZP1 ■ Z 1 -OK *CG0 153 

CAP • CBS*CGO-CGS*C0D 154 

CtiP « CGS*CAD-CAS*CGO 155 

CGP ■ CAS*CBO-C0S*CAO 156 

PI • C AP* ( XP1-XZ ) ♦C0P*(YP1-YZJ ♦CGP*! ZP 1 ~ ZZ ) 157 

T 1 - Pl/SS 150 

SI - T1*CC-SZ 159 

CALL GGMM ( S 1 . S 1 »0S , T l , T l »0 T , OK , CGOS . SGO S t S GOT , C C . E T A ,G AM , P 1 l , PI 2 , 160 

IP21.P22) 161 

RETURN 162 

163 

} AMS * A M* AM 164 

RG = (Xl-XA) MXl-XAJMYl-VA l»lYl-YA) 155 

IF (RG.LT.AMS) PG • AMS 166 

UG * SUAT1(Z1-ZAIMZI-ZAI*AG1 167 

CPH • ARS(Zl-ZA)/OG 16B 

SSPH ■ RG/DG* DG 169 

PR1 ■ CSCRT(ERR-SSPH) 170 

R V 1 - - (EKR*CPH-RR1 )/ (ERR*CPH*RRl ) 171 

Pil * PU*RV1 172 

RG = (Xl-XBI*(Xl-X0Jf IY1-YR IMYI-Y0) 173 

IF (RG.LT.AMS) RG = AMS 174 

UG * SQRTU Z1-Z0 )»1Z1-ZBI»RGI 175 

CPU = A0S( Zl-ZBJ/OG 176 

SSPH * RG/DG*OG 177 

PR 1 = CSORMERR-SSPH) 170 

RVl = - (ERMCPH-RR1 )/(ERR*CPHFRRl ) 179 

P 1 2 - P 1 2 * RV 1 100 

RG = (X2-XA)MX2-XA)MY2-YA)MY2-YA) 101 

IF (kG.LT.AMSI RG * AMS 102 

OG - SOKTUZ2-ZAJMZ2-ZAJFRG) 103 

CPH - ABS1Z2-ZAJ/0G 104 

SSPH • RG/DG*DG 105 

RR1 - CSQRT(ERR“SSPHl 106 

RVl • - ( ERR*CPH“RR1 )/ ( ERR*CPH*RR 1 ) 107 

P21 • P214RVI 100 

RG - 1X2“X0|*(X2-XBJ*1Y2-YB)*1Y2-YB) 109 

IF (RG.LT.AMS) RG * AMS 190 

OG • SQRT1(Z2-Z0I*(Z2-Z0JFRG) 191 

CPH - A 0 Si Z2-ZBJ/0G 192 



SSPH • RG/OG*OG 

RR 1 • CSQRT(ERR-SSPH) 

RVl • -1ERR*CPH-RRI )/ 1ERR*CPHFRRI) 

P22 • P22*RVl 

RETURN 

ENO 



193 

194 

195 

196 

197 
190 
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GNF 



PURPOSE: to calculate the near-zone electric field of 

a sinusoidal electric monopole. 



z E 




1 



METHOD: An electric line source is located on the z 

axis with endpoints at z and z as shown in the above 

12 

figure. Let the electric monopole have the following current 
distribution: 



I(D 



I sinh /(d - 1) + I sinh yi 

, 1 2 

sinh yd 



where 1^ and 1^ are the endpoint currents, y is the complex 



propagation constant of the medium, d = z ? - Z ( is the 

source length. The cylindrical components of the field are 
E (9) =0 and 



-yR -yR 

E(e) = <2= [ (I e 1 - I e 2 ) sinh yd 

4nS sinh yd » 2 

-y R 



+ (I cosh yd - I ) e 1 cos 9 
12 1 

-yR 

+ (I cosh yd - I ) e 2 cos © ] 
2 1 2 
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-Y R 

I cosh yd) e 
2 S 

z 

-Y R 

I cosh yd) e * ] 

t 

where is the intrinsic impedance of the medium and where 
(e,4>*z) denote the cylindrical coordinates in a coordinate 

system centered at the endpoint of z . 

These expressions exclude the field contributions from 
the point changes at the endpoints of the line source, since 
these charges disappear when two monopoles are connected to 
form a dipole. 

Let the coordinate s measure distance along the test 
monopole with the origin at (XA,YA,ZA) . From any point 
X,Y,Z, a line is constructed perpendicular to the monopole. 
SZ denotes the s coordinate of the intersection of this line 
with the monopole. The length of the line is the radial 

2 

coordinate 6, and BS denote 6 . R1 and R2 are the distances 
from (XA,YA,ZA) and (XB,YB,ZB) to the point (X,Y,Z) . 

In the statements above statement 1, the above equations 
are solved; ai*d after statement 1, the cartesian components 

(E ,E , E ) of the field are determined. If a ground plane is 
x y z 

present (IGRD>0) the reflection coefficients are applied to 
the cartesian components before returning tp the calling 
program. 

CALLED BY; GNFLD 

CALLS TO; NONE 



E (z) = [ (I 

4u sinh /d 1 



+ <1 - 
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SUBROU T ( NE GNF < X A , Y A , Z A. X B , Y B . Z 6 , X , Y , l , AH , OS , C GOS , SGO S , E f A , GAM , E X 1 

1 l ,t'Yl ,6 Z1 ,tX2,EY2. EZ2, lGAO.ERRl 2 

CCMPLEX ERR.RV1 ,RH1 ,RV2 f RH2 ,RRl ,RR2,EE 3 

COMPLEX E J A , E JB p EJ1 .EJ2.ER1 , E R2 , 6 S l , E S2 , S GOS . GAM , C S T , CGOS . E T A A 

COMPLEX EX1.EY1 • E Z 1 ,EX2,EY2,EZ2 5 

OATA PI /3. 1 A 1 59/ 6 

C AS = ( X0-XA )/OS 7 

CBS * ( YB-YA I /OS B 



COS = ( ZB-ZAI/OS 


9 


SZ = < X-XAI*CASMY-YA)*CBS* (Z-ZA)*CGS 


10 


ZZ1 = SZ 


11 


112 * SZ-OS 


12 


XXZ - X-XA-SZ*CA$ 


13 


YYZ = Y-YA-SZ*CBS 


14V 


III = Z-ZA-SZ*CG$ 


15 


RS * XXZ**2*YYZ**2*ZZZ **2 


16 


R 1 * SlJRT 4 R S ♦ Z Z 1**2 ) 


1 7 


E J A * CEXPI -GAM*R1 J 


IB 


EJ1 = EJA/Pl 


19 


R 2 * SORT(RS»ZZ2**2) 


20 


EJ8 = CEXP|-GAM*R2) 


21 


E J 2 = EJB/R? 


22 


ESI = E J2-E J ( *CGOS 


23 


ES2 = EJ1~EJ2*CG0S 


2 A 


E« 1 * ( .0, .0) 


25 


ER2 = ( .0, .0) 


26 


AMS = AH* AM 


27 


IF (RS.LT.AMS) GO TO 1 


28 


CIM1 * ZZ1/R1 


29 


CTH2 ZZ2/R2 


30 


ER1 * (EJA*SG0$*EJA*CG0S*CTH1-EJB*CTH2)/R$ 


31 


EK2 =* i -E JB* SGOS *E JB*CG0S*CTH2~E J A*C TH l 1/RS 


32 


1 CST = ETA/14,. •Pi^SGOS) 


33 


E X 1 * CST*(ESl*CAS*ER 1*XXZ ) 


3A 


EY1 - CSr*(FSI*CB$*ERl*YYZ) 


35 


EZ1 - C$T*tE$l»CGS»ERl»ZZ2> 


36 


EX2 - C$T*1ES2*CAS*ER2*XXZ| 


37 


E Y 2 - CST* (ES2*CBS*ER2*YYZ 


38 


EZ2 - CST* (ES2*CG$*ER2*ZZZ 1 


39 


IF UGRO.IE.O) RETURN 


AO 


RV1 - 1 -1. ,0) 


A 1 


R HI - 1 -1 . ,01 


A 2 


R V 2 * 1 - 1 . ,0) 


A3 


RH2 - 1 - 1 . ,0) 


AA 


IF 1 1GR0.EQ. 1 1 GO TO 2 


A5 


HI - S0RT1 1 XA-X> • ( XA-X) ♦! YA-Y1 * 1 YA-YM 


A6 


R 2 - SORT! 1XB-X)*(X8-X| MY8-Y)*(Y8-Y1 1 


A 7 


TH1 - ATAN1RI/1 ZA-Zll 


A8 



A T AN ( R 2/ | Z B - Z ) ) 

C l g R l' ! E ? H : S } N 1 { H J J : S S ! N ! { H J i ! 



TH2 
RR 1 
R R 2 
R V 1 
KM 1 
R V2 
RH 2 

2 RG 
C PM = 6 
S PH * 0 

IF IRC. LU 1 .E-321 GO TO 3 
CPH = ( X B- X A ) / R G 
SPH - ( YB- YA)/RG 

II! : !Ii:?- v i;!f* CPH 
II! : iII2i-ZwS!* CPM 

RETURN 

ENO 



I YA- Y B 1 1 



A9 

50 

51 

52 

53 
54. 

55 

56 

57 

58 

59 

60 

n 

63 

6A 

65 

66 

67 

68 

69 

70 

71 
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GNFLD 



PURPOSE: to calculate the near-zone electric field 

intensity at a given point. 

METHOD: This subroutine calls GNF for the near-zone 

field of each wire segment, and sums over all segments to 
obtain the near-zone field of the wire antenna. FI is used 
in a manner similiar to FI of subroutine SGANT. CJ(I) is the 
loop currents calculated by subroutine GANTl. 

The use of the variables JFLAG and KFLAG are described 
in subroutine SGANT. 

CALLED BY: MAIN 

CALLS TO: GFF 



SUBROUTINE GNF10 «l*. IB. INM. 11.12. I 3.H0. N.NO.NN, RH.CGO. SGO ,E TA.CRN 
complex err 

S?«N!*oN J !i!li? 0 iellw 0 illn. ‘3“'* °< 1 '* *“»• v ‘ n ' m 

. EX - ( .6* .0) 

EY * ( .0* .01 

e z * c * d,.d> 

00 2 K* ( » NH 
K A * IA 1 KI 
KB * 1 B € K I 
NCKO - I GKO 

1 GKO - NGRO 



I I * € • NOK 
HD( Kill) 



NOK - ND( K ) 

00 . 

I - MD( 

IF ikb'.eq.uiiii GO TO l 
| F (KB* EQ. II ( 111 { {■•l • 
EX • EX*F 1*EX1*C J ( I I 
EY • EY*F1*EY1*CJI1I 
EZ - EZ *F l *E Z 1 *C J ( I I 

IF ( KX • EQ* I 34 1 1 I FI--I. 
EX - EX*FI*EX2*CJj I J 
E Y - E Y ♦ F I *EY2*C J 4 I I 
EZ - EZ *F I •EZ2*CJI I I 
► CONTINUE 

RETURN 

ENO 



1 

2 
3 
* 

5 

6 
7 

e 

9 

10 

u 

12 

13 

1 * 

15 
6 
7 

16 

19 

20 
21 
22 

ii 

26 

27 

28 

29 

30 

31 

?! 

36 

35 

36 

37 

38 

39 
60 
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LEFT 



PURPOSE: to determine position (location) of the left 

paren symbol on the input data card. 

METHOD: The character search begins in the column 

passed to the subroutine. On returning to the calling 
program the argument passed is the column following the left 
paren symbol. 

CALLED BY: READ 

CALLS TO: NONE 



c 



c 



SUBROUT IN E LEFT ( N » 

COMMON /A/ A ( 60 I 
OAfA PLEFT/M '/ 

K • N 

00 l l *K t 00 

7f"(A(| I.EQ.PLEFT) GO TO 2 

1 CONIINUE 

N • l 

2 RETURN 
ENO 



1 

2 
3 
A 

5 

6 
7 
0 
9 

10 

n 

!i 
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LINECK 



PURPOSE: 



to insert grid characters on the polar plot. 



METHOD: The peroid character (ISM (2)) is inserted in 

the proper position in the statements above statement 4. In 
the statements after statement 4, the grid numbers labels 
are inserted on the horizontal axis. 



CALLED BY: 



POLPLOT 



CALLS TO: 



NONE 



SUBROUTINE LINECK IX,Y) 

THIS SUBROUTINE INSURES ALL GR10 CHARACTORS LIE ON THE POLAR GRID 
COMMON 1SYM.L1NE 

integer y 

DIMENSION ISYMU4). LlNE>130> 

IF IY.fcQ.0) GO 10 3 



K - 0 

IF IX.LT.10 


.0) GO TO S 




9 

10 

11 

12 

13 

14 


SET UP AREAS OF «PER100 N POLAR GR10 POINT 


CHARACTERS 


1 - INF I X 1 






1 - 1 A 8 S ( 1) 






15 


L « A8SIX1 






16 


IF IIZ-ll.GT.0.5) !•!♦! 




l 7 


IF I2.LT.10 


•0.3R.2.GT. 111.0) GO TO 2 




1 8 


LiNcm - lsvMm 




19 


LINEI60I - 


I SYMI 31 




20 


L1NEI62) - 


ISYMI 3 I 




21 


K ■ K ♦ 1 






22 


IF IK.EQ.21 


GO TO 2 




23 


l ■ 122-1 






24 


GO TO 1 






25 


L 1NE <61 ) - 


ISYMI2) 




26 


IF IY.NE.0) 


GO 10 5 




27 

28 
29 


00 * K= 11.111 




LINEIKI - 1 
CONTINUE 


SYMI2 ) 




30 

31 

32 

33 


FILL IN GR10 NJM8ER LABELS ON HORIZONTAL 


AXIS 


34 

35 

36 


L INEI 111 - 


1 SYMI 7 ) 




L1NCI20I - 


1 SYMI 10) 




37 


LINE! 2 II - 


ISYMI 5) 




3 8 


LINEI22) • 


ISYMIU) 




39 


LINEOO) ■ 


1 SYMI 9) 




40 


L INEI 31 1 - 


1 SYMI 51 




41 


LINE1321 - 


1 SYMI 1 l ) 




42 


L 1NE ( *0 ) - 


ISYMI 81 




43 


L <Nt< *1 1 - 


ISYMI 5 1 




44 


L 1 NE ( 4? ) - 


1 SYMI 1 1 ) 




45 


L1NEI50) - 


ISYMI 7) 




46 


LINE! 51 I « 


ISYMI 5 1 




47 


LINEI52) ■ 


ISYMI 11) 




48 


L 1 NE ( 6l 1 « 


ISYMIl) 




49 


L 1 NE 1 70 1 - 


ISYMI 7) 




50 


LINE! 71) • 


1 SYMI 5 ) 




51 


LINE! 721 - 


ISYMI U) 




52 


L1NFI80) - 


1 SYMI 8) 




53 


LINEI81 1 • 


ISYMI 5 1 




54 


L 1 NE 1 82 1 - 


ISYMI l | ) 




55 


LINE 190 I - 


| SYM | 9 1 




56 


L INEI 91 ) « 


ISYMI 5 1 




57 


L1N6192) - 


ISYMI 11) 




58 


LINE1100) - 


ISYMUOI 




59 


L1NEI101I - 


1SYM 51 




60 


LINE! 1021 - 


ISYMI 1 I) 




61 


L1NEI111I - 


ISYMI7) 




62 


CONTINUE 






63 


RETURN 






64 


END 






65 



65 



NUMB 



PURPOSE: to place degree numbers on the polar plot. 

METHOD: The current line which is being printed is 

passed to the subroutine in the calling argument. If this 
line contains degree numbers, these numbers are placed in 
the correct position by the IF statements. 

CALLED BY: PTPLOT 

CALLS TO: NONE 



E 



SUBPUUT INE NUM8 IY) 

THIS SUBROUTINE PUTS DEGREE NUMBERS ON POLAR GRID 

COMMON ISYM.L I N E 
INTEGER Y 

OIHENSION ISYMII4), L I NE 11301 
IF IY.Nf.J/1 GO TO 1 
l I NE l 3 3 1 = I S Y M ( 7 I 
l I NE I 3 4 ) = ISYH18I 
l I NE l 3b I = ISYMI6I 
LINE! 87 ) * ISYH16) 

L INE I 88 ) ^ | syH| | Z J 
l I NE I 89 I * I SYM ( 6 I 

1 IF IY.NE.2II GO TO 2 
l INE< l? I • I SYM ( 7 | 

L INt II 3 I ■ I $YM( | | ) 

IINIIKI - ( S Y M ( 6 I 
l INE ( 1081 = I SYM< 61 
L i NE < 10 91 = I S YM ( 9 ) 

LINE(IIO) = I SYM ( 6 I 

2 IF IY.NE.0I GO TO 3 
LINE! 71 * I SYM l 7 ) 

LINfclfll * I S Y M ( | 3 ) 

L I NE I 9 ) = I SYM ( 6 | 

L INE l I I 3) * | SYMt 61 
l I NE l 1 I 4 ) * | SYMl 6 I 

l I NE l l I 5) * | SYM ( 6 I 

3 I F IY.NE.-2I I GO TO A 
L INF II 2 I - I SYMl 8 | 

l INE II J I - I SYMl 7 I 
L INE I I 4 I ■ | SYMl 6 ) 

LINE! 1081 - ISYMI9I 
LINE! 109) - ISYMI9) 

L INEI I I 01 • ISYMI 61 
A IF IY .NE.-37) GO TO 5 
L I NE j 3 3 1 - I SYMl 8 1 
l I NE I 3 4 j - 1 SYMl | 6 I 



l INf | 3b I 
l I NE ( d 7 ) 
L INE I 88 ) 
L INE I 89 ) 
5 CONT INUE 
ME TURN 
END 



ISYM|6{ 
i SYMiel 



I SVMj 9 j 
I SYMl 61 



3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
lb 
16 

17 

18 

19 

20 
21 
22 
23 
2* 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3* 

35 

36 

Vs 

39 

40 

41 

42 

43 

44 

45 
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NUMBER 



PURPOSE: to convert alpha-numeric numbers to floating 

or fixed point numbers. 

METHOD: After initially determining the sign of the 

number, the DO LOOP ending at statement 6 scans each 
character beginning at HI. The DO LOOP ending at statement 
3 terminates the outer DO LOOP if the character being 
compared is not an alpha-numeric number. The DO LOOP ending 
at statement 5 converts the alpha-numeric number to an 
actual number. Below statement 7, the multiplier correction 
is applied to the floating point number before returning to 
the calling program. 

CALLED BY: READ 



CALLS TO: NONE 



SUBROUTINE NUMBER lNl,N2tX,lX) 

COMMON /A/ A I 60 I 
DIHFNS (ON HI 10 > 

OATA B/ , 0», , I , 1 »2', , 3* I *4*.*5'. , 6S , 7* I # B , . # 9»/ 
OATA AMNUS.PLUS* POINT/'-* t • ♦ • ?'•*/ 

OATA AK.AM,AU/'K' «'M V «’U V / 

N = N 1 
NS 1 GN • 0 
1 I - -l 
IX ■ 0 
1SET - 0 

IF IAIN) .EO.PLUS ) N*NM 
IF |A1N).NE.AMNU$I GO TO I 
NS 1 GN = 1 
N ■ NM 

1 00 6 I r N » 60 

IF I A | I ).NE. POINT) GO TO 2 
1 SET • 1 

GO TO 6 

2 IF 11SET.E0.1) 11 - 11M 

00 3 K=1 , 10 

IF I A 1 I J.EQ.BUII GO TO A 

3 CONTINUE 

GO TO 7 

♦ 00 5 K s 1 * 1 0 
KK = K-l 

IF I A I I l.EQ.BUI I NUMB-KK 

5 CONT INUE 

IX • NUMBMOMX 
N2 - 1»1 

6 CONTINUE 

7 IF INS1GN.E0.1I IX - -IX 
Y - IX 

IF I 1 1. IT. 01 11-0 

X - Y/| 10**11) 

IF I AIN2) .EO. POINT) N2-N2H 
IF (A1N2I.E0.AK) X • XMOOD 
IF AlN2l.E0.AM 1 “ ‘ * 

IF iAIN2j.E0.AU 
ifuaIn2I.eo.ak 

N 1 « N2 
RETURN 



x • x*o. 666001 

OR. (A <N2). E0.AHI.0R.1AIN2I.E0.AU) I N2-N2*l 



X *0. 001 
00 



ENO 



1 
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A 

5 
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7 
B 
9 

10 
1 1 
12 
13 
1* 

! I 

17 

IB 

19 

20 
21 
22 
23 
2* 

25 

26 
27 
26 

29 

30 

31 

32 

33 

34 

35 

36 

37 
3B 

39 

40 
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42 

43 

44 
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47 

48 
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POLPRT 



PURPOSE: to control the plotting of the polar plot. 

METHOD: This subroutine is the main subroutine in the 

polar plot package and is responsible for calling the 
various subroutines of the package. 

The scale factor, S, must be changed according to the 
printer characteristics. The scale factor in this subroutine 
is set for ten, 10, characters per inch for the abscissa and 
eight, 8, characters per inch for the ordinate axis. 
Therefore S = 10. /8. 

After initializing DATAX, DATAY, and X, the input data, 
Y, is scanned to determine the normalizing factor. If this 
normalizing factor is less than 1.E-32, an error statement 
is printed and the plotting is abortted. 

In the DO LOOP ending with statement 8, each line of 
the polar plot is printed after a call is made to PTPLOT to 
establish the ploar grid information. The variable, DIM, is 
used to as a scaling factor for the polar plot. The value of 
1.0 will cause all of input data to be plotted, however, if 
only the values less than one-half of the normalizing factor 
are of interest, then DIM can be set to .5. This will 
enlarge of the center of the polar plot. 

CALLED BY: MAIN 

CALLS TO: PTPLOT 

SART 
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noooo ooooonn o ooooo noon non on noon on o noon noon nnnn 



1NE POLPRT 4 NAME « Y ) 

ISYM. LINE 

ON X4 360), Y ( 360) « DATAXI36D), DA T AY ( 360 ) t LINEI13D), 1$YM4 

ON TITIAI2). TITL242) 

TLA/*PH1 • » • THE T * / 

.0 



SUBROUT 
COMMON 
0 1 ME NS 1 
ihl 
0 4 ME NS 4 
OATA U 
N - 360 
DIM - I 
NS T - I 
KST « I 

S IS SCALE FACTOR OF PRINTER* 

ABSCISSA CHAR. PER INCH / ORDINATE CHAR. PER INCH 

S - 10.0/8.0 

ZERO OATAX ANO OATAY 



00 I 1 A * 1 , N 
0 * 1 A - 1 

DATA XI I A) * 0.0 
OATA YIIAI = 0.0 
I XIIA1 = 0*3.1415927/180.0 



FACTOR IS THE NORMALIZING DIVISOR 
FACTOR «= Y I I ) 

00 2 1A*2,N 

2 IF 4FACT0R.LT.YI I AM FACTOR *Y I 1 A ) 



IF INAME.EO.il TITLMTITLAI I) 

IF INAME.EQ.2) T I TLI *= T 1 TLA < 2 > 

IF 4 (NAME.EQ.3) .OR. INAME.EQ.4) .DR.INAME.EQ.7).OR.4NAHE.EQ.8 )) TlTL 
12(1) *T1TL All) 

IF I INAME.E0.5) . OR . 4 NAME. EO . 6 ) . OR . I NAME . E 0 . 9 ) .OR . 1 NAME. EQ. 10 )) TIT 
IL2I 1) =T 4TLAI 2) 

IF 4 I NAME.EQ.3 I.OR.I NAME. EQ.5) .OR. I NAME. €0.7) .OR . 4 NAME . EO. 9 )) TlTL 
1 2 I 2 M T 1 TLA I 1 ) 

IF l <NAME.E3.4).0R. I NAME. EO .6) . OR . I NAME . EO . 8) .DR . I NAME . EO. 10 1 ) TIT 
IL2I2MT 1TLAI2) 

IF IFACTDR.GT.I.E-32) GO TO 3 
IF INANE. LE. 2) WRITE 16,9) T I TL 1 
IF I NAME. GE • 3 ) WRITE 46,10) TITL2 
RETURN 
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9 

10 

i) 
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14 

15 

\ 6 i 

18 
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28 
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32 
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37 

38 
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41 

42 
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44 

45 

46 

47 

48 



NORMALIZE OATA TO ONE 

3 00 4 1A»| ,N 

4 Yl 4 A ) - Y 1 I A ) /FACTOR 

IF INANE. LE. 2) WRITE 16,11) TITLl, FACTOR riTl , 

IF ( (NAME. GF. 3). ANO. INANE. LE. 61) WRITE 16,13) TITL2,FACT0R 
IF INANE. GE. 7) WRITE 16,12) T1TL2, FACTOR 
FILL OATAX ANO OATAY ARRAY FROM X AND Y ARRAY 

00 5 IAM.N 

OATA XI 1 A I = Y<1AMC0S<X<1A ) 

5 OATA YIIAI- Yl 1A MSINlXI I A) I 

SORT OATA 8 Y ORDINATE MAGNITUOE 
CALL SART I OA T AX , 0 AT AY , N ) 

OATAX ANO OATAY ARE SORTEO 8Y OESENOING MAGNITUOE ON THE OATAY VAL 
SET UP FOR PLOTTING POLAR GRID WITH OATA 

00 8 I Y Y M ,81 
CALL PTPLOT I 1 Y Y , S ) 

LINE IS RETURNED WITH POLAR GR10 INFORMATION 

SET UP * Y 1 8 1 N SIZE UPPER ANO LOWER LIMITS 
ULL IS THE LOWER 8IN LIMIT 
UL IS THE UPPER 8 I N LIMIT 

8 IN - 0 IM/80.0 

ULL - 0 IM-I2MYY-1 M8IN 

UL - UL L* 2*8 IN 

CYCLE THROUGH OATA TO FIND WHICH ONES FALL IN *Y* 81NS 

IF INST.GT.N) GO TO 7 
00 6 JJ-NST.N 



49 
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non 



l 

C 

C 

c 

c 



c 



c 

c 



IF 10ATAylJJI.lT.ULl) GO TO 7 
KST • JJ 

AMAG « SQRTIOATAX(JJI*OATAXl JJ) ♦0ATAYIJJ1*0ATAYI JJ) ) 

CHECK THAT MAGN1TUOE IS NOT OVER OIM 

IF 1 AMAG. GT .0 1M) GO TO 6 

OK IS THE FINAL LINE POSITION FOR THE 

OK • OA TAXI JJI*S**0.0/01N*61 .0 
IF tOK.LT. 10.01 GO TO 6 
K - | N T I OK 1 

K - lABS(K) 

OK - ABS10K) 

IF ( (OK-Kl ,GT .0.5 ) K * K * 1 

IF (UK. LT. 10.0. OR .OK.GT .1 I I .0) GO TO 6 

LINE(K) * 1SYMI*) 

CONI INJE 

CONTINUE 
N ST - K S T ♦ l 

PRINT OUT ONE LINE OF PLOT 



WRITE 16, 1*1 LINE 
8 CONTINUE 



RETURN 



•9 FORMAT ( l OX , 1 A* • * COMPONENT OF THE ELECTRIC FIELO IS LESS»/IOX. 

I 'THAN 1 • E -6* , tHEREFORE THIS FIELO WAS NOT • /I OX , • PLOT T EO. EXEC 
2UT10N WILL CONTINUE AS NORMAL. •//) 

10 FORMAT I 1 OX • ' I HE MAXIMUM VALUE OF THE BISTATIC PATTERN FOR •/ 

1 10X, |A*. ,1A*, • I INC IOENT-SCATTEREOI IS LESS THAN •/ 

2 10X, • 1.E-30.I POLAR PLOT NOT CALLEO.*///) 

11 FORMAT I'l'.lA*.’ ELECTRIC FIELO ANTENNA PATTERN FOR SPECIF1E0 PL A 
INE.* ,9X .• NORMAL 1 1 1 NG FACTOR- '.E10.5) 

12 FORMAT I* 1BISTAT1C SCATTERING PATTERN F OR • . 1 A* , • - • , I A* , •( | NC I OEN T- 
1SCATTERE0) POLAR U AT ION. • ,9X. 'NORMAL 11 I NG t AC TOR- • , E 10. 5 I 

13 FORMAT 1 MBACKSCATTER1NG PATTERN FOR * , 1 A A, • -• , 1 A A , • I INC 1 OENT- SC AT T 
1ERE0I POLARIZATION. • ,9X, »NORMAL WING F AC TOR -• • E 1 0 . 5 I 

ib FORMAT 1 IX , 1 30 A 1 ) 

ENO 
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1 12 
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117 
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12* 
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128 
I 29 
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I 31 
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1 33 
I 3* 
1 35 
1 36 

137 

138 
I 39 
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PTPLOT 



PURPOSE: to establish the grid information for the 

polar plot. 



METHOD: In the DO LOOP ending at statement 1 the 

alpha-numeric characters are transferred to ISYN in order 
to pass via COMMON to other subroutines. In the statements 
following statement 2, the equations for the plotted 
concentric circles are established. Below statement 7 the 
grid marks on the 090-270 axis are inserted. 

CALLED BY: POLPRT 

CALLS TO: LINECK 

NUMB 
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SUBROUTINE PTPLOT <lYY t S> 

C 

C THIS SUBROUTINE SETS UP POLAR GR10 INFORMATION 

C 

COMMON ISYM.IINE 

OIMENSION L INEI 130) , ISYMIlAl, ISYNI1AI 

OAT A ISYN/ IN* , IH. , IH . 1H* , I H/ , 1 HO , 1H I , 1 H 2 , 1 H 3 , 1 H A t 1 H5, 1 H6 , 1 H B t 1 H9/ 
INTEGER Y,YY,W 

SET UP 1SYM FROM ISYN FOR COMMON 



00 I K=1,1A 
ISYMIKI * 1SYNIKI 
I CONTINUE 



CLEAR LINE ANO SET TO BLANK 



00 2 1=1.130 

2 LINE! II - ISYMI3I 

Y • A|- I YY 
If IY.EQ.OI GO TO 7 

SET UP EQUATIONS FOR CONCENTRIC CIRCLES 
YY * Y • Y 

1 = I Y7*2.5/2I*S 
X = 61.0*SQRTI2500.0-ZI 
CALL LINECK IX, Yl 
IF I Y.GT.32.0R.Y .LT.-32I GO TO 3 
X - 6 I . 0* SUR T I 1600. O-ZI 
CALL L I NECK I X, Y | 

3 IF I Y.GT.2A.0K. Y.LT.-2A) GO TO A 
X ■ 61 ,0»SQRT(900.0-2 I 
CALL LINECK (X,YJ 

A IF I Y.GT. 16.0R.Y.LT.-16) 

X - 61 .0*SQKTlA00.0-2» 

CALL LINECK (X,YI 

5 IF IY.GT.8.PR.Y.LT.-8) GO TO 6 
X ■ 61 .0*SQRT( 100-2) 

CALL LINECK IX, Y| 

SET UP EQUATIONS FOR MULTIPLES OF 30 DEGREES 

6 X - 61.0* l.732051*Y*S 
CALL LINECK IX, Y| 

X • 61. 0»Y»S/l. 732051 



GO TO 5 



1 

2 

3 

A 

5 

6 

7 

8 
9 

10 
11 
12 
13 
1 A 

n 

ie 

19 

20 
21 
22 
23 

2 A 

25 

26 

27 

28 

29 

30 

31 

M 

3 A 

35 

36 

37 

38 

39 
AO 

i! 

A3 

AA 

A5 

A6 

A7 

AB 
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CALL LINECK 


IX, Yl 




A9 
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PUT IN POLAR PLOT NUMBER LABELS 


51 










52 




CALL NUMB (Yl 




53 




W - lAbSIYI 






5A 

55 

56 




F ILL IN POLAR PLOT 


AT 000, 090, 180, ANO 270 










57 




IF IW.NE.A0I GO TO 


8 


58 




LINFI55) • 


ISYMI2I 




59 




L 1 NE I 5 7 1 * 


I SYMI 2 1 




60 




LINE159I - 


1 SYMI 2 | 




61 




L I NE ( 6 3 1 - 


1 SYMI 2 1 




62 




L INE I 65 1 - 


ISYHI 21 




63 




L INM 67 | * 


I SYMI 21" 




6 A 


8 


IF (W.NE.32I GO TO 


9 


65 




LINEI56I - 


ISYil 21 




66 




L I NE I 58 1 » 


I SYMI 2 1 




67 




LINEI60I - 


I SYMI 2) 




68 




LINU62I » 


I SYMI 21 




69 




LINEI6AI • 


I SYMI 2 1 




70 




LINE I 66 1 » 


I SYMI 2 1 




7 1 


9 


IF IW.NE.2AI GO TO 


10 


72 




L I Nf. ( 5 7 1 • 


1 SYMI 2 1 




73 




LINE I 59 1 - 


1 SYMI 2 1 




7 A 




L I NE I 60 1 - 


1 SYMI 2 1 




75 




L I NE I 62 1 • 


I SYMI 2 1 




76 




L INEI 631 - 


1 SYMI 21 




77 




L I NE I 65 1 • 


ISYMI 2 | 




78 


10 


IF IW.NE.16I GO TO 


11 


79 




L I NE I 58 1 • 


ISYM12I 




80 




L I NE I 60 1 • 


ISYMI 21 




81 




LINEI62I • 


1 SYMI 2 j 




82 




LINEI6AI - 


I SYMI 2 ) 




83 


11 


IF 1W.NE.0BI GO TO 


12 


BA 




LINE 159 1 • 


ISYMI 21 




85 




LINE (631 • 


1 SYMI2I 




86 


12 


CONTINUE 






87 




RETURN 






88 




ENO 






89 
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BEAD 



PURPOSE: to interpret and translate the input data 

cards. 

METHOD: The program utilizes free format for the data 

cards, that is, the program uses character recognition to 
determine which parameters are being read. In the IF 
statements containing A(1), A (2), A (3) , and A (4), the first 
four characters on the data card are compared to the first 
four letters of the key words. This will determine the type 
of parameters that card contains. The other IF statements 
determine which parameters are being read. 

Subroutine BLNK is called to remove the blank spaces on 
the parameter cards. Subroutines EQUAL and LEFT are called 
to determine the position of the equal character and the 
left paren, respectively. Subroutine NUMBER is called to 
convert the alpha-numeric characters to numbers, either 
fixed or floating point. This numerical value is assigned 
to the parameter just determined. 

A detailed explanation of the data cards is found in 
appendix II titled "USERS MANUAL". 

CALLED BY: BLNK 

EQUAL 
LEFT 
NUMBER 



73 



c 

c. 

c 



SUBROUTINE READ 1 IA. IB, IBISC.ICARD, (GAIN* I GRO. 1NE A R. INT f I SCAT . |WR* 
l iflag.kflag.kgen.load .LZD.MSG.NBAP.NBIP .NFFP.NGEN.NM.NP , ABAP, ABAT. 
2AFFP.AFFT fABIPtAB I T , AM . B N , C MM . E R 2 , E R 3 . E R 6 . F MC • HGT . PHAF , PHA I .PHIFtP 
3HU , PHSF V PHSI t tHAF . fHAl 9 TH(F 9 f HI t . THSF V THSI 9 siG2»SIG3« SlG* t TD2«T63 
6 .VOL T . X » XNP i Y * YNP ,Z,ZLLD,ZNP.STEP) 

CCHMON /A/ AlBOl 
C CM PL E X VOL T 4 1 ) , ZLLDl I I 

DIMENSION 4 A ( II* I B 4 Ll» XIII. Y 1 I I 9 Z 4 I ) 9 KGENl I ) 9 KFLAGIII 
DIMENSIC1N X N P ( 1 I 9 YNPlI). ZNPlI), LZDIII 

OATA AA , AB, AC, AO, AE. A F. Afc.AH. A l . AK . AL . A M A . A N . AO . A P . AQ.AR.AS.AT.AU, 
lAW.AX/'A* . * B * » * C ' * , D , »*E , * , F*. , G , »*H*.*I , » , K’.*L*»*M , .*N*» , U , **P*. 
2' Q* . *P* .'SS'T'.'U'.'H'.’X'/ 

DATA BLANK, COMMA, M I NUS*PLEFT,P01NT.R IGHT, SLANT/* • 9 • 9 * 9 * - • , • I • 9 * . • 

I . * I * t */ * / 

RAO * 5 7.295 7 79 
INI - 6 
I B 4 SC = * I 
IGAIN < -I 
INEAR r -I 
ISCAT - -1 
I HR = - I 

IF IIFLAG.EQ.6) GO TO 2 
IF I MSG.NE.OI GO TO 6 

1 READ 15.76. E NO = 7 2 I A 

2 IF 1 1 A I I I .NE.ACI -UR. I A I 21 .NE.8LANKI.OR.I A I 3 ) . NE . BL ANK ) .OR. I AI6I.NE 
I .BLANK | | GO TO 3 

WRITE (6*761 A 
GO TO 1 

3 WRITE 16.751 

GO TO 5 « 

6 READ 15.76. ENO c 7 2 ) A 

5 ICARO = ICARDM 

WRITE (6.771 ICARO. A 

IF I I MSG.NE.OI .AND. ((All) .EQ.AE) . ANO . I A I 2 ) . EQ. AN ) . ANO. I A I 3 ) . E Q. AO ) 
111 GO TO 70 

IF (IMSG.NE.01.AN0.il A I l I .EQ . AS I .ANO. I A I 21.EQ.ATI . AND. I A I 3 I . EQ. AO I 
I • ANO. I A I 6 | .EQ.API I I GO TO 69 

IF II A I I I .EQ.ACI .ANO. I A I 2 1 . EQ . 8 L ANK I . ANO . I A I 3 I . EQ . BL ANK I .ANO. I A I 6 I 
1 • EQ. BL ANK I | GO TO 73 
IF (MSG.GT.OI GO TO 6 
CALL BLNK (A) 

N - 6 

INSULATION 

IF I I A I II .NE.AI1 .OR. IAI2I.NE.AN1 .OR.IAI 3) .N E . AS I . OR . I A I 6 I .NE.AUll 
IGO TO 10 
KFLAGI20I - 1 



1 

2 
3 
6 

5 

6 
7 
B 
9 

10 
I I 
12 
13 
16 

15 

16 
1 7 
IB 

19 

20 
21 
22 
23 
26 

25 

26 
27 
2B 

29 

30 

31 

32 

33 
36 

35 

36 

37 
3 B 
39 
60 
41 
62 
63 
66 

65 

66 

67 

68 



CALL LEFT INI 

6 . l *. 1 1 U 1 ;U E ; ARI .NE.AAI.OR. IA<Nf2).NE.AO).OR.IAIN»3l .NE 

i ♦ « i # f on i o i 

KFIAGI6I - 1 

call EQUAL INI 

CALL NUMBER IN.N2 .XI , IX) 

BM • XI 

IF I A I N 2 I .EQ.R (GHTI GO TO 6 
I F 1 A( N2I .NE. SLANT 1 GO TO 71 
N • N 2 ♦ I 
GO TO 6 

7 , 1 tt 1 •oa.iainmi.ne.au. or. i a<n*2» .ne.aei.or.iainoi.ne 

l.ALII GU TO B 
KFLAG16I • I 

call equal INI 

CALL NUMBER (N,N2,XI 9 1XI 

ER2 • XI 

IF 4 A 4 N 2 I .EQ. RIGHT) GO TO 6 
IF I A| N2 I .NE. SL ANT I GO TO 71 
N - N2* 1 
GO TO 6 

0 . •°R*<AIN^|| *NE»A0I.0R.IA<N^2I.NE.ANI .OR . IA(N^3I.NE 

I * A 0 1 I GO TO 9 
KFLAGI5I - I 
CALL E jUAL INI 
CALL NUMBER 1N.N2.XI.IX1 
SIG2 * XI 

IF 1 AIN21 .EQ. RIGHT I GO TO 6 
IF 1 AIN2 | .NE .SLANT | GO TO 71 
N * N2> I 
GO TO 6 

9 . *^ , ' { ** * •OR.UINMI .NE.*0).0A.< A IN* 2). NE. AS > .0 A . < A I N ♦ 3 ) . NE 

l.Abll GU TO 71 
KFIAGI7! - I 
CALL EQUAL INI 
CALL NUMBER IN.N2 f Xl. 1X1 
T02 c XI 

IF 1 A1 N2) . f Q .R I GMT I GO TO 6 
IF I AIN2I .NE. SLANT I GO TO 71 
N - N2M 
GO TO 6 

WIRE » 
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51 

52 

53 
56 

55 

56 

57 

58 

59 

60 
61 
62 
63 
66 

65 

66 

67 

68 

69 

70 

r! 

73 

76 

75 

76 
7 7 
7 B 
79 
BO 
81 
82 
83 
B6 
85 
B6 
87 
BB 

89 

90 

91 

92 

93 
96 

95 

96 
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noon 



l AI2I.NE.AI ) .OR.UI 3I.NE.AR) , 
< A (NH< .NE. AA ) .OR. < A < N+2 < .NE. 
• I X I 

GO TO 4 
GO TO 71 



1*1X1 

GO TO 4 
GO TO 71 



OR.UIH .NE • AE ) 1 
AO).OR.(A(N*31 .NE 



10 IF HAill.Nf .AW) .OR. 

1GU TO 13 

CALL LEFT ( N | 

11 IF ( ( A I N) .NE.ARJ .OR. 

1 • A I ) ) GO TO 12 

• KFLAGI2I = 1 
CALL EQUAL INI 
CALL NUMBER (N,N2.X1 
AM = XI 

IF ( A ( N 2 ) .EQ.R1GHTI 
IF ( A ( N2 ) • NE • SL AN T ) 

N = N2» 1 
GO TO 11 

12 IF ((AIN) .NE.AC) .0R.(A(N+1< .NE.A0).0R*(A(N+2(.NE.AN( .0R.(A(Nt3(.NE 
1 * AO I ) GO TO 71 

KFLAGI3I = 1 
CALL EQUAL (N| 

CALL NUMBER (N f N2«Xl« 

CMM - XI 

IF < A(N2I .EQ.P 1GHT) 

IF ( A( N2I .NE. SLANT ( 

N = N2* 1 
GO TO l 1 

EXTERNAL ME01UM 

13 )F (<A( 1I.NE.AE(.0R.(A(2) .NE . AX ( .OR. ( A( 3 ( . N E . AT I . OR . ( A ( A ( .NE . AE ) ( 

1 GO TO 1 7 

KFLAG<0) «= 1 
CALL LEFT (N) 

1* IF (<A<NI.NE.AC( .OR. IAIN* 1( .NE.AOl.OR.I A(N*2< .NE.AN( .OR . ( A ( N *31 . NE 
1 .AO) I GO TO 15 
KFLAGI91 * 1 
CALL ECUAL ( N I 
CALL NUMBER (N.N2.X1. 

S1G3 * XI 

IF IAIN2I.EQ. RIGHT) 

(F (A(NZ).NE. SLANT) 

N - N2* 1 
■ GO TO 14 

15 IF ( ( A (N) .NE.AO) .OR.IA(N* 1) . NE . A I ( .OR . < A ( N* 2 < . NE . AE I .OR • < A < N* 3 ) . NE 
1 .AL1 I GO TO 16 
KFLAGIIOI * 1 
CALL EQUAL I N ) 



} 1 X ( 



GO TO 4 
GO TO 71 
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98 

99 
1 00 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 
1 16 
117 
1 Id 

119 

120 
121 
122 

123 

124 

125 

126 

1 27 

128 

129 

130 

131 

132 

133 

134 

135 

136 
13 ? 
138 
1 39 

140 

141 

142 

143 

144 



CALL NUMBER (N,N2.X1,IX( 

ER3 ' XI 

IF ( A(N2) -EQ. RIGHT I GO TO 4 
IF (A(N2) .NE. SLANT) GO TO 71 
N ■ N2 ♦ 1 
GO TO 14 

lfc |f AS! » A GO , fo 6 71 L * * 0R * -NE.ACM.OR.I A€N*2* . NE . AS ) .OR . f A « N . 3 1 . N E 

kfiagmij • i 

CALL EQUAL (N) 

CALL^NUfSER ( N * N2 « XI « IX) 

IF ( A(N2) .EQ.R 1GHT ( GO TO 4 
IF ( A<N2| .NE. SLANT! GO TO 71 
N * N2» ) 

GO TO 14 

LOAO 

• l Vo TO*! J ,,ne ** lkor *‘*‘ 2 '-ne.*0( .OR.IAOI.NE.AAI.OR.UUI.NE.AOM 
KFlAGim • l 
GO TO 19 

I9 1 l GO J To' ' .OR.(A(2).NE.AMA).0R.<A<3).NE.AP).0R.U<4).NE.AEH 

. KFLAG(24( - 1 
)9 I - I 

CALL LEFT (N) 

20 CALL NUM8EP (N,N2.XI.IX( 
cA tA 

N ■ N2»l 

CALL NUMBER IN, N2, XI, 1X1 
R M AG * X) 

N - N2* 1 

CALL NUMBER IN V N2.X1 9 1XI 
ROEG • XI 

RREAL - RM AG*COS ( ROE G/R AO ) 

• Rmag*siniroeg/raoi 

/LLO< I ( - CMPLX|RREAL,RIHAG( 

LOAO - I 

IF (A(N2(.EQ. RIGHT) GO TO 4 
IF <A(N2).NE. SLANT ( GO TO 71 
1 - I ♦ l 
N - N2M 
GO TO 20 



145 

146 

147 

148 

149 

150 

151 
1*52 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 
) 65 
166 

167 

168 

169 

170 
| 71 



174 
1 75 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 
187 
1 88 

189 

190 

191 

192 



75 



non 



21 KFLAGI24) • -I 
LQAO - -( 

GO TG * 

F OfcOUENCY 

22 JF UAl(J.NE.AF).OR,(AI2).NE.AR|.OR.(A( 3 ) . NE . AE ) . OR. I A ( 4 ) . NE . AQ ) ) 

I GO TO 2 3 

KFLAGI1) « 1 
CALL L F T T IN) 

CALL NUMBER (N,N2,Xl t (X) 

FMC • XI 
GO TG 4 

plot 

2 3 IF KAI 1) .NE.AP) .OR. ( A 12) .NE.AL ) .OR.( A( 3 ) .NE.AO) .ORi I At 4) .NE. ATI ) 
IGG TO 31 
KFLAGI22I * 1 
CALL LEFT IN) 

24 IF I IAIN) .NE.AFJ.OR. IAINM) .NE.AAJ.OR.I A ( N* 2 I «NE • AR J .OR • I A I N* 3) . NE 
1 • AF ) ) GO TO 25 

IGA1N • I 
NFFP • I 
GO TO 2 7 

25 IF II AIN) .NE.AB) .OR. I AINU) .NE. A! ) .OR. I AINU) .NE. AS) .OR. I A (NO) .NE 
1 .ATM GO TO 26 

(BISC = ) 

NBIP - 1 
GO TO 27 

26 IF ( (AIN). NE.AB) . OR. I A I NH ) .NE.AAJ.OR.I A (Nf2).NE.AC) .OR.(A(Nf3)*NE 
I • AK ) ) GO TO 7) 

(SCAT « ) 

NBAP • ( 



27 DO 2B I -N, BO 
K * 4 ♦ I 
IF (All) 

20 CONTINUE 



? F * | A ( | I.EQ. SLANT) GO TO 29 



GO TO 7 I 
29 N • K 

) F (I AIN) .NE.AT) .OR.IAINM ) .NE.AHI.QR.I A 1 N* 2 ) «NE • *E ) .OR. I AINOI.NE 



193 
I 94 

195 

196 

197 
I 9B 

1 99 
200 
201 
202 

203 

204 

205 

206 
207 
20B 

209 

210 
211 
2(2 

213 

214 

215 

216 
217 

2 l B 

219 

220 
221 
222 
223 
2?4 

225 

226 
227 
22 B 

229 

230 

231 
2 32 

233 

234 

235 

236 

237 
23B 

239 

240 



) . ATM GO TO TO 
CALI EQUAL IN) 

CALL NUMBER <N,N2*X1,IX) 

IF (NFFP.EQ.1J AFFMX1 
IF INBIP.tQ.il A 3 I T r X 1 
If (NBAP.EO.l) ABATrXl 
IF IA(N2) .f Q.RIGHI ) GO TO 4 
IF IAIN2I.NE. SLANT) GO TO 71 
N ■ N2M 

GO TO 24 , 

30 IF IIA(N).NE.AP) .0R.(A(N+l).Nr.AH).0R.(A(N+2).NE.A) )) GO TO 71 
CALL EQUAL INI 

CALL NUMBER IN,N2.X| t (X) 

IF INFFP.EQ.I) afFp-xi 
IF INBI P.EU. I I AB IP-X I 
IF (NBAP.EQ.l) ABA P»X I 
IF I A1N2) .EU.R 1GHT) GO TO 4 
If I A | N 2 ) .NE. SLANT ) GO TO 7 1 
N « N2* ( 

GO TO 24 

OUTPUT 

31 IF I I A1 1) .NE .AO) .OR. ( AI2) .NE.AU) .OR. I A I 3) .NE. AT I .OR. I AI4I .NE . API I 
1 GO TO 44 

KFLAGI22) - 1 
CALL LEf T IN) 

3 2 IF II Al Nl .NE.A0 | ,0R.( AlNf) ) .NE.AI I.OR.I A(N»2) .NE.ASI .0R.IAINA3) .NE 
I . AT) ) GO TO 33 
KFLAGUB) • 1 
I BISC = 1 
CALL ECUAL IN) 

CALL NUMBER IN # N2*Xl , (X) 

PH S I • XI 
N - N2» I 

CALL NUMBER I N t N2 r X I « (X ) 

PH SF - XI 
N - N 2» 1 

CALL NUMBER I N » N 2 » X I « (X) 

THSI • XI 
N - N2» I 

CALL NUMBER I N » N2 » X I » IX) 

TH SF • XI 

IF (A|N2).EQ. RIGHT) GO TO 4 
IF I A( N2) .NE. SLANT I GO TO 71 
N - N 2* 1 

GO TO 32 



241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 
2 5B 

259 

260 
261 
262 

263 

264 

265 

266 
267 
26B 

269 

270 

271 

272 

273 

274 

275 
2 7e> 
277 
2 7 B 
279 
2 BO 
2BI 
262 
2 B 3 
2B4 
2B5 
2B6 
2 B 7 
2BB 



76 



c 



c 



c 

c 

c 



33 1 !af! ) A G0 , T0 E 3* F> • 0R * <a<n# ' 11 *NE.AA| .or. «A(N^2> .NE. AR ) .OR . (A I NO » . NE 
KF L AG I 16> - l 
I GAIN = | 

CALL EQUAL « N I 

CALL NUMBER IN, N2, XI, IX) 

PHA 1 = X 1 
N = N2» 1 

CALL NUMBER IN.N2.X1, IX) 

PHAF - XI 
N ■ N2f I 

CALL NUMBER lN,N2,X|,|X> 

T HA I » X 1 
N = UZ * 1 

CALL NUMBER IN, N2, XI, IX) 

IF ( A I N 2 ) • EQ. RIGHT ) GO TO 4 
IF UIN2) .NE. SLANT) GO TO 71 
N ■ N 2* I 
GO TO 32 

34 | IaR) ) A GO* TO E 40 N> ,0R * U<NM * -NE.AEl.OR. I AIN+2) .NE . A A > .OR. I A IN Ol .NE 
KFLAGU9) - l 
INEAR * 2 
CALL EQUAL IN) 

IF (A(N).EQ.PLEFT) GO TO 35 
INEAR - | 

I * I 
GO TO 36 



35 00 37 LM ,50 
I - L 

N = N ♦ | 

36 CALL NUMBER IN,N2,X|,|X> 
XNPI1) - XI 

N - U2* | 

C^NUMBER I N t N2 t X i t IX) 

N - N2» 1 

2NPI 1? U * 8 X I <N ’ N2fXl ’ m 
IF I1NEAR.E0.1) GO TO 39 
1 NE AR - L* 1 

IF IAIN2) .EQ. RIGHT ) GO TO 38 



289 

290 

291 

29 2 

293 

294 

295 

296 
29/ 

298 

299 

300 

301 

302 

303 

304 

305 

306 

30 7 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 
32B 

329 

330 

331 

332 

333 

334 

335 

336 



C 

C 

C 



C 



I 



C 



N - N2 
37 CONTINUE 



GO TO 71 

38 N2 - N2 ♦ 1 

I NEAR = INEAR-1 

39 IF I AIN2) .EQ. RIGHT) GO TO 4 
IF IAIN2) .NE. SLANT) GO TO 71 
N ■ N 2* 1 

GO TC 32 f 

40 IF IIAI N).NE.A8) .OR. I A( NM ) .NE.AA) .OR.IAIN*2).NE.AC).OR.lAINF3) .NE 
I . AK) ) GO TO 41 

KF l AG I 1 7 ) - | 

1 SC AT = 1 
CALL EQUAL IN) 

CALL NUMBER I N » N2 * X l , IX) 

PHU - XI 
N * N2 ♦ 1 

CALL NUMBER I N , N2 , X 1 , IX) 

PH 1 F - XI 
N = N 2 ♦ 1 

CALL NUM8ER IN.N2.Xl, IX) 

TH1I - XI 
N - N2 ♦ 1 

CALL NUMBER (N.N2.XI . IX) 

• TH1F « XI 

IF I A I N 2 ) .EQ.R 1GHT) GO TO 4 
IF IA( N2) .NE. SLANT ) GO TO 71 
N - N 2 ♦ 1 
GO TO 32 

41 IF I IAIN) .NE.AC) .OR.IAINM 1 .NE.AUI.OR.I A I Na 2 » .NE . AR ) .OR . I A I N ♦ 3 ) . N E 
I.ARil GO TO 43 

KFLAGI 15) - I 
I WR * 1 



00 42 K-N.80 

IF I A I K ) • t Q • R 1GH T ) GO TO 4 
N - KM 

IF IAIKI.EQ. SLANT) GO TO 32 
42 CONTINUE 

GO TO 71 



337 
3 38 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 
? 5 7 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 
3 70 

371 

372 

373 

374 

375 
3 76 

377 

378 

379 
3 BO 

381 

382 

383 
3 84 



77 



non 



4 3 ) F I I Af N) .NE.ASI .0R.)A)N*1) .NE.AT).OR.) AIN*2).NE.AE) .OR.) A)N»3).NE 
1.APII GO TO 71 
CALL EQUAL IN) 

CALL NUMBER 1N,N2,X1 , IX) 

SIEP - XI 

IF )A(N2).EQ. RIGHT) GO TO 4 
IF ) A ) N 2 ) • NE .SLANT ) GO TO 71 
N * N2 ♦ I 

GO TO 32 

FEEO POINT 

44 IF (I A( O.NE.AF) .OR. I A( 2) .NE.AE ) .OR. ) Al 3 ) .NE.AE ) .OR. ) A ) 4 ) .NE .AO) ) 

1 GO TO 45 
KFLAGII3I • I 
GO TO 46 

4 5 ) F I ( All) .NE.AG) .OR. ( A I 2) .NE.AE ) .OR.) A) 3).NE.AN) .0R.IAI4 ).NE .AE) I 
) GO TO 49 
KFLAG)2J) * ) 

46 NGEN * 0 
CALL LEFT IN) 

47 CALL NUMBER IN, N2, XI, 1X1 
NGEN - NGEN* 1 
KGENINGEN) = IX 

IF (AIN21 .EQ.R)GHT 1 GO TO 4 
N = N2* l 

CALL NUMBER IN,N2,X) , IX) 

V MAG * XI 
N - N2* 1 

CALL NUMBER IN, N2, XI, IX) 

VOEG = XI 

VREAL * VMAG*COSI VOEG/RAO) 

VIMAG = VMAG*SINI VOEG/RAO ) 

VOLTINGEN) * CMPLX (VREAL i VI MAG) 

IF (AIN2I.EQ. RIGHT) GO TO 4 
) F I AlN2I .NE. SLANT) GO TO 7) 

IF I I A( N2I .EQ.SL ANT) .ANO. I AIN2M) .EQ.BL ANK) ) GO TO 48 
N • N2* 1 ’ 

GO TO 47 

4 8 RE AO 15,76) A 
ICARO • 1C ARO ♦ ) 

WRITE (6,77) ) C AR 0 , A 
N « I 

CALL BLNK (Al 
GO TO 47 



3 85 

386 

387 

388 
309 
390 
39| 

392 

393 

394 

395 

396 

397 
390 

399 

400 

401 

402 

403 

404 

405 

406 

407 
400 
409 

4 10 

411 

412 

413 
4|4 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 
*31 
432 



0 E SCR 1 P T I ON 

49 IF ( ( A( 1 ) .NE .AO! . OR . 1 A ) 2 ) . NE . AE 1 . OR. ( A I 3 I . N E. AS I . OR . ( A I 4 I .NE . AC I ) 

IGO TO 52 
KFLAGI12I - ) 

J * 0 

CALL LEFT IN) 

50 CALL NUMBER (N,N2,XI,1X) 

J - J ♦ 1 

NH = J 

1 A ( J I * IX 
N = N2* 1 

CALL NUMBER )N,N2,X|,tX) 

IBIJ) = IX 

. IF (A(N2).EQ. RIGHT) GO TO 4 
IF ( A(N2) .NE .SLANT ) GO TO 71 

IF )(A|N2).EQ. SLANT). ANO. 1A)N*)).EQ. BLANK)) GO TO 51 
N - N2* 1 
GC TO 50 

5) R E AO (5,761 A 
ICARO - ICAROM 
CALL BLNK (Al 
WRITE 16,771 ICARO, A 
N = I 
GO TO 5 0 

G ECME TRY 

52 IF ( ( A( 1 I .NE .AG) .OP. ) A( 21 .NE. AE ) .OR. ( Al 3) .NE.AO) .OR. IAI4 ) • NE . AHA ) ) 

1 GO TO 55 
KF L AG ( I 2 ) - | 

JJ * 0 

CALL LEFT IN) 

53 CALL NUMBER )N,N2,Xt,tX) 

JJ - JJ*1 

NP * JJ 
X ( J J ) - X) 

N - N2*l 

CALL NUMBER IN,N2,Xl , ) X I 
Y(JJ) * XI 
N « N2* I 

CALL NUMBER (N,N2,X| , )X) 



2 ( J J ) - X) 



IF I AIN2) .EQ. RIGHT) GO TO 4 
IF )A|n2) .NE. SLANT) GO TO 7) 

IF ((A(N2I. EQ. SLANT). ANO. 1A1NO l.EQ. BLANK!) GO TO 54 



N * N2* I 

GO TO 53 



433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 
44 8 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 
4 69 
4 70 

471 

472 

473 

474 
4 75 
4 76 

477 

478 
4 79 
480 



78 



non o o o no non non non non 



54 READ 15,761 A 
1CAR0 * 1CAR0U 
WRITE 16,771 I CARO , A 
CALL 8LNK C A I 

N * 1 
GO TO 53 

INTERVAL FOR CALCULATION 

55 IF I I A( 1) .NE.AI I .OR. I Al 2 1 .N 6.ANI .OR. I At 3 I .N E. AT I .OR . I A ( 4 | .NE . AE II 
1 GO TO 56 

KFLAGI21I - 1 
CALL LEFT ( N » 

CALL NUMBER (N,N2,X1,1X| 

1 NT * IX 

IF IAIN2I.EQ. RIGHT) GO TO 4 
GO TO 71 

GROUND 

56 IF I I All) .NE.AGI .OR. I AC 2 1 .NE.A R| .OR. I Al 3 I .NE.AOI 
1GU TO 66 

KF LAG 12 51 « 1 
KFLAGI26I - 1 
1GRD * 2 
CALL LEFT INI 

57 IF I (AIN) .NE.API .OR.lAlNMI .NE.AEI .OR.IAIN+21 .NE 
' 1 • AF I I GO TO 58 

1GR0 * 1 
GO TO 64 

58 IF (IAIN). NE.AGI. OR. IAINMI. NE.AOI. OR. (AIN+2I.NE 
l.AOll GO TO 59 

ER4 * 30. 

SIG4 * .02 
GO TO 64 

59 IF I (AIN) .NE.API .OR.lAlNMI .NE.AOI. OR. I A1N»2) 

I • AR » | GO TO 60 

ER4 * 4. 

S1G4 = .001 
GO TO 64 

60 IF l(AIN| .NE.ASI .OR.lAlNMI .NE.AEI. OR. I AIN42I 
ER4 • 80. 

SIG4 - 4. 

GO 10 64 

61 IF ( (A(N).NE.AH) .OR.IA(NM) .NE.AE I.OR.I AIN+2I.NE 
l.AGII GO TO 62 

CALL EQUAL INI 
CALL NUMBER (N,N2,X1,1X| 



0R.IAl4l.NE .AUI I 

AR| .0R.(A(N+3I .NE 
. AO ) .OR .( A(Nf 3I.NE 

NE.AOI .OR.IAI N+3I.NE 



.NE . 



A A I I GO TO 61 
Al I.OR.I AINf 31. NE 



481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

51 7 

518 

519 

520 

521 

522 

523 

524 

525 

526 

52 7 
528 



HGT * XI 

IF (AIN2I.EQ.RIGHTI GO TO 4 
IF (AIN2I.NE.SLANTI GO TO 71 
N * N2» 1 
GO TO 57 

6 2 IF I I AIN) .NE.AC1 .OR.IAINM I .NE.AOl.OR.l A|N*2| .NE.ANI .OR. IAINOI.NE 
l.AOII GO TO 63 
CALL EQUAL (Nl 
CALL NUMBER IN, N2, XI, 1X1 
S1G4 * XI 

IF I AIN2) .EQ. RIGHT I GO TO 4 
IF I A I N2) .NE • SLANT ) GO TO 71 
N « N2M 
GO TO 57 

6 3 IF ( ( Al KM .NE.AOI .OR. I A I N* 1 I . NE . A 1 I .OR . I A I N«- 2 1 .NE . AE I .OR . I A I N* 31 . NE 
l.ALI I GO TO 71 
CALL EQUAL INI 
CALL NUMBER (N,N2,X1,1XI 
ER4 « XI 

IF IAIN2I .EQ. RIGHT) GO TO 4 
IF I Al N2I .NE .SLANT I GO TO 71 
N - N2M 
GO TO 57 



64 00 65 K C N, 80 

IF (AlKl.EQ.RlGHTI GO TO 4 
N * KM 

IF IAIKI.EQ. SLANT! GO TO 57 

65 CONTINUE 



GO TO 71 

66 IF l(A( II .NE.ASI. OR* I Al 21. NE .AT I .OR. I Al 3 I .NE.AOI .OR. I Al 41. NE.API I 
1 GO TO 67 

1FLAG - 2 
RETURN 

67 IF (I A I ll.NE.ACI .OR. I A I 2 1 .NE.AHI .OR.IAI 31 .N E . AA I . OR . I A I 4 I . NE . ANI I 
1 GO TO 68 

1FLAG - 3 
RETURN 

68 IF (I Al 1 I. NE.AEI .OK.I Al 21. NE. AN) .OR.IAI 3). NE.AOI I GOTO 71 



529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 
5 \2 

543 

544 

545 

546 
54 7 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 
5 70 

571 

572 

573 
5 74 

in 



1 FLAG ' 1 
RETURN 

69 1FLAG » 5 
RETURN 

70 (FLAG - 4 
RETURN 

71 MSG » 1 

KF LAG I 3 0 I - 1CAR0 

72 ^ F M FL AG .NE .51 WRITE 16,781 
1 F LAG ' 5 

RETURN 

73 1 FLAG » 6 



5X.80A1 | 





l C ar o - 
RETURN 


74 


FORMAT 1 


75 


FORMAT ( 


76 


FORMAT j 


77 


FORMAT ( 


78 


FORMAT I 
END 



/ / 7 / 5X , • 0 AT A CAROS*//) 

""A 1 | 

,l2 i»lil°ENO CARO/STOP CARO MISSING ♦*••* I 

79 



577 
5 78 

579 

580 

581 

582 

583 
5H4 

585 

586 

587 
5 88 

589 

590 

591 

592 
59 3 
59 4 

595 

596 

597 

598 

599 

600 
601 



RITE 



PURPOSE: to generate a list of branch currents from 

the input loop currents. 

METHOD: The generation of branch currents is 

accomplished in the DO LOOP ending at statement 2. The 
branch currents are stored in CJ (I) by the latter part of 
the DO LOOP ending at statement 3. If the branch currents 
are reguested for output (IWR positive) , the DO LOOP ending 
at statement 5 accomplishes this. 

CALLED BY: GANT1 

• GFFLD 

CALLS TO: NONE 



80 



o oo no no on 



SUBROUTINE RITE ( I A , I 0 , 1 NM, IWR » 1 1 » 1 2 , I 3 » HO* NOt NM t C J » CG • I GRO l 
COMPLEX C J ( l ) ,CG ( l ) ■ C J&,C JB 

0 I ME NS I ON lAlIl, IBU), 11(1), 1211), I 3 1 1 1 . M0(1NM,41, N01I1 

AMAX - .0 



00 3 K = I ( NM 
KA = I A I K } 

KB - IBIK) 

C J A * ( .0, .0 

CJB * I .0,.0 
NOK * ND( K } 



I 



00 2 I I 'I ,N0K 

1 - M0( K , I I } 

FI * I. 

IF (KB. FQ. 12(1)1 GO TO I 
IF (K8.E0.lIin) FI— I. 
CJA = CJA*FI*CJ( 1 ) 

GO TO 2 

1 IF (KA. EO. I3( I) ) FI — 1. 
CJB - CJB*FI*CJ( 1 ) 

2 CONTINUE 



CG ( K ) - CJA 
KK = K ♦ NM 
CG(KK) = CJB 
AC J * CABS(CJA) 

BCJ = CABS(CJB) 

IF (ACJ.GT.AMAX) AMAX-ACJ 
IF (BCJ.GT.AMAX) AMAX-BCJ 
3 CONTINUE 



IF ( IWR.GT.O) GO TO 4 
RETURN 

4 IF (AMAX.LE.O.) A MAX * 1 • 
WRITE (6, B) 

NMG * NM 

IF ( IGRO.GT.O) NMG - NM/2 

00 5 K= I , NMG 
CJA - CG(K) 

KK ■ K *■ NM 
CJB - CG(KK) 

CCJA • CABSICJA) 



1 

2 

3 

4 

5 

6 
7 
B 
9 

10 

11 

12 

13 

14 

15 

1 ? 

IB 

19 

20 
21 
22 

n 

25 

26 
27 
2B 

29 

30 

31 

32 

33 

34 

35 

36 

37 
3B 

39 

40 

41 

42 

43 

44 

45 

46 

47 

40 



CCJB - CABSICJB) 
AC J * CCJA/AMAX 
BCJ * CCJB/AMAX 
PA * .0 
PB * .0 

IF (ACJ.GT.O.) PA 



5 7.29 57B*ATAN2IAIHAG(CJA) f REAL(CJA) ) 



IF (BCJ.GT.O.) PB - 57.295?B*ATAN2(AIMAG(CJB» ,REAL(CJB) ) 

5 WRITE (6,7) K,1A(K), CJA, CCJA, ACJ, PA, IB(K),CJB, CCJB, BCJ, PB 

WRITE (6,6) 

RETURN 

6 FORMAT ( I HO ) 

7 FORMAT (2X,I2,2<2X,I2,2X,EI 1.5, 1X,EI1.5, IX, E 1 1 . 5 , 1 X , E 1 1 . 5, 1X,F6.1) 

0 FORMAT (/2(46X. , NCRMALlZE0 f ,5X )/ f SEG , ,2(» NOOE * , 4X, ’REAL* ,6X,*1MA 
1GINARY* ,3X,*MAiNl TUOE * , 3X , * MAGN 1 TUO E * , 3X , * PHASE • ) 1 
ENO 



49 

50 

51 

52 

53 

54 

55 

56 

57 
5B 

59 

60 
61 
62 

63 

64 

65 

66 
67 
60 
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SART 



PURPOSE 



to sort data for polar plot 



METHOD 



This subroutine sorts the values of the 



points to be plotted by the polar plot package starting with 

the greatest positive value of y to the greatest negative 

value. In the DO LOOP ending at statement 1, the value of 

(x ,y ) is interchanged with the value of (x ,y ) if y is 
i i 3 3 3 

greater than y . 

i 

CALLED BY: POLPRT 

CALLS TO: NONE 



JWHUUlINt 11/ A! A A . UA I AT i N I 

DIMENSION DA T AX ( 5001 . OAT AY I 500 ) 

THIS ROUTINE SORTS OATA IN OATAY BY MAGN 1 TU 0 E 

NN - N- l 

00 2 I * I » NN 
NM • | + I 

nn i i - um iu 



SUBROUTINE SART ( DA T A X.0A T A Y . N ) 



c 

c 

c 



3 

A 

5 

6 
T 



C 



8 



C 




9 

10 
I I 



C 

C 



2 CONTINUE 



RETURN 

EN0 
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SGANT 



PURPOSE: to calculate the mutual impedance between 

filamentary monopoles. 




1 

2 



SOURCE 

DIPOLE 



RECEIVING 

DIPOLE 



Z = P11 + P12 + P21 + P22 



METHOD: In the induced emf formulation, the mutual 

impedance of coupled dipoles is 

Z = - J I 2 ( t) E t (t) dt 

where I^(t) denotes the current distribution (normalized to 

unit terminal current) on dipole 2, and E (t) is the field 

1 

of dipole 1 when it transmits with unit terminal current. 
Distance along the axis of dipole 2 is denoted by the 
coordinate t. E^ may be expressed as the sum of the fields 

from each of the monopoles comprising dipole 1. Furthermore, 
the integral is the sum of the integrations over each of the 
monopoles comprising dipole 2. Thus, the dipole-dipole 
mutual impedance may be expressed as the sum of four 
monopole-monopole impedances. 

It may be convenient to draw the above figure in terms 
of monopoles with the current distribution shown as dotted 
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lines. (The monopole letters remain the same.) 



\ / 4 

\ P11 / 


c 


A 


\ PI 2 \ 

1 \ 

l \ 

r 


J j 






’ 1 /' 
; P2i i 


c 


B 


; P 22 > 


• i 

' i 

S l 






\ 

V >1 



The surface impedance is calculated just above statement 

2. B01 denotes J / J where J and J are the Bessel 

0 1 0 1 

functions of order zero and one with complex argument, ZARG. 
It is assumed that all the wire segments have the same 
radius, conductivity and surface impedance. 

In the DO LOOP ending with statement 3, SGANT calculates 
the segment lengths D(J). DMIN and DMAX denote the lengths 
of the shortest and longest segments. If the wire radius or 
the segment lengths are clearly beyond the range of 
thin-wire theory, N is set to zero at statement 4 followed 
by RETURN to the main program to abort the calculation. 

At statement 5, the program selects a segment K, and a 
few statements below this it selects another segment L. K is 
a segment of test diple I, and L is a segment of expansion 
mode J. The mutual impedance between segments K and L is 
obtained by calling subroutine GGS or GGMM . In statement 18, 
this impedance is lumped into C (WMM) . The mutual impedance 

Z between dipoles I and J is the sum of four 

ij 

segment-segment impedances. 

The variables IFLAG and JFLAG are used if a ground 
plane is present for the calculation of the mutual impedance 
elements. If IFLAG is equal to JFLAG, the mutual impedance 
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terms will not have the effects of a ground plane since both 
monopoles lie on the same side of the ground interface. If 
the monopoles are on the opposite sides of the interface 
(IFLAG not equal to JFLAG) , the reflection coefficient 
correction must be applied to the mutual impedance elements. 
This same technique is applied in subroutines GNFLD and 
GFFLD. 

In SGANT, segment K has endpoints KA and KB, and segment 
L has endpoints LA and LB. It is convenient to think of KA 
and KB as points 1 and 2 on segment K, and LA and LB as 
points 1 and 2 on L. The four segment-segment impedances 
can be defined as P(IS,JS). The first subscript IS refers 
to the terminal point on segment K, and the second subscript 
JS refers to the terminal point on L. Thus IS=1 or 2 if 
dipole I has its terminal point 12(1) at KA (point 1) or KB 
(point 2) , respectively. Similarly, JS=1 or 2 if mode J has 
its terminal point 12 (J) at LA or LB. The impedances 
P (IS, JS) are defined with the following reference directions 
for current flow: from point 1 toward point 2 on each 
segment. If dipole I has this same reference direction on 
segment K, FI=1; otherwise FI=-1. Similarly FJ=1 or -1 in 
accordance with the reference direction for mode J on 
segment L. In statement 18, P(IS,JS) is multiplied by FI and 
FJ before its contribution is added to Z 

ij 

Subroutine GGMM calculates the impedances Q (KK, LL) which 
are like the P(IS,JS) but have different conventions for 
reference directions and subscript meaning. The 
transformation from the Q impedances to the P impedances is 
accomplished in the DO LOOP ending with statement 13. 

If the wire has finite conductivity, the appropriate 
modification is applied to the impedance matrix just above 
statement 15. The terms arising from the dielectric shell 
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on an insulated segment are obtained from subroutine DSHELL 
just above statement 16. Finally, the lumped loads, ZLD, are 
added to the diagonal elements of the impedance matrix in 
the DO LOOP ending at statement 23. 

K is a segment of test dipole I, and L is a segment of 
expansion mode J. When the segment numbers K and L are 
egual, SGANT calls GGMM to obtain the mutual impedance 
between two filamentary electric monopoles. These monopoles 
are parallel and have the same length. Monopole K is 
positioned on the axis of the wire segment, and monopole L 
is on the surface of the same wire segment. Thus, the 
displacement is equal to the wire radius. The two monopoles 
are side-by-side with no stagger. 

When segments K and L intersect, SGANT again calls GGMM 
for the mutual impedance between the two filamentary 
monopoles. Monopole K is situated on the axis of wire 
segment K, and monopole L is on the surface of wire segment 
L. The axes of segments K and L define a plane P,and 
monopole K lies in this plane. Monopole L is parallel with 
plane P and is displaced from it by a distance egual to the 
wire radius. 

CALLED BY: MAIN 

CALLS TO: CBES 

DSHELL 

GGMM 

GGS 
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SUBROUTINE SGANT ( 1 A , 1 8 . 1 NM , 1 , ( SC , ( 1 , I 2 » 1 3 , JA , J 8 , MD , N , NO . NM, NP , A 
IH .BM.C.CGD.CHM, D,EP2,EP3,6TA # FHZ ,GAH .S&O.X, Y, Z, ZLO.ZS.ERR, iGRO) 
cchplex ERR 

COMPLEX ZG.ZH.ZS.EGD.GD.CGDS.SGOS.SGOT.8D1 

CUHPLEX PIi.Pl2,P21,P2?,QlI,Q12,Q21,Q22,EP2,EP,ETA, GAM, EP3 
COMPLEX EPStLA.CWEA.8ETA.ZARG 

CDHPl EX P(2,2),Q(2.2 ).Cg6<1 ) , SGDM ) ,C(1 ) , ZLDM ) 

0 I ME NS I DN XM). Y t I I . ZM), 0(1). 1AM). 18(1). Mn(lNM,4) 

DIMENSION mil, 12(1). (ill). JA(1), J B ( 1 I « NO I 1 T . I SC < 11 
DATA E0,TP,U0/8.854E-12,6.283l8, I .25666 -6/ 

EP ■ E P 3 

ICC - ( N* N *N ) / 2 

DO 1 1*1, ICC 
I CM) - ( .D, .0) 

ZS * I . D, .D) 

IF (CMM.LE.O. ) GD TO 2 
OMEGA = T P *F HZ 

EPSILA - CHPLXIED.-CMHM.E6/0NEGA) 

C WE A ■ ( .0,1. )*0HEGA*EPS1LA 

BETA - 0MEGA*SQRT (UD) *CSQRTI EPS1LA-EP) 

ZAPG - 8E T A • AH 
CALL C8FS (ZARG.80I ) 

ZS ■ 8f TA*BOI /CWEA 
l ZH - ZS/ ( TP*AH*GAM ) 

UHIN - 1.E3D 
DMAX - .D 

00 3 J- 1 ,NH 
K - I AUl 
L - I8( J ) 

DIJ) - SQRT ( (X(K)-X(L) ) **2 ♦ ( Y( K )- Y( L) ) • * 2 ♦( Z( K ) - Z ( L ) )** 2 ) 

IF (DU1.LT.DHIN) OMIN-DU) 

IF ( D( J ) • GT. DMAX ) DHAX*D(J) 

EGD * CEXPIGAM*D( J ) ) 

CGD(J) - ( EGDM ./EGD)/2. 

) SGD( J ) • ( EGD- 1 • / EGD )/ 2 • 

IF (DMIN.LT.2.*AM) GO TO 6 
(F (CA8SI GAM*AM) .GT.0.D6) GO TO A 
IF ( C ABS ( GAM* OH AX ).GT.3.) GO TO 4 
IF (AM.GT.O.) GO TO 5 
I CONTINUE 
N =0 

WRITE (6.24) AM , OMAX , OM IN 
WRITE (6,25) 



1 

2 

3 

4 

5 

6 

7 

8 
9 

ID 
1 1 
12 

13 

14 

15 

16 

17 

18 
19 
2D 
21 

I! 

24 

25 

26 
27 

i? 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 
4D 

41 

42 

43 

44 

45 

46 

47 

48 



OD 19 K-l.NM 
1 FLAG ' D 

IF ( ( IGRD.GT.D) . ANO. (K.GT.NM/2) ) 1FLAG* 1 
NDK ■ N 0 ( K ) 



K A 
K 8 
DK 

CGDS 

SGDS 



I AIM 
( 8 ( K ) 

D( K) 

■ CGD(K) 

■ SGD(K) 



DO 19 LM.NM 
JFLAG = D 

IF (( (GRD.GT.D). ANO. (L.GT .NM/2) ) JFLAG* 1 
NUL ■ ND ( L ) 

LA - 1 A ( L ) 

L 8 - (0U ) 

DL ■ DILI 
SGDT ■ SGD(L) 

N(L - D 

00 19 ( IM .NDK 

1 - MDI K. ll) 

MM - ( I -l )*N-( 1* 1-1 )/2 
FI * 1. 

IF ( K8. EQ • 12 ( I)) GO TO 6 
IF (K8. EQ. 11 M) ) FI — 1 . 

IS ■ 1 
GO TO 7 

IF ( K A . EQ . 1 3 ( 1)1 FI— l. 

IS - 2 

DO 19 JJ-l.NOL 
J * MD(L.JJ) 

MMM ■ MHfJ 

IF M .GT. J) GO TO 19 
F J ■ 1 . 

IF (L8.EQ.12U) 

IF (L8.EQ.IKJ) 

JS ■ 1 
GD TO 9 

IF (LA. EQ . 13 ( J ) ) FJ--1. 

JS • 2 

IF (N1L.NE.0) GO TO 18 
N (L - 1 

IF (K.EQ.U GD TO 14 

1 ND ■ ( LA-KA)MLB-KA)M LA-K8 )*( LB-K8) 

NGRO ■ 1GRD 

IF (1FL AG. EQ. JFLAG) IGRO — l 



GO TO 8 
FJ—1. 



49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
6D 
61 
62 

63 

64 

65 

66 

67 

68 
69 
7D 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 
89 
9D 

91 

92 

93 

94 
93 
96 
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r>r> 



IF ( INO.EQ.OI GO TO 10 
SEGMENTS K AND L SHARE NO POINTS 

CALL GGS (X(KA(,Y(KA!,Z(KA(|X(KBI,Y<K0(,2(KBI.X(lA|,Y(LA|*Z(LA|fXl 
1LBI ,V(LB) , Z ( L B) • AM. OK, CGOS, SGOS • OL t SGOT , INT , ETA,GAM , P( l , I I , P( 1 1 2 I . 
2P(2.1 ( ,P(2,2 I«ERR, IgRD) 

1GF0 *= NGRO 
GO TO (B 

SEGMENTS K AND l SHARE ONE POINT (THEY INTERSECT! 

10 KG - 0 

JM « KB 
JC « KA 
KF - 1 

(NO * ( KB-IA )• ( KB-L8 ) 

(F ( 1N0.NE .01 GO TO II 
JC « KB 
KF « -1 
JH * KA 
KG - 3 

11 LG * 3 
JP - LA 
LF - -1 

IF (LB. EQ. JCI GO TO 12 
JP - LB 
LF r 1 
LG * 0 

12 SGN • KF*LF 

CPS( - K X( JP(-X( JCI IMXl JN»-X( JCI I*1V( JP)-Y( JCI l*( Y(JM|-Y( JC( >♦( 2 
1 ( JP)-Z (JCI (♦( a JMl-Zl JC(I 1/ (OK*Ol) 

CALL GGMM ( .O.DK, .0,0L .AM.CGOS, SGDS.SGDT.CPSl ,ETA,GAM,Q( I. I I ,Q( 1.2 
11.0(2.1 I .0(2.21 ) 
lF( (GRD.GT .Ol SGN--SGN 

DO 13 K K = 1 ,2 
KP « (ABS(KK-KGI 

00 13 LLM.2 

LP * I A BS ( LL -LG I 

P ( KP • L P I • SGN*Q1KK,LLI 

13 CONTINUE 

(GRD-NGRD % 

GO TO 18 % 

K-L (SELF REACTION OF SEGMENT Kl 
1 * 0(1 - ( . 0 .. 0 ) 

012 » ( . 0. .01 
IF (CMH.Lt.O.I GO TO 15 
GD « GAH*OK 
ZG « ZH/ ( SGOS* • 2 I 



97 
9B 
99 
1 00 
1 0 ( 
102 
103 
10 * 
(05 
(06 
107 
I OB 
109 
1(0 
( 1 1 
1(2 
I 13 
l (* 
1(5 
I 16 
1(7 
(IB 
I 19 
1 2D 
121 
122 
123 
12* 
1 25 
1 26 
127 
1 2 B 

129 

130 

131 

132 
1 33 
13* 

135 

136 

\n 

139 
) *D 

1*1 

1*2 

1*3 

1 ** 



01( « ZG*( SGDS*CGOS-GD (/2. 

012 * ZG*(GD*CGDS- SGD S I / 2 • 

15 1 SCK - 1SC(K ( 

PI 1 = ( . 0 , . 0 ( 

P 1 2 = ( .0, .01 

(F (1 SCK. EQ. 01 GO TO 16 
(F (bM.LE.AM( GO TO 16 

CALL DSHELL ( AM , BM , OK , CGOS , SGD S , E P2 , EP , E T A , GAM , P 1 I , P 1 2 I 

16 Q ( ( « P(l»OI( 

C ALL ’ GGMM *?!d, 3K, . D, OK, AM, CGOS, SGOS, SGOS ,1. .ETA, GAM, PI 1 , PI2.P2l.P2 

12 ) 

011 * P((*Q1I 

012 * P ( 2*012 
P ( 1 .1 I - Oil 

P ( 1 . 2 ( « 0(2 

P( 2,(1 - 012 

P ( 2 , 2 I = 01 ( 

IF (KA.NE.LA) GO TO 17 
GO TO 18 

17 P(l,l) - 
P ( 1,21 - 
P(2, II * 

P ( 2 , 2 I * 

18 C (MMM) = 

19 CONTINUE 



-012 
-Oil 
-OK 
-01 2 

C(MMM) »F)*FJ»P( 1 S, JSI 



DO 23 I r 1 , N 

MM * ((-II *N-( 1*1-11/2 

l J * MM ♦ 1 
J JA - J A | 1 | 

J1 « J J A 
( (2 - 12(11 
( (1 - 11(0 

IF ( 112.EQ.10UI I I Jl-JKNM 
JJB ■ JB( 1 I 
J 2 - JJB 

IF ( 112.EO.10( J2I 1 J2-J2*NM 
C( 1JI « C ( (JI*ZL0IJ1I*ZL0(J2I 
JJJ - J J A 

DO 22 K«l,2 
NO J « N D( J J J I 

00 21 JJ-l.NOJ 
J • M0( J J J * J J I 
IF (J.EQ.lI GO TO 2) 



1*5 

1*6 

1*7 

1*B 

1*9 

150 

151 

152 

153 
15 * 

155 

156 

157 
15B 

159 

160 
161 
162 
163 
( 6 * 

165 

166 
167 
I 6 B 
169 
I7D 
I 71 

172 

173 
) 7* 
1 75 

176 

177 
17B 
I 79 
I0D 
101 
182 
| 03 
IB* 
185 
166 
187 
IBS 
1 B9 
19D 

191 

192 
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c 

c 

c 

c 



If < 121 JI.NE.II2I GO TO 21 
I J * MM ♦ J 
FI = I. 

IF IK.E0.2) GO TO 20 
IF (III JI.NE.III) FI — |. 

C C I J1 - Cl I J ) +F I • ZLDI J 1 ) 

GO TO 21 

20 IF I 131 JJ.NE.l3t III FI — 1. 
Cl I J! * Cl I J > ♦ F I *ZLD( J2 I 

21 CONTINUE 

22 JJJ - JJB 

23 CONTINUE 
RETURN 



2* FORMAT l‘3X,*AM - • . E I 0. 3 , 3X , * DM AX ■ • .E I 0. 3 , 3X , • DM IN ■ *.EI0.3) 

25 FORMAT (• WARNING ••••••••♦•••••••••••••••••••••••••••••••••♦••••/ 

I,' THIS PROBLEM EXCEED LIMIT OF THIN WIRE CONDITION, THE RESULTS 
2 ARE NOT CORRECT* ) 

END 



193 

19* 

195 

196 

197 

98 

99 
200 
201 
202 
203 
20 * 

205 

206 
207 
200 

209 

210 

M 

213 

21 * 
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SORT 



PURPOSE: to define the set of dipole inodes. 

METHOD: In the DO LOOP ending at statement 3, the set 

of dipoles is defined by filling the vectors II (I) and 13 (I) 
(the endpoints of dipole I) ; 12 (I) ( the terminal point of 

dipole I); and the vectors JA(I) and JB(I) (the monopoles 
comprising dipole I) with the node numbers and segment 
numbers, respectively. The DO LOOP ending at statement 8 
determines MD(J,K) (the list of dipoles sharing segment J) 
and ND (K) (the number of dipoles sharing segment J) . 

CALLED BY: MAIN 

CALLS TO: NONE 
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SUBSQUT INE SORT I I A , I B , i 1 , I 2 . I 3 . J A, JB, HO . NO* NH, N P. N , MAX , M I N , I C J , 1 N 
i M 1 

DIMENSION JSPC20) 

OIHFNSION llCn. 12(1 I • 13(1). J A ( I ) « J B ( 1 1 

0 I HE NS I ON 1 A ( 1 1 « 1 B I I I . NOC 1) , MD(INH,4) 

1 - 0 

00 3 K - 1 . NP 
NJK * 0 

00 I J e l i NH 

1 NO * ( 1 A ( J )-K )* ( IB(J)-K) 

IF I 1N0.NE.0) GO TO I 
NJK * NJK * I 

JSP(NJK) ■ J 
l CONTINUE 

HOD = NJK-1 

IF ( HOD* LE .01 GO TO 3 

00 2 1 HO* 1 « HOD 

1 * hi 

IF (1.GT.1CJ) GO TO 2 
1 PO - 1 HO* 1 
J A 1 * J SPC IMO) 

JAM! * J A 1 
JB I « JSPC1PD) 

JBI11 - JB 1 
III I 1 - lA(JAl) 

IF (1A( JAD.EO.K) 

12(11 - K 
13(11 * lA(JBj) 

IF (1A( JBII.EO.K) 
l CONTINUE 



IK D-IBIJAI) 
I 3( 1 1 -I B( JBI 1 



CONTINUE 
N - 1 



00 4 J« 
ND( J) - 



'jNH 



00 4 K* 1.4 
HOU.Ki • 0 



I II - N 
IF (N.GT.ICJ) 



III 



ICJ 



1 

2 

3 

4 

5 

6 
7 
B 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IB 

19 

20 
21 
22 

23 

24 

25 

26 
27 
20 

29 

30 

31 

32 

33 

34 

35 

36 

37 
3B 

39 

40 

41 

42 

43 

44 

45 

46 

47 
4B 



00 B 1-1. I I 1 
J - J A ( I J 

00 7 L* I i 2 . . 

NO ( J 1 * NO ( J ) ♦! 

K * I 
H ■ 0 

5 HJK - HOU.KI 

IF (HJK.NE.OI GO 'TO 6 
H * 1 

HOU.KJ « I 

6 K « KM 
IF (K.GT.4) 

IF (H.EQ.O) 

7 J - JBI II 

B CONTINUE 



GO TO 7 
GO TO 5 



HIN 

MAX 



100 

0 



00 9 J ■ I , NH 
NO J • NDt J ) 



IF INOJ.GT.HAXI HAX-NOJ 
IF (NOJ.LT.MINI HIN-NOJ 

RETURN 

ENO 
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52 

53 

54 

55 

56 

57 
5 B 

59 

60 
61 
62 

63 

64 

65 

66 
67 
6B 
69 

i? 

n 

74 

75 

76 
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SQROT 



PURPOSE: to solve the set of simultaneous equations to 

determine the currents on the thin wire structure. 

METHOD: This subroutine considers the matrix equation 

ZI = V which represents a system of simultaneous linear 
equations. NEQ denotes the number of simultaneous equations 
and the size of the matrix Z. 

On entry to SQROT, S is the excitation column V. On 
exit, the solution I is stored in S. Z(I,J) denotes the 
symmetric square matrix. Also on entry, the upper-right 
triangular position of Z(I,J) is stored by rows in C(K) with 

K = (1-1) * NEQ - (I * I) / 2 + J . 

If 112 = 1, SQROT will transform the symmetric matrix 
into the auxiliary matrix (implicit inverse) , store the 
result in C (K) and use the auxiliary matrix to solve the 
simultaneous equations. If 112 = 2, this indicates that C(K) 
already contains the auxiliary matrix. 

The transfromation from the symmetric matrix to the 
auxiliary matrix is accomplished in the DO LOOP ending at 
statement 5. The solution of the simultaneous equations is 
accomplished in the remainder of the program. 

CALLED BY: GFPLD 

CALLS TO: NONE 
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on 



C 



SUBROUTINE SQROT ( C . S , I MR . 1 12 . NE Q I 
COMPLEX Cl 11. Sill, SS 
N * NE 0 

IF ( 112. EQ. 2) 30 TO 6 
C(ll - CSORMCt 1 I I 

00 1 K- 2 • N 

1 C ( K 1 - cUi/cu I 



C 

C 

c 

c 



c 

c 

c 

c 

c 

c 

c 



00 5 1 * 2 • N 

1 NO ■ l-l 
IPO - 1*1 

10 ■ 11-11 •N-l 1*1-1 1/2 

11 - 10*1 

00 2 L-I.IHO 

tl - 1L-1 I*N-IL»L-U/2M 

2 CUll - Cllll— C(Ll 1*CILU 

cuu - csoRUcmn 

IF ( 1P0.GT.N) GO TO 5 

00 4 J-1P0.N 
1J ■ 10*J 

00 3 H- 1.1 HO 
HO - 1H-1 I *N — I H* M-H 1/2 
HI - HO* 1 
HJ - HO * J 

3 Cl 1 J) - Cl 1 J 1 -C I H J ) *C 1 H I ) 

4 ciiji - ciui/cim 

5 CONTINUE 

6 sm - sm/cm 

00 8 1-2. N 

1 HO - (-1 

00 7 L >1.1 HO 

LI - Tl-I l*N*l L*L-Ll/2*1 

7 SM) - S ( 1 l-CI L 1 I * S 1 L ) 

11 - (l-ll*N-(l*MI/2U 

asm - si i »/c« m 



1 

2 

3 

4 

5 

6 
7 

e 

9 

10 

U 

II 

l* 

15 

If 

ie 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

M 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 



C 

C 

C 

C 

C 

C 

C 



C 



C 



NN - < l N*1 I*NI/2 
S1N1 - S1NI/C1NNI 
NHO - N-l 

00 10 1*1. NHO 
K = N-l 
KPO * KM 

KO - 1K-I |*N-( K*K-K) /2 

00 9 L-KPO.N 
KL - KO* L 

9 S ( K I - SIKI-C1KLMS1LI 
KK - K0*K 

10 S1K1 - S(K1/C(KK) 

IF 1 (WR.LE.OI GO TO 13 
CNOR - .0 

SSi'ciMm 

11 IF ISA. GT. CNOR) CNOR- SA 

IF 1 CNOR. L E • 0» 1 CNOR- 1 • 

00 12' 1-l.N 
SS - SI 11 
SA - CABS 1 SSI 
SNOR - SA/CNOR 

IF "si?GT.0.l PH • 57.29578*ATAN2(AlMAG( SS) .REAL (SSI) 

12 WRITE 16.141 1 .SNOA.SA.PH.SS 

WRITE (6.15) 

13 RETURN 

14 FORMAT UX.ll5,lFl0.3f 1F15.7.1F10.Q.2FI5.6I 

15 FORMAT 1 IHdl 
ENO 
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TO 
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• BLNK 








.READ • I 



.EQUAL 

? • T * * • ‘LEFT 

* * NUMBER 



. SORT 
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*: 



. 

• • 
• • 



^ . .CBES 
. • DSHELL 

• • GGS • • 




• • • GANT1 « 

• • . 

* GNFLD* 
GFFLD* 



© 



GNF 



GFF 

© 







*POLPRT« 



•POLPLOT* 

•SART 










LINECK 

NUMB 



EX P J 

GGMM 

. . GDISS 

® . • • 

• | ••• • .RITE 

• . . 

* SQROT 

CALLING SEQUENCE OF THE SUBROUTINES 
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SYMBOL DICTIONARY 



A 

ABAP 

ABAT 

ABIP 

ABIT 

ACSP 

ACST 

AFFP 

AFFT 

AM 

BM 

C 

CG 

CGD 

CJ 

CMM 

D 

ECSP 

ECST 

EFF 

EP 

EPP 

EPPS 

EPTS 

EP2 

EP3 

EP4 

ERR 

ER2 



characters of the input data cards 

backscattering phi plane angle for plotting 

backscattering theta plane angle for plotting 

bistatic scattering phi plane angle for 
plotting 

bistatic scattering theta plane angle for 
plotting 

absorption cross section for phi polarization 

absorption cross section for theta 
polarization 

far-zone phi plane angle for plotting 

far-zone theta plane angle for plotting 

radius of the thin wire of the structure 

outer radius of the dielectric shell of the 
insulation of the wire 

elements of the open-circuit impedance matrix 

branch currents for the structure 

cosh yd for a given segment 

loop currents for the structure 

conductivity of the wire 

length of a given segment 

extinction cross section for phi polarization 

extinction cross section for theta 
polarization 

radiation efficiency 

loop currents induced by a phi polarized wave 

phi-polarized far-zone feild of the dipole 
mode 

scattered electric field in the phi direction 
due to a phi polarized wave 

scattered electric field in the theta 
direction due to a phi polarized wave 

complex permittivity of insulation 

complex permittivity of ambient medium 

complex permittivity of ground 

EP4/EP3 

relative dielectric constant of insulation 
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ER3 

ER4 

ET 

ETA 

ETPS 

ETT 

ETTS 

EX 

EY 

EZ 

EO 

FHZ 

FMC 

GAM 

GG 

GPP 

GTT 

HGT 

IA 

IB 

IBISC 

ICARD 

ICJ 

IFLAG 

IGAIN 

IGRD 

INC 

INEAR 

INM 

INT 



relative dielectric constant of the ambient 
medium 

relative dielectric constant of the ground 

loop current induced by a theta polarized 
wave 

intrinic impedance of ambient medium 

scattered electric field in the phi direction 
due to a theta polarized wave 

theta polarized far-zone field of the dipole 
mode 

scattered electric field on the theta 
direction due to theta polarized wave 

near-zone electric field in x direction 

near-zone electric field in the y direction 

near-zone electric field in z direction 

8.854E-12 

frequency in hertz 
frequency in megahertz 

intrinic progration constant of the ambient 
medium 

time-average power input 

power gain associated with the phi polarized 
component 

power gain associated with the theta 
polarized component 

height of the structure above ground plane 

first node of a given segment 

second node of a given segment 

indicator for bistatic scatter calculations 

indicator for the data cards 

dimension corresponding to the number of 
simultaneous linear equations 

indicator for program termination 

indicator for antenna gain calculations 

indicator for presence of the ground plane 

indicator for the type of far-zone 
calculations 

indicator for near-zone calculations 

dimension corresponding to the number of 
mon poles 

number of integration steps 
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ISC 

ISCAT 

IWR 

11 
112 

12 
13 
JA 
JB 

KFLAG 

KGEN 

LOAD 

LZD 

MAX 

MD 

MIN 

MSG 

N 

ND 

NGEN 

NM 

NPL 

OMEGA 

PH 

SCSP 

SCST 

SGD 

SIG2 

SIG3 

SIG4 

SPPM 

SPTM 

STPM 

STTM 



indicator for the insualtion 
indicator for backscatter calculations 
indicator for current distribution output 
endpoint node of a given dipole 
indicator for auxiliary matrix 

terminal node number of a given dipole 

endpoint node number of a given dipole 

first segment number of a given dipoile 

second segment number of a given dipole 
print indicator 

list of generator/feed locations 

indicator for structure load 

list of impedance/load locations 

maximum of the number of segments connected 
to any one given node 

list of dipoles sharing a given segment 

minimum of the number of segments that 
connected to any one given node 

indicator for error printout 

number of simultaneous linear equations 

total number of dipoles sharing a given 
segment 

indicator for antenna calculations 
number of segments 
indicator for polar plot 
angular frequency 

phi angle for far-zone calculations 

scattering cross section for phi polarization 

scattering cross section for theta 
polarization 

sinh yd of a given segment 
conductivity of insulation 
conductivity of the ambient medium 
conductivity of ground 

echo area phi incident-phi scattered wave 
echo area phi incident- theta scattered wave 
echo area theta incident-phi scattered wave 
echo area theta incident-theta scattered wave 
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TD2 


loss tangent of 


the insulation 


TD3 


loss tangent of 


ambient medium 


TH 


theta angle for 


far-zone calculations 


TP 


2n (6.28318) 




00 


1 . 2566E-6 




VG 


antenna complex 


driving voltages 


VOLT 


list of VG's 




X 


x-coordinate of 


each node 


XNP 


list of XP 1 s 




XP 


x-coordinate for 


near-zone calculations 


Y 


y-coordinate of 


each node 


YNP 


list of YP 1 s 




YP 


y-coordinate for 


near-zone calculations 


Y 1 1 


complex power input 


Z 


z-coordinate of 


each node 


ZLD 


complex load at 


a given node 


ZLLD 


list of ZLD's 




ZNP 


list of ZP's 




ZP 


z-coordinate for 


near-zone calculations 


ZS 


surface impedance of the wire 


Z 1 1 


antenna input impedance 
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USER'S MANUAL 



The Antennas-Scatterers Analysis Program (ASAP) for thin 
wire structures in a homogenous conducting medium performs a 
freguency domain analysis of antennas and scatters. The 
program is applicable in the presence of either a perfect or 
a finite ground. This appendix will describe and explain 
the data cards necessary to execute the compute program. 
Although the program was written for the IBM 360 computer 
system, it can be executed on another system with minor 
modifications. 

The program utilizes piecewise sinusoidal expansion for 
tbe current distribution with Kirchhoff Current Law enforced 
everywhere on the structure. If the structure contains end 
points, the currents at these points are assumed to vanish. 

I. Program Limits 

The thin wire assumptions are questionable and the 
accuracy and convergence deteriorate if the radius of wire 
utilized for the structure exceeds 0.01 of a wavelength, if 
the longest segment is greater than one-fourth of a 
wavelength, if the length ratio of the longest and shortest 
segments exceeds 100, or if the total wire length is less 
than 30 times the wire diameter. If a wire is bent sharply 
to form a small acute angle (less than 30 degrees) , the thin 
wire model is questionable. It is assumed that the wire 
conductivity greatly exceeds the conductivity of the ambient 
medium. For insulated wires, the dielectric layer is assumed 
to be electrically thin. 

II. Minimum Data 

The minimum data necessary to execute the program is: 
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a. description of structure 



b. radius of wire used for the structure 

The program will default to the other parameters necessary. 
The default parameters are: 

a. wire for the structure is copper 

b. frequency of operation is 300 mhz 

c. homogneous medium is free space 

A more detailed explanation of the defaults will be 
discussed when the data card for the parameter is described. 

III. Outputs 

In antenna problems, the output includes structure 
currents, impedance(s) of feed(s), gain, polar radiation 
plots, and near field calculations. In bistatic scattering 
problems, the output includes structure currents, complex 
elements of the polarization scattering matrix, polar 
reradiation pattern plots, and echo areas produced by a 
plane wave. For backscatter ing problems the output includes 
absorption, scattering and extinction cross sections in 
addition to the outputs of bistatic scattering. Most of the 
outputs are suppressed and must be requested. Since the 
program can produce a large volume of output, care should be 
exercised until the user is familiar with the outputs. 

IV. Data Cards 

The Analysis Program utilizes free format for the data 
cards, that is, the program utilizes character recognition 
to determine which parameters are being read. Data placement 
(location) on the input card is not critial. Blank 
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characters, on all input cards but the COMMENT data card, 
are ignored and may be used at the discretion of the user. 
Since character recognition is used, only the first four 
characters of the key words must be present and correct. 

The format for the COMMENT CARD utilizes standard 
FORTRAN format (i. e. *C* in column 1 followed by at least 
four blanks) . The COMMENT CARD is the only type of input 
card that position in the data block is critial. This 
(these) card (s) must be placed at the beginning of a data 
block. A data block is a series of related data cards. 
Several data blocks may be used to define an analysis 
problem. This will become clear when the termination cards 
(END, STOP, or CHANGE) are discussed. There is no limit to 
the number of comment cards that may be used. As a check for 
the user, all input data cards will appear on the output as 
they appear in the input deck. 

The format of other data can be of one of two forms: 

a. type of card (option 1/option 2/ ) 

b. parameter (value) . 

The type of format to use will be apparent as the individual 
data cards are discussed. 

The numerical values for the parameters may be stated in 
any one of the following forms. The program will translate 
the number to the proper form for the specified parameter, 
either fixed or floating point. All of the following 
examples have the same value. 

0.0001 or .0001 or 100. U or 100U or .1M or 0.1M or .0000001K 

-6 -3 3 

0 = 10 M = 10 K = 10 
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1. WIRE This card is used to define the parameters 
associated with the wire utilized by the thin wire 
structure. Two options are available and are defined as: 

RADIUS=value of the radius of the wire in meters 

CONDUCTIVITY=value in megamhos per meters . 

The wire data card must appear in the first data block to 
define wire radius. The default value of the conductivity 
is 50 megamhos/meter (copper) . 



IRE ( RADI US= . 00 1/ CONDU=28.5) 



2. INSOLATI ON This card is utilized to define the parameters 
associated with the insulation of the wire used for the 
structure to be analyzed. If this card is omitted, the 
program assumes that the structure is uninsulated. Four 
options are available and are defined as: 

RADIUS=value of outer radius in meters 

CONDUCTIVITY=value in micromhos per meter 

DIELECTRIC=value of relative dielectric constant 

LOSS TANGENT= value . 

The conductivity and either the relative dielectric constant 
or the loss tangent (but not all three) options may be 
stated. 
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3. EXTERIOR MEDIUM This card is utilized to describe the 
homogeneous medium surrounding the structure. If the medium 
is free space, this card may be omitted. Three options are 
available and are defined as: 

DIELECTRIC= value of relative dielectric constant 

CONDUCTIVITY = value in micromhos per meter 

LOSS. TANGENT=value . 

As with INSULATION card state either conductivity or loss 
tangent. 
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4. DESCR I PTION This card is utilized to describe the shape 
of the wire structure to the program. The user must divide 
the wire structure into segments of the appropriate length 
and number each node starting at one. A node is a point 
where a segment begins or ends. A maximium of four segments 
can meet at any given node. An isolated wire must contain 
at least two segments and three nodes. The structure is 
described by stating the node numbers that each segment 
connects. The description of a square loop might appear as: 

DESCRIPTION (1-2/2- 3/3-4/4-1) . 

The description of a dipole and reflector might appear as: 

DESCRIPTION (1-2/2-3/3-4/4-5/6-7/7-8/8-9/9-10) . 

If the description will not fit on one data card continue on 
the next card as if the previous card were longer. The 
dipole example might appear as: 

DESCRIPTION ( 1-2/2 -3/3- 4/4- 5/ 

6-7/7-8/8-9/9-10) . 

Note that the last character on the card to be continued is 
a slant (/) . As many cards as necessary may be used. The 
maximum number of nodes permitted is fifty. If ground plane 
is present, the maximum number is twenty-five. If a ground 
plane is present and the structure touches the ground plane, 
the lowest node numbers MUST be used for the touching nodes. 
That is, if the structure touches the ground plane at two 
points, node numbers 1 and 2 MUST be assigned to these 
nodes. 
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5. GEOME T RY This card is used to state the physical location 
in rectangular coordinates of each node of the DESCRIPTION 
CARD . The rectangular grid is in units of meters. If node 1 
is located at x1,y1,z1 and node 2 at x2,y2,z2 and node 3 at 
x3, y3 ,z3 ,etc. , the GEOMETRY CARD might appear as: 

GEOMETRY (X 1 , y 1 , z 1/x2 , y2 , z 2/x3 , y 3 , z3/ ) 

As with the DESCRIPTION CARD, continuation cards are 
permitted. 



'GEOM ( . 1,0,. 1/-. 1,0. 1/-. 1,0-. 1/. 1,0,-. 1) 
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6. FEED For antenna analysis the feed point (s) and 
voltage (s) must be stated. In the foremention dipole and 
reflector example if the feeds were at node 2 with a voltage 
source of .5 at an angle of -90 degrees and at node 4 with a 
voltage source of .5 at an angle of +90 degrees the FEED 
CARD might appear as: 

< 

FEED (2,. 5, -9 0/4,. 5 ,+9 0) 

The order of the information for each voltage source is node 
number, magnitude, and phase angle. This order is repeated 
until all sources are stated. If the source information will 
not fit on one card, use another card similiar to the 
initial one; that is, repeat the word "FEED”. If only one 
voltage source is applied to the structure, only the node 
number must be stated. In the dipole example, if the drive 
is at node 3, the FEED CARD might appear as: 

FEED (3) 

A default source of one voltage at zero degree phase is 
assumed. Voltage sources should only be stated for nodes 
with only two segments. 



FEED (2,. 5, -9 0/4, .5, +90) 



7. LOAD This card is used to describe the loads to be placed 
at various locations on the structure. The format for this 
card is similiar to that of the FEED CARD, that is, the 
word "LOAD” is used in the place of "FEED". The order of 
the information on the card is the same. Since this card is 
frequency dependent, it must be changed if the frequency of 
operation is changed. No default parameters are available. 
The structure is assumed unloaded unless this card is used. 
Once the structure is loaded, it will remain loaded for the 
remainder of the data block series. To unload the structure 
the following card may be used: 

LOAD (-1) 
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8. O UTPUT This card is used to request output data. Most of 
the output is in tabular form. More than one OUTPUT CARD 
is permitted per data block, but not for the same type of 
output. If only the antenna input impedance, antenna 
efficiency, or time-average power input is of interest, no 
OUTPUT CARD is necessary. These parameters are automatically 
printed if a FEED CARD or GENERATOR CARD is utilized. One 
or more of the following options may be used to request the 
various outputs available. 

FAR FIELD=phi initial, phi final, theta initial, theta 

final 

This option gives the components of the electric field 
intensity in the far field as phi and theta varies between 
limits specified in one degree divisions. 

BACKSCATERING=phi initial, phi final, theta initial, 
theta final 

This option gives the absorption, scattering, and extinction 
cross sections, and the complex elements of the polarization 
scattering matrix for an incident plane wave illuminating 
the structure from the spherical direction of phi, theta as 
both vary between limits specified in one degree divisions. 

BISTATIC=phi inital, phi final, theta initial, theta 

final 

This option gives echo area and the complex elements of the 
polarization scattering matrix for an incident plane wave 
illuminating the structure from the spherical direction phi, 
theta final of the backseat tering output option, reradiated 
in the phi, theta direction as both vary between limits 
specified in one degree divisions. A bistatic output request 
must be accompanied with a backscattering request in the 
same data block. 

STEP=value in degrees 
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This option will cause any of the above output options to be 
stepped at a different interval size. That is, if one of the 
above options is to be stepped at ten degrees intervals, use 
this option. This option overrides the one degree stepping. 

CURRENT 

This option gives the currents on the structure which are 
produced by the feed/generator voltages and/or the incident 
plane wave of the backscattering request. 

NEAR=x1,y1,z1 

or 

NEAB=(x1,y1,z1/x2,y2,z2/x3, y3,z3/etc ) 

This option gives the value of electric field components in 
the near field for the antenna at the point or points 
specified. 



/OUTPUT (FARF = 45, 50,25, 50) 
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9. PLOT This card will produce normalized polar plots in the 
specifed plane for the stated option. The plane is 

specified by stating either ”PHI= " or »THETA= " . The 

PLOT CARD overrides the limits of the OUTPUT CARD for the 
same option. If only a normalized pattern is of interest, 
only a PLOT CARD is necessary. If a table of values and a 
normalized pattern is desired, both a PLOT CARD and OUTPUT 
CARD must be used. Only one PLOT CARD is permitted per 
data block. The following pattern plots are available: 

FAR FIELD/plane 

This option will plot the far field intensity for each 
component of the electric field. 

BACKSCATTERING/plane 

This option will plot the normalized magnitude of each of 
the elements ot the polarization scattering matrix. 

BISTATIC/plane 

This option will plot the normalized magnitude of each of 
the elements of the polarization scattering matrix produced 
by the incident plane wave stated by final limits of the 
backscattering option of the output request. 



PLOT (FARF/THET=90) 



10. G RO UND This card is used to describe the ground 
parameters if a ground plane is present. If no ground plane 
is present, the structure is assumed to be in free space or 
the homogeneous medium of the EXTERIOR MEDIUM data card. 
Seven options are available and are defined as: 

PERFECT 

This option will analyze the structure over a perfect ground 
plane. 

GOOD 

This option will analyze the structure over a good ground 
plane where the conductivity of the ground is .02 mhos/meter 
and the relative dielectric constant is 30. 

POOR 

This option will analyze the structure over a poor ground 
plane where the conductivity of the ground is .001 
mhos/meter and the relative dielectric constant is 4. 

SEA 

This option will analyze the structure over salt water where 
the conductivty of the water is 4. mhos/meter and the 
relative dielectric constant is 80. 

HEIGHT=value in meters 

This option will analyze the structure with origin of the 
GEOMETRY card this height above the ground plane. The 
lowest point of the structure must not lie below the ground 
plane. It may lie on the ground plane. 

CONDUCTIVITY = value in mhos/meter 

This option is used to state the value of conductivity of 
the ground plane if the default values mentioned above are 
not utilized. 
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DIELECTBIC= value 



This option is used to state the relative dielectric 
constant of the ground plane if the default values mentioned 
above are not utilized. 



/ground (HEIG=1 0/COND=. 002/DI EL= 10) 
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11. INTE RVAL FOR C ALCULA TION This card is used to state the 
number of intervals to be used for calculating the elements 
of the impedance matrix with Simpson ' s-rule integration. A 
large value for the number improves the accuracy at the 
expense of greater execution time. For most problems a 
suitable combination of speed and accuracy is obtained with 
a value of four, the default value. If the rigiorous 
closed-form impedance expressions in terms of the 
exponential integrals is desired, set this value to zero. 

INTERVAL=value 



INTE (6) 



12. GENERATO R This card is similiar to the FEED CARD in use, 
except that the segment numbers are stated instead of the 
node numbers. This is useful if tl\ree or four segments meet 
at a node. The positive terminal of the generator is 
connected to the specified segment such that current is 
forced in the the positive direction. The positive direction 
of current flow is from the first stated node number of that 
segment toward the second stated as ordered on the 
DESCRIPTION CARD. 



GENE (2,. 5, -9 0/4, .5, +90) 



13. I MPEDANC E This card is similiar to the LOAD CARD in use, 
except that the segment numbers are stated instead of the 
node numbers. As with the GENERATOR CARD, this is used if 
three or four segments are connected to a node. The 
impedance will be connected to the positive terminal of the 
specified segment. The format of this card is the same as 
the LOAD CARD. 
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14. CHANGE This card at the end of the data block signals 
the program that the following data cards are changes to the 
previously read data, for the next run. If a "CHANGE CARD" 
is used, the outputs must be requested again in the next 
data block. 



15. END This card signals the program that this is the end 
of a data block series and to reinitialize data for the next 
problem. An "END CARD" cannot be used with a "CHANGE CARD". 



16. STOP This card signals the program that all of the data 
cards have been read and to terminate itself when execution 
is completed. This card must be used as the last card in 
place of the "END CARD" of the last data block series. A 
"STOP CARD" cannot be used with an " END CARD" in the same 
data block. 
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C AN EXAMPLE PROBLEM 

C 

C V ANTENNA 

C 

WIRE (RADI US = 1M) 

GE0M(0,~. 18, +.18/0, -.09,*. 09/0, 0,0/0, 0.09, .09/0 , . 18,. 18) 
DESC(l-2/2-3/3-4/4-5) 

FEED( 3 ) 

OUTPUT! FARF=45, 50,65, 80/S TEP=5) 

CHANGE 

OUT PUT (BIST=45»45,45,^5/BACK=0»0» 10, 12) 

OUTPUT! CURRENT) 

CHANGE 

C 

C CHANGE STRUCTURE SHAPE TO DIPOLE 

C 

GE0M(0,-. 25, 0/0, -.125, 0/0, 0,0/0, .125, 0/0, .25,0) 
PL0T(FARF/PHI=90) 

GROUND (HEIGHT = . 2 5/ GOOD ) 

STOP 



THE ABOVE DATA DECK WILL PRODUCE THE OUTPUT ON THE 
FOLLOWING PAGES. 
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FAR FIELOS FOR PHI VARYING FROM 45.0 TO 50.0AN0 THETA VARYING FROM 65 

IN STEPS OF 5.0 OEGREES. 
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FAR FIELD 
CALCULATIONS 
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CONTINUE EXECUTION WITH THE FOLLOWING AOOITIONS AND/OR CHANGES 
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STRUCTURE CURRENTS 

BACK SCATTERING FOR PHI VARYING FROM 0.0 TO 0.0 ANO THETA VARYING FROM 10.0 TO 

IN STEPS OF 1.0 OEGREES. 

BISTATIC SCATTERING FOR PHI VARYING FROM *5.0 TO *5.0 AND THETA VARYING FROM *5 

IN STEPS OF 1.0 OEGREES. 
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CURRENTS INDUCED BY THE PHI POLARIZED WAVE 
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CHANGE STRUCTURE SHAPE TO OIPOLE 
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OUTPUT REQUESTEO 
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PHI ELECTRIC FIELO ANTENNA PATTERN FOR SPECIFIED PLANE. NORMALIZING FACTOR* .17362E-06 
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THET ELECTRIC FIELD ANTENNA PATTERN FOR SPECIFIED PLANE. NORMALIZING FACTOR* .42801E 00 



o 

o 



o 

• m 
o 

* * 

• •* # 

• ft * 

• •## 

• « 

• * • 

• ft 

* • 

* * 

* • 



ft * ft * # 
* * * 
ft ft ft 



•»»«**« • • • 



in 

• • • • »s 
• • ^ 



####» ••• 

ft * ft • • • 

ft ft ft • • 

* ft ft • • 

* # • 
* * 

« * 

♦ ft 



O 

m • 
m 



ft • 
* 
• ft 



* 

* 

* * 
* 
* 



• 




ft • • 








• • • • s • • 


• • 


• 




ft • 






* • t 


m 


• • 


• o • 




• • 




• • 




• 




•vO 




* 


• • 






• 




o • • 




ft • 


• 






• 




• • 




ft • 








• 




• • 




ft • 


• 






• 




• 




• ft • 








• 




• 


• • 


ft • 


• 






• 








• • ft • 








in 








• ft • 








• • *N, • • 


• 


• 




• • • ft 




• 


• • 


• no 


• • • 






• ft • 






• 


• 








• •# • 




• • 




• 




• 




ft 


• 






• 








ft • 


• 




• 


• 






• 


ft • • 








• 








ft • 


• 




• 


• 


• 






• ft • • • 








• 




• 




ft ft • 






• 


in 


• 




• 


ft • 








• • • 








• ft • 


• 






• • *-4 


• • 






ft 


• 




• 


• 


• 






ft ft 


• 




• 


• 


• 






ft 




• • 




• • 


• 


• 




ft ft 






• 


• 








• 


ft ft 










• 




• ' 


ft ft 






• • 


• 






• 




ft ft 




• • • • 


• 


• 




• • 






ft ft ft 


• 


ft ft* 
















• 




• • 








• 








• 


1 ' 


• 




• 1 • 1 


• 


• 




• 








• • 


• 






• 












• 










• 


• 












• • 




• • 


• 






• 


• 




• 


• 


• 








• 




• 


• 


• 






• • 


• 






• • in 


• • 




• 


• 








• IS • • 




• 




• 






• 


-4 


• 



• * 
* 

* 

• •# 

• ft 
« • 
• * 



* # 



« * 



o 

• <M 



in 

• • »s • 



in 

•s • • 

m 



in 

•s, • • • • 



• • 

• ♦ • 



O 

•in 



• • 

o 

—4 

<M 



O 

CO 



13 ? 



****** 



LIST OF REFERENCES 



2. (a) 



(b) 



3. 



4. 



5. 



Richmond, J.H., "Radiation and Scattering by 
Thin-Wire Structures in the Complex Frequency 
Domain," Report 2902-10, July, 1973, The Ohio State 
Unversity ElectroScience Labortory, Department of 
Electrical Engineering: prepared under Grant NGL 
36-008-138 for National Aeronautics and Space 
Administration, Langley Research Center, Hampton, 
Virginia 23365. 



Richmond, J.H., "Computer Program for Thin-Wire 
Structures in a Homogeneous Conducting Medium," NASA 
Contractor Report CR-2399, June 1974, for sale by 
the National Technical Information Service, 
Springfield, Virginia, 22151, Price $3.75. 



Richmond, J.H., "Computer Program for Thin-Wire 
Structures in a Homogeneous Conducting Medium," 
Report 2902-12, August 1973, The Ohio State 
University Electroscience Laboratory, Department of 
Electrical Engineering prepared under Grant NGL 
36-008-138 for National Aeronautics and Space 
Administration, Langley Research Center, Hampton, 
Virginia 23665. 



Richmond, J.H. and Geary, N.H., "Mutal Impedance of 
Nonplanar-Skew Sinusoidal Dipoles," Report 2902-18, 
August 1974, The Ohio State university 
ElectroScience Laboratory, Department of Electrical 
Engineering. 



Miller, E. K. , Poggio.A.J., Burle, G.J., and Selden, 
E.S., "Analysis of Wire Antennas in the Presence of 
a Conducting Half Space: Part I. The Vertical 

Antenna in Free Space," Canadian Journal of Physics, 
50, pp 879-888. 



Miller, E.K., Poggio.A.J., Burle, G.J., and Selden, 
E.S., "Analysis of Wire Antennas in the Presence of 
a Conducting Half Space: Part II. The Horizontal 

Antenna in Free Space, Canadian Journal of Physics, 
50 pp 2614-2627. 



133 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 2 

Cameron Station 
Alexandria VA 22314 



2. Library, Code 0212 

Naval Postgraduate School 
Montery CA 93940 



3. Asst Professor R W Adler, Code 52 5 

Department Of Electrical engineering 
Naval Postgraduate School 
Montery CA 93940 



4. Asst Professor J B Knorr, Code 52 1 

Department Of Electrical Engineering 
Naval Postgraduate School 
Montery CA 93940 



5. Lt Jerry W McCormack 1 

4315 Ashworth Street 
Lakewood CA 90712 



6. Dr. J H Richmond 1 

Electroscience Lab 
Ohio State University 
1320 Kinnear Road 
Columbus OH 43212 



7. J CP McEachen 1 

Naval Sea Systems Command 

Code SEA 06T 
Washington DC 20360 

8. Dr. CM Butler 1 

Electrical Engineering Department 

University of Mississippi 
University MS 38677 

9. Lt. R B Birchfield 1 

Naval Material Command 

Code Mat 034 
Washington DC 20360 

10. Dr. R C Hansen 1 

Suite 218 

17100 Ventura Blvd 
Encio CA 91316 

11. Stewart J Allen 1 

School of Aerospace Med 



Brooks AFB Texas 78235 



134 



12 



1 



. Dr. R C Baird 

National Bureau of Standards 
Electro-magnetics Division 
Boulder COLO 80302 



13. Phillip Blacksmith 1 

Code LZR 

AF Research Lab L G Hanscom FD 
Bedford MASS 01730 

14. Commander OSACEEIA 1 

Attn: ACCC-CED-RP 

(Edwin F Brameli 
Ft Huachuca AZ B5613 

15. Russel M Brown 1 

Code 5252 

Naval Research Laboratory 
Washington DC 20390 

16. John A Downs 1 

Code 3556 

Naval Weapons Center 
China Lake CA 93555 

17. Dennis E Fessenden 1 

Code SA32 

Naval Underwater Systems Center 
New London CONN 06320 

18. Frederick Fine 1 

Code 7944F 

Naval Research Laboratory 
Washington DC 20390 

19. Richard G Fitzserelle 1 

01 ITS 

Boulder COLO 80302 

20. Boaz Gelernter 1 

Attn: AMSEL CT R 

OS Army Electronics Command 
Ft Monmouth NJ 07703 

21. George T Gobeaud 1 

OSA Satellite Communications Agency 

Ft Monmouth NJ 07703 

22. Joseph Halberstein 1 

Code FVN 

Naval Weapons Laboratory 
Dahlgreen VA 22448 

23. Herbert Heffner 1 

Code 20411 

Naval Air Development Center 
Warminster PA 18974 



135 



1 



24. Whilden G Heinard 
Branch 150 Bldq 92 
Harry Diamond Laboratory 
Washington DC 20438 

25. Douglas Hess 1 

USA MERDC 

Ft Belvoir VA 

26. John M Horn 1 

Naval Electronics Laboratory Center 

San Diego CA 9 2152 

27. Egbert H Jackson 1 

Code 423 4 

Naval Ordnance Laboratory 
Silver Springs MD 20910 

28. R H Jones ' 1 

Code 6173 

Naval Ship Engineering Center 
Hyattsville MD 20782 

29. Dr Rudy Kalafus 1 

Transportation Systems Center 

Cambridge MASS 02143 

30. Cyril M Kaloi 1 

Naval Missile Center 

Pt Mugu CA 93041 

31. Abe Kampinsky 1 

Code 110 

NASA Goddard Space Glith Center 
Greenbelt MD 20771 

32. Fan King 1 

Attn: AMSMI REG 

Army Missile Command Redstone Arsenal 
Huntsville ALA 35809 

33. Paul A Lantz 1 

Code 811 

NASA Goddard Space Flight Center 
Greenbelt MD 20771 

34. Maj Anthony Martinez 1 

(LZRj 

Air Force Communications Research Laboratory 
Hanscom Field MASS 01730 

35. Dr E K Miller (L158) 1 

Lawrence Livermore laboratory 

PO Box 808 
Livermore CA 94550 



136 



1 



36. James Muro 

Attn: AMXST SR Z 
FSTC 

220 7th NE 

Carlottesville VA 22901 



37. Martin Natchipolsky 

Federal Aviation Agency 
800 Independence Avenue 
Washington DC 20590 



38. C L Pankiewicz 
OCTS 

Rome Air Develop Center 
Griff iss AFB NY 13441 



39. Sol Perlman 

Attn: AMSEL VL G 
Army Electronic Command 
Ft Monmouth NJ 07703 



40. John Potenza 
(OCTS) 

Rome Air Development Center 
Griffiss AFB NY 13441 



41. Victor W Richard 
Attn: AMXBR CA 

Army Ballistic Research Laboratory 
Aberdeen Proving Ground MD 21005 



42. Emilio Rivera 
Code Air 53356B 
Naval Air System Command 
Washington DC 20360 



43. Dr John W Rockway 

Naval Electronics Laboratory Center 
Code 2120 

San Diego CA 92152 



44. Julius Ross 
Code IBC ET 
OS Information Agency 
Washington DC 20547 



45. F Rouffy 

Army Safeguard System Command 
Attn: SSC HR 
PO Box 1500 
Huntsville ALA 35807 



46. Gaylon E Ryno 
Code 5013 

Naval Weapons Center 
China Lake CA 93555 



1 



1 



1 



1 



1 



1 



1 



1 



1 



1 



137 



1 



47. Dr Felix Schwering 
Attn: AMSEL NL CR 1 
Army Electronics Command 
Ft Monmouth NJ 07703 



48. James R Seale 1 

Naval Air Test Center 

Patuxent River MD 20670 

49. John P Shanklin Jr 1 

Attn: AFAL TEM 

Air Force Avionics Laboratory 
WPAFB OHIO 45433 

50. Carlyle J Sletten 1 

Code LZ 

Air Force Cambridge Research Laboratory 
Bedfored MASS 91730 

51. Joseph M Smiddie 1 

Code 3083 

Naval Ammunition Depot 
Crane IND 47522 

52. A A Strejeck 1 

Code K311 

National Security Agency 
Ft George G Meade MD 20/55 

53. Gottfried Vogt 1 

Attn: AMSEL WL S USAEC0M 

Ft Monmouth NJ 07703 

54. Alvin C Wilson 1 

National Bureau of Standards 

Boulder COLO 80302 

55. Dr Don Dudley 1 

Electrical Engineering Department 

University Of Arizona 
Tuscon AZ 87521 

56. A1 Mink Code 6179 1 

Naval Ships Engineering Center 

Hyattsville MD 20784 

57. Fran Prout Code 11 '74 1 

Naval Ships Engineering Center 

Hyattsville MD 20784 

58. Tony Testa Code 6174 1 

Naval Ships Engineering Center 

Hyattsville MD 20784 

59. Dr S Siahatgar Code 6174 1 

Naval Ships Engineering Center 

Hyattsville KD 20784 



138 



60 



1 



. Dr A Sankar 
MS H-1 1144 

TRW Systems, 1 Space Pack 
Redonao Beach CA 90278 



61. Ronald Prehoda 
Code FVR 

Naval Surface Weapons Center 
Dahlgren VA 22448 



62. Cdr EG Neely III 

Naval Electronics Systems Command 
Code ELEX 094 
Washington DC 20360 



63. Dr R Tanner 

Technology For Communication International 
1625 Stierlin Road 
Mt View CA 94043 



64. Dr Fred Tesche 

Science Applications, Inc 
PO Box 277 
Berkeley CA 94701 



65. Walter Curtis 

Boeing Aerospace Co 
PO Box 3999 
Seatle WA 98124 



66. Dr Raj Mittra 

Electrical Engineering Department 
University Of Illinois 
Drbona IL 61801 



67. Cdr Russ Shields 

Naval Electronics System Command 
Code PME-107 
Washington DC 20360 



68. Donn Campbell 
COMM/ADP 
USAECOM 

Ft Monmouth NJ 07703 



69. Dr B Strait 

Electrical Engineering Department 
111 Link Hall 
Syracuse University 
Syracuse NY 13210 



70. Dr Art Sindoris 
USAECOM 

Ft Monmouth NJ 07703 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



139 



* 








Thesis 156502 



M181S5 McCormack 
c.l "ASAP" Antennas- 

Scatterers Analysis 
Program, a general pur- 
pose user-oriented com- 
puter program for analy- 
sis of thin-wire struc- 
tures in the presence. 




DUDLEY KNOX LIBRARY 




3 2768 00018677 9 



